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Pesquisadores na агеа de controle ie 
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Introdução - 


> Teoria de controle clássico; 
+ Teoria de controle moderno; 
< Teoria de controle robusto; 


*• Controle automático 


+ Engenharia e ciência; 


4$ Sistemas robóticos; 
* Veículos espaciais; 
4$ Manufatura; 


* Operações industriais; 
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Breve história do controle automático a. 


* 1? trabalho de controle automático (James Watt): regulador centrífugo 


para controle de velocidade de uma máquina a vapor; 


EL) 


d TV 7 
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Breve história do controle automático a. 


** Minorsky (1922): controladores automáticos para pilotagem de 
embarcações. Demonstrou como a estabilidade poderia ser determinada a 


partir de equações diferenciais; 
4$ Nyquist (1932): determinou a estabilidade de sistemas de malha 
fechada com base na resposta de malha aberta a excitações senoidais 


estacionárias; 


* Hazen (1934): introduziu o termo servomecanismo para sistemas de 


controle de posição, apresentou o projeto de servomecanismo a relé; 
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и "A" zn ui 
Breve historia do controle automatico EE. 


Nicolas Minorsky Harry Nyquist Graham J. Hazen 
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Breve história do controle automático " 


« Década de 1940: diagrama de Bode; 

1940 e 1950: controladores PID em sistemas de controle industrial; 
4 Início de 1940: regras de sintonia de Ziegler-Nichols; 

“ Final de 1940 (Evans): Lugar Geométrico das Raízes (LGR); 


* LGR e resposta em frequência são a essência da teoria clássica de 


controle; 
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Breve história do controle automático a. 


« Teoria de controle moderno, com base nas variáveis de estado, foi 
desenvolvida para lidar com a crescente complexidade dos sistemas 
modernos e seus rigorosos requisitos relativos à precisão, à importância e 


ao custo em aplicações militares, espaciais e industriais; 


4 1960 e 1980: controle de sistemas determinísticos e estocásticos, 


controle adaptativo, controle de aprendizagem de sistemas complexos; 


* 1980 a 1990: controle robusto e tópicos associados. 
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Breve história do controle automático Ë 
Teoria de controle moderno 


мттитонр 
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4$ Ваѕеіа-ѕе na análise no domínio do tempo em sistemas de equações 
diferenciais. Esta teoria simplifica o projeto de sistemas de controle porque 


se baseia no modelo de um sistema de controle real; 


+ Estabilidade do sistema é sensível ao erro entre o sistema real e seu 


modelo; 


* O controlador é projetado a partir do modelo e pode alcançar um 


sistema instável; 


* O método de projeto para que o sistema náo se torne instável é 


chamada teoria de controle robusto. 
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Breve história do controle automático a. 


* Variável controlada: é a grandeza ou a condição que é medida e controlada; 


% Variável manipulada (sinal de controle): é a grandeza ou condição modificada 


pelo controlador, de modo que afete o valor da variável controlada; 


4 Plantas: pode ser uma parte de equipamento ou apenas um conjunto de 
componentes de um equipamento que funcione de maneira integrada, com o 


objetivo de realizar determinada operação; 


4 Processos: operação natural de progresso contínuo ou um desenvolvimento 
caracterizado por uma série de modificações graduais que se sucedem umas às 


outras de modo relativamente estável, avançando em direção a um dado objetivo; 
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Breve história do controle automático a. 


*$* Sistemas: ё a combinação de componentes que agem em conjunto para atingir 


determinado objetivo. Sistemas físicos, biológicos, econômicos e outros; 


% Distúrbios: é um sinal que tende a afetar de maneira adversa o valor da variável 
de saída de um sistema. Se um distúrbio for gerado dentro de um sistema, ele será 
chamado distúrbio interno, enquanto que o distúrbio externo é aquele gerado 


fora do sistema e que se comporta como um sinal de entrada no sistema; 


«Controle com realimentação: refere-se a uma operação que, na presença de 
distúrbios, tende a diminuir a diferença entre a saída de um sistema e alguma 


entrada de referência e atua com base nessa diferença. 
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Exemplos de sistemas de controle Ë 
Velocidade i 
* A quantidade de combustível fornecida ao motor é ajustada de acordo 


com a diferença entre a velocidade esperada e a velocidade efetiva do 


motor 


Cilindro 
de potencia 


E 
Aceite а 


baja tensión 
| iem 


Motor 


Combustible —= = 


de control 
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Exemplos de sistemas de controle Ë 
Velocidade ui 


* A planta (sistema controlador) é o motor e a variável controlada é a 


velocidade do eixo do motor; 


* A diferença entre a velocidade desejada e a velocidade real é o sinal de 


erro, 


* O sinal de controle (quantidade de combustível) a ser aplicado à planta 


(motor) é o sinal atuante; 


* A grandeza externa que perturba a variável controlada é o distúrbio 


(mudança inesperada na carga). 
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Exemplos de sistemas de controle Ë 
Temperatura i 


Termómetro 


Conversor 


АЛ» 


| 


Controlador 


Horno 
eléctrico 


| DM 
| cador | 


Calefactor 
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Exemplos de sistemas de controle Ë 
Empresariais i 


> É um sistema de malha fechada; 
* О projeto reduzirá o controle administrativo necessário; 


4 Os distúrbios correspondem à carência de mão de obra ou matéria-prima, 


a interrupção de comunicação, a erros humanos entre outros. 


Producte 
requerido 


Diseño de Fabricación | Producto 
producto v y Verificación 
delineación | ensamblaje 
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Controle de malha fechada e malha aberta 


мттитонр 


Malha fechada (controle com realimentacao 


** O sinal do erro atuante, realimenta o controlador, de modo a minimizar o erro 


e acertar a saída do sistema ao valor desejado; 


+ O termo controle de malha fechada sempre implica a utilização do controle 


com realimentação para reduzir o erro do sistema. 


Distúrbio 


+ | Variável de 


Variável de 
Controle 


Variável de 
Referência 


Elemento 
Sensor 
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Controle de malha fechada e malha aberta É 
Malha aberta ulli 


5$ Os sistemas de malha aberta são aqueles que o sinal de saída não exerce 
nenhuma ação de controle no sistema; 

“ O sinal de saída não é medido nem realimentado para a comparação com a 
entrada; 

* Máquina de lavar roupas: as operações são executadas em uma sequência 
baseada no tempo. A lavadora não mede o sinal de saída, isto é, não verifica se 
as roupas estão bem lavadas; 


“ Controle de tráfego: operado em função do tempo. 


Distúrbio 


Variável de 
Saida 


Variável de 
Controle 
Controlador Processo 


Variável de 
Referência Erro 
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Controle de malha fechada e malha aberta Ë 
Malha fechada x Malha aberta 


Malha aberta Malha fechada 


Sensível a distúrbios externos е а | Relativamente insensível a distúrbios 
variações internas do sistema. externos e a variações internas do 
sistema (devido a realimentação). 


Em relação a estabilidade, é mais fácil | Estabilidade её um problema 
de ser construído, pois a estabilidade é | importante. 


um problema menos significativo. 


Sistemas onde as entradas são | Vantajosos nos casos em que houver 

conhecidas com antecipação e são |distúrbios e/ou alterações пао 

isentos de distúrbios. previsíveis nos componentes do 
sistema. 


Menor número de componentes | Maior número de componentes 
utilizados. utilizados. 
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Controle de malha fechada e malha aberta É 
Malha aberta (vantagens) dl 


* São simples de serem construídos e possuem fácil manutenção; 
* São menos dispendiosos que um sistema correspondente de malha fechada; 
* Não apresentam problemas de estabilidade; 


* São adequados quando existem dificuldades de medição da saída ou quando 


a medição precisa da saída não é economicamente possível; 
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Controle de malha fechada e malha aberta Ë 
Malha aberta (desvantagens) ili 


% Distúrbios e mudanças na calibração causam erros, e a saída pode 


apresentar diferenças em relação ao padrão desejado; 


4 Para que a saída mantenha a qualidade requerida, é necessária uma 


regulagem periódica; 
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Projeto e compensacao de sistemas Ë 
de controle 


* Compensação: é a modificação da dinâmica do sistema para satisfazer 


as especificações dada; 
* Em projetos reais de sistemas de controle, a utilização de um 
compensador eletrônico, pneumático ou hidráulico, é uma questão que deve 


ser decidida em parte na natureza da planta a ser controlada; 


% Se a planta a ser controlada inclui líquido inflamável, é necessário 


escolher componentes pneumáticos; 
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Projeto e compensação de sistemas de controle Ë 
Especificações de desempenho mes 


% Requisitos importantes em um sistema de controle; 
4$ As especificações podem ser em termos de requisitos de resposta 
transitória e de regime estacionário ou em termos de resposta em 


frequéncia; 


* As especificacóes de desempenho nào devem ser mais restritivas que o 


necessário para a realização da tarefa em questão; 
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Projeto e compensacao de sistemas Ë 
de controle 
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* Compensação do sistema: 
“ ajustar o ganho é o primeiro passo no ajuste do sistema para um desempenho 
satisfatório; 
% Reprojeto ou acréscimo de um dispositivo adequado; 
4$» Compensador: dispositivo inserido ao sistema com o propósito de satisfazer às 


especificações de projeto; 


“ Procedimento de projeto: desenvolve-se um modelo matemático do sistema de 
controle e ajusta-se os parâmetros de um compensador; 
* Verificação do desempenho do sistema por meio da análise de cada ajuste dos 
parâmetros; 
< Software de simulação; 
+ Implementação real = simulação: deve-se ajustar. 
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F- AGENDA 


+ Análise por Espaço de Estados 
+ Estabilidade de Sistemas Dinâmicos 


% Critério da estabilidade de Routh 


% Erros estacionários em sistemas de controle 
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E. BIBLIOGRAFIA 
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- ANÁLISE POR ESPACO DE ESTADOS 


* Considere a equação de estado 


x = Ах+Ви 


* е a equação de saida 


y = (x + Du 


* Aplicando-se a Transformada de Laplace, tem-se: 


sX(s) — x(0) = AX(s) + BU(s) Ө 
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F. ANÁLISE POR ESPAÇO DE ESTADOS 


*• |solando-se X(s) (equação de estado), tem-se: 


(sI — A)X(s) = x(0) + BU(s) 
X(s) = (sI — AJTIx(0) + (sI — A) 1 BU(s) 


. adj(sl — A) 


"ue det(sl — A) 


[x(0) + BU(s)] 
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- ANÁLISE POR ESPACO DE ESTADOS 


* Os polos da função de transferência determinam a natureza de 


resposta transiente do sistema. 


** Os polos de um sistema representado em Espaço de Estado são iguais 


aos autovalores da matriz A. 


** Os polos são obtidos a partir da expressão det(SI - 4) = 0. 


Y(s) _ .|adj(sl — А) 
U(s) | in — A) lids 
Y(s) Cadj(sI — A)B + Ddet(sI — A) Ө 


U(s) | det(sI — A) 
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F. ANÁLISE POR ESPAÇO DE ESTADOS 


* Exemplo 1: Dado o sistema representado em EE, faça: 
* Resolva a equação de estado precedente e obtenha a saída para a 
entrada exponencial fornecida. 


> Obtenha os autovalores e os pólos do sistema. 


0 1 0 0 
x=|0 0 1]х+|0|е* 
—24 —26 —9 1 

у= |1 1 olx 


q» Ө 
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- ANÁLISE POR ESPAÇO DE ESTADOS 


* Exemplo 2: Dado o sistema representado em EE, faça: 
“ Resolva para y(t) utilizando as técnicas de EE e da TL. 


* Obtenha os autovalores e os polos do sistema. 
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- ANÁLISE POR ESPACO DE ESTADOS 


* Exemplo 3: Dado o sistema representado em EE, faça: 
* Resolva para y(t) utilizando as técnicas de EE e da TL. 
> Obtenha os autovalores e os polos do sistema. 


* Matriz de transição de estados. 


Me 
| 
N 
E 
У 
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- ESTABILIDADE DE SISTEMAS DINÂMICOS 


Estabilidade: é especificacáo mais importante do sistema; 


Se um sistema é instável, a resposta transiente e os erros do estado 


estacionário deixam de ter significado; 


Um sistema instável nào pode ser projetado para atender uma resposta 


transiente específica ou um requisito de erro de estado estacionário; 
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F. ESTABILIDADE DE SISTEMAS DINÂMICOS 


« Sistema Estável: se a resposta natural tende a zero quando o tempo 


tende a infinito (sistema linear invariante no tempo). 


« Sistema Instável: se a resposta natural aumenta sem limites, na 
medida que o tempo tende a infinito (sistema linear invariante no 


tempo). 


“ Sistema Marginalmente Estável: quando a resposta natural não 


apresenta aumento, nem atenuação, mas permanece constante ou 


oscila quando o tempo tende a infinito (sistema linear rvariant m 
tempo). 
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- ESTABILIDADE DE SISTEMAS DINÂMICOS 


* A definição de estabilidade implica a permanência da resposta forçada 


quando a resposta natural tende a zero. 


% Estabilidade BIBO (Bounded-Input, Bounded-Output): 
* Um sistema é estável se toda entrada limitada gerar uma saída 
limitada. 
% Um sistema é instável se alguma entrada limitada gerar uma 


saída ilimitada. 
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- ESTABILIDADE DE SISTEMAS DINÂMICOS 


“+ Como determinar se um sistema é estável? 


% Os sistemas estáveis possuem função de transferência de malha 
fechada com polos somente no semiplano esquerdo (Respostas 


naturais exponenciais decrescentes ou senóides amortecidas). 


* Os sistemas instáveis possuem função de transferência de malha 
fechada com pelo menos um polo no semiplano s da direita e/ou 


polos de multiplicidade maior que um no eixo imaginário. 
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F. ESTABILIDADE DE SISTEMAS DINÂMICOS 


4 Como determinar se um sistema é estável? 


% Os sistemas marginalmente estáveis apresentam função de 
transferência de malha fechada somente com polos de 
multiplicidade 1 no eixo imaginário e polos no semiplano s da 


esquerda (resposta natural, oscilações senoidais puras); 
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= ESTABILIDADE DE SISTEMAS DINÁMICOS 


R(s) = 1/s 3 Y(s) 


s(s + 1)(s + 2) 


Step Response 
1.6 Root Locus 
6 T T T T 
1.4 | / 
4 System: Y 
| Gain: 0.0558 7 
1.2 ^ System Y Pole: -0.153 + 1.07i 
Gain: 0 Damping: 0.141 
2. Pole: -2.67 Overshoot (%): 63.9 || 
1 Damping: 1 Frequency (rad/sec): 1.08 
T; E Overshoot (96): 0 = 
5 & Frequency (rad/sec): 2.67 
> = 
E 0.8 5 0 " 
E 5 
[e] 
Ë 
0.6 = 
о. 
0.4 
-AL 
0.2 
-6 l l l 
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0 5 10 15 20 25 30 35 40 —€— 
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CRITÉRIO DE ROUTH-HURWITZ 


* Por este método, pode-se dizer quantos polos do sistema em malha 
fechada estão no semiplano esquerdo, no semiplano direito e sobre o 


eixo jo. 
% Não precisa-se calcular os polos do sistema em malha fechada. 


* Etapas: 
< Gerar uma tabela de dados denominada tabela de Routh. 


€ Interpretar a tabela de Routh para definir quantos polos do 


sistema em malha fechada se situam no semiplano esquerdo 
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semiplano direito e sobre o eixo jo. 


m _ CRITÉRIO DE ROUTH-HURWITZ 
ETAPA 1: TABELA DE ROUTH 


* Observe a função de transferência de malha fechada. 


К(5) N(s) Ү (5) 


dust + ass3 + ass? + ast + ags? 


% Gera-se a Tabela de Routh nomeando as linhas com potência de s а 


partir da potência mais alta no denominador da função de 
transferência de malha fechada até sº. 


Engenharia de Controle | - Rejane de Barros Araújo 


m _ CRITÉRIO DE ROUTH-HURWITZ 
ETAPA 1: TABELA DE ROUTH 


« Entradas Remanescentes: cada entrada é igual ao valor negativo 


dos determinantes formados com os elementos das duas linhas 


anteriores, dividido pelo elemento da 1° coluna diretamente "Ө 
linha que está sendo calculada; 
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m _ CRITÉRIO DE ROUTH-HURWITZ 
ETAPA 1: TABELA DE ROUTH 


* A coluna à esquerda do determinante é sempre а 1º coluna das duas 
linhas anteriores, e a coluna à direita é constituída dos elementos da 


coluna a cima e à direita. 
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- CRITÉRIO DE ROUTH-HURWITZ 
ETAPA 1: TABELA DE ROUTH 


E = | Qo ju 0 
me аз 04 LE аз 0 bs = аз 0 0 
аз аз аз 

" "X as 0 jaz 0 
lb, b B b O 
c, = 1 2 р; 0 1 
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- CRITÉRIO DE RoUTH-HURWITZ 
ETAPA 1: TABELA DE ROUTH 


db b; “lb o dbi | 
c 0 c4 O "M cı 0|_ 
1 C 2 C C, 
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CRITÉRIO DE ROUTH-HURWITZ 
ETAPA 2: INTERPRETAÇÃO DA TABELA DE ROUTH 


мттлонр 


Os sistemas seráo interpretados com polos reais nos semiplanos 


esquerdo e direito. 


O critério de Routh-Hurwitz estabelece que o número de raízes do 
polinômio que se situam no semiplano direito é igual ao número de 


mudanças de sinal na 1º coluna. 


Se a função de transferência de malha fechada possui todos os polos 


na metade esquerda do plano s, o sistema é estável; 


Um sistema é estável se não ocorrerem mudanças de sinal 9 
coluna da tabela de Routh; 
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E CRITÉRIO DE ROUTH-HURWITZ 


* Exemplo 4: Construa a tabela de Routh e interprete para o sistema da 
1000 
(s + 2)(s + 3)(s + 5) 


% Exemplo 5: Construa uma tabela de Routh e determine quantas 


figura. 


raizes do polinômio a seguir estão no semiplano da direita e no 


semiplano da esquerda; 
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m _ CRITÉRIO DE ROUTH-HURWITZ 
CASOS ESPECIAIS 


* Zeros apenas na 1º coluna: 


% Se o 1º elemento de uma linha é igual a zero, seria necessário 


uma divisão por zero para formar a próxima linha. 


* Para evitar a divisão por zero, atribui-se um valor épsilon, e, em 


substituicáo ao zero na 1? coluna. 


* Faz-se o épsilon tender a zero por valores positivos ou negativos, 


após, os sinais dos elementos na 1° coluna, podem ser 


determinados; o 


Engenharia de Controle 1— Rejane de Barros Araüjo 


m _ CRITÉRIO DE ROUTH-HURWITZ 
CASOS ESPECIAIS 


* Zeros apenas na 1º coluna: 


% Outro método que pode ser utilizado é substituir o polinômio 


procurado, pelo polinômio original com seus coeficientes escritos 


na ordem inversa; 


“ Exemplo 6: Determine a estabilidade da função de transferência de 


malha fechada 


10 
5° + 25* + 353 + 652 + 55 + 3 


0 
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T(s) = 


m _ CRITÉRIO DE ROUTH-HURWITZ 
CASOS ESPECIAIS 


% Exemplo 6: Solução 


i 

кише ТГ 

"T ma. 

И ЭМИ ЖШ 
S 


(42e — 49 — 6є2)/(12є — 14) 


s| 083 oO jo 
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- CRITÉRIO DE ROUTH-HURWITZ 
CASOS ESPECIAIS 


* Exemplo 6: Solução 
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m _ CRITÉRIO DE ROUTH-HURWITZ 
CASOS ESPECIAIS 


* Exemplo 6: Solução 
% Sendo €*, ocorrerá uma mudança de sinal da linha sz para a linha 
S2, e haverá uma mudança da linha s; para a linha s, > sistema 


instável com dois polos no semiplano da direita. 


% Sendo €”, haverá uma mudança de sinal da linha s, para a linha 
s, € uma mudança da linha sz para a linha s; Ə sistema instável 


com dois polos no semiplano da direita; 
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m _ CRITÉRIO DE ROUTH-HURWITZ 
CASOS ESPECIAIS 


* Exemplo 7: Usando a função de transferência de malha fechada do 
exemplo 3, determine a estabilidade do sistema para o método da 


inversão dos coeficientes do polinômio característico; 


* Linha Completa de Zeros: 


* A linha completa de zeros na tabela de Routh é função da 
existência de um polinômio estritamente par ou estritamente 


impar, como fator do polinômio original. 
** Um polinômio par possui apenas potências pares de s; 


* Os polinômios pares possuem apenas raízes simétricas é 


relação a origem; 
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m _ CRITÉRIO DE ROUTH-HURWITZ 
CASOS ESPECIAIS 


% Linha Completa de Zeros: 


4 Е o polinômio par o responsável pelo aparecimento da linha de 


Zeros. 


*% Alinha anterior a linha de zeros, contém o polinômio par. 


Ф, 


* А simetria pode ocorrer em três condições de posição das raízes: 
* As raízes são simétricas e reais. 
* As raizes são simétricas e imaginárias. 


* As raízes são quadrantais. Ө 
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CRITÉRIO DE ROUTH-HURWITZ 
CASOS ESPECIAIS 


a 


emgernmuna uv curvo а Jane de Barros Araújo 


E CRITÉRIO DE ROUTH-HURWITZ 


* Exemplo 8: Determine o numero de polos no semiplano da direita 


referente à funcáo de transferéncia. 


T(s) 10 

^] =Z ——— É—— ES 
S? + 754 + 6s3 + 4252 + 8s + 56 

* Exemplo 9: Determine o número de polos no semiplano da direita, no 
semiplano da esquerda e sobre o eixo imaginário da função de 


transferência. 


20 


T = сс. 
(5) s8 + s? +1256 + 22s? + 395“ + 59s? + 4852 + 385 + 20 Ө 
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E CRITÉRIO DE ROUTH-HURWITZ 


* Exemplo 9: Solução 


1/3 


= й 2 dP (s 
Р(5) = s^ + 35 + 2 XO) — Ast + 6s+0 
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üjo 


E. CRITÉRIO DE ROUTH-HURWITZ 


% Exemplo 9: Solução 
* Mudança de sinal da linha s^ para a linha sê e outra mudança da 


linha sé para s? Ə dois polos no semiplano direito (raízes 


simétricas). 
% Dois polos no semiplano esquerdo (raízes simétricas). 


% Quatro polos sobre o eixo jo (raízes quadrantais). 


* Polinómio par. D 
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E CRITÉRIO DE ROUTH-HURWITZ qe 


* Exemplo 10: Determine o número de polos no semiplano da direita, 


no semiplano da esquerda e sobre o eixo imaginário da função de 


200 
s(s? + 6s? + 11s + 6) 


* Exemplo 11: Determine o número de polos no semiplano da direita, 


transferéncia. 


R(S) 


no semiplano da esquerda e sobre o eixo imaginário da função de 


transferéncia. 1 


Т(5) = 


25° -3s* + 253 + 352 + 25 + 1 
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E CRITÉRIO DE ROUTH-HURWITZ qe 


* Exemplo 12: Determine a faixa de valores de ganho K, para o 
sistema mostrado na figura, que fará com que o sistema seja estável, 


instável e marginalmente estável. Admita K > 0. 


* Exemplo 13: Determine a faixa de valores de ganho K, para o 


sistema mostrado na figura, que fará com que o sistema seja estável. 


K(s + 20) 


D s(s + 2)(s + 3) ü 
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E CRITÉRIO DE ROUTH-HURWITZ je 


* Exemplo 14: Dados os sistemas em Espaço de Estado. Determine a 


quantidade de polos no semiplano da esquerda, no semiplano da direita 


e sobre o eixo jo. 


0 3 1 10 

х=| 2 8 1I|x+l0Olu y=|1 0 olx 
—10 —5 -2 0 
2 1 1 0 

х=|1 7 1|x>+)0|r у= [0 1 olx 
—3 4 —5 1 © 
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ERROS ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 


* Alterações na entrada de referência causaráo erros inevitáveis durante о 


regime transitório, podendo causar também erros estacionários. 


* Qualquer sistema de controle físico apresenta, inerentemente, erros 


estacionários na resposta a certos tipos de entradas. 


* O erro estacionário que um sistema apresenta em relação a determinado 


tipo de entrada depende do tipo de função de transferência de malha 


aberta desse sistema. 
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ERROS ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 


* Considere o sistema de controle com realimentação unitária, com a 
seguinte função de transferência de malha aberta: 


K(T4s + 1)(Tys + 1) -*(Tins + 1) 
(5) == 
sN(T,s + D(Tos + 1) =- (Tps + 1) 


* s" polo de multiplicidade N. 


* 


« Número de integrações indicadas pela função de transferência de malha 


aberta. 


€ Sistemas tipo 0, 1, 2,...., ѕе N = 0, N = 1,.... 
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ERROS ESTACIONÁRIOS EM SISTEMAS ОЕ CONTROLE 


R(S) + š K С(5) 


C(s) G(s) a C(s) | 1 
R(s) 14G(s) R(s) ^ R(s 1-4G(s) 
1 
E (s) j r 56) 

sR(s) 


= = lim e(t) = lim s E(s) = = lim 


$201 + G(s) 
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ERROS ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 


Physical Time Laplace 
Waveform Name interpretation function transform 
кіт) Step Constant position I 1 
Kp 5 
I 
# Ramp Constant velocity t > 
Kv 
| 
rit) 1 1 
Parabola Constant acceleration =P -= 
Ka 2 55 


" e 
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Ёввоз ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE POSIÇÃO - K, 


« Erro estacionário do sistema para uma entrada em degrau 


=p $ 1_ 1 
“ss S01+G(s)s 1+6(0) 


* A constante de erro estático de posição é definida por: 
K, = lim G (s) = G(0) 
S0 


“+ Então o erro estacionário é: 
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tERROos ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE POSIÇÃO - K, 


* Рага um sistema do Tipo 0 


"C K(T4s + 1)(Tys + 1) -- X 1 
= = = е = — 
P^ 550 (Tis 1)(Ту& + 1) --. — iut. 


* Para sistemas do Tipo 1 ou maior 


т K(Tas + 1) (75 +1) -- 


— me = o0 е = 0 
P' soos (T,s + 1)(T>s + 1) ~ os 
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tERROos ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE VELOCIDADE - K, 


* Erro estacionário do sistema para uma entrada em rampa 
I S 1 I 1 1 
е. = llm ———- = 11M ——— = —— TT 
% s>014+G(s)s? s>0s+G(s) lim s G(s) 
S 


* A constante de erro estático de velocidade é definida por: 
K, = lims G(s) 
S—0 
* Então o erro estacionário é: 


a 
E © 
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tERROos ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE VELOCIDADE - K, 


* Para um sistema do Tipo 0 


sK(T4s + 1)(T,s + 1) -- 1 


K, = lim =0 — e. = — = о 


S—0 (Т, 5 + 1)(T5s T 1) “т 
* Para sistemas do Tipo 1 


. SK(T,s + 1)(Tys + 1)-- 
K, = li 


DOC. > l= 
о s(T,s + 1)(T;js +1) ·- с 
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tERROos ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE VELOCIDADE - K, 


* Para um sistema do Tipo 1 


r(t) 
c(t) 


c(t) 


0 t 


* Para sistemas do Tipo 2 ou maior 


ms 1 
SK(T4s + 1)(Tys + 1) m" pousse =й 


” Ө 
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E o sN(T,s + 1)(T5s + 1) т 


tERROos ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE ACELERAÇÃO - K, 


* Erro estacionário do sistema para uma entrada em parábola 


2] S 1 — 1 
"ss = S301 G(s)s ` lim s2 G (s) 
S> 


* A constante de erro estático de aceleração é definida por: 
— 1; 2 
Ka = lim s“ G(s) 
S—0 
* Então o erro estacionário é: 


1 
= @ 
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tERROos ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE АСЕГЕВАСАО - K, 


* Para um sistema do Tipo 0 


. s^K(T,s + 1)(T,s + 1) -- 1 
Ка = lim—————M————————-—0 — ls, = 00 
590 (Ts + 1)(Tjs +1) Ка 
* Рага sistemas do Tipo 1 
s^K(T,s + 1)(T,s + 1)-- 1 
Ka = lim (Ta XT, ) = 0 —— ĉl, = = Oo 


| $0 s(T4s + 1)(T5s + 1) ::: Ka 
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tERROos ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE АСЕГЕВАСАО - K, 


* Para um sistema do Tipo 2 


— 5 К(Таз + 1)(Ть5 +1) _„ o 1 1 
a $50 s2(Tis+1)(T,s+1)-- — " Ka 


rt) 
c(t) 


r(t) 
e(t) 
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"ERROS ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 
CONSTANTE DE АСЕГЕВАСАО - K, 


* Para um sistema do Tipo 3 ou maior 


"- s^K(T4s + D)(Tps + 1) -- 1 А 
= ————— — OO = — E= — 
а $50 s"(T,s + 1)(T>s + 1) --- — “fs 


Entrada em Entrada em Entrada em 
degrau rampa асе!егасао 
r(t) = 1 r(t) = t r(t) = (1/2)t2 
Sistema do 1 
tipo 0 1+K 


Sistema do 1 
tipo 1 K 

Sistema do 1 
tipo 2 K 
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ERROS ESTACIONÁRIOS EM SISTEMAS ОЕ CONTROLE 


* Exemplo 15: Determine os erros de regime permanente para as 
entradas 5u(t), 5tu(t) e 5t“u(t) no sistema mostrado na figura. A função 


u(t) é um degrau unitário. 


R(S) K) E(S) 


100(s + 2)(s + 6) 
s(s + 3)(s + 4) 


G(S) = 
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ERROS ESTACIONÁRIOS EM SISTEMAS DE CONTROLE 


* Exemplo 16: Um sistema com realimentação unitária possui a seguinte 


função de transferência no caminho direto: 


o 1000(s +8) 
(5- GY5G-9 


* Determine o tipo de sistema Kp, K, e Ka. 


* Utilize as repostas anteriores para determinar os erros de regime 


permanente para os sinais de entrada padronizados: degrau, rampa; 


parabólico. 
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ENGENHARIA DE CONTROLE I 


UNIDADE Ill: CONTROLADORE PID 


Engenharia de Controle e Automação 
Prof.? Dr.? Rejane de Barros Araújo 


- SISTEMA DE CONTROLE EM MALHA FECHADA 


GJ . | 1] 1 OD. 


“ Controladores geralmente utilizados na indústria: 
m Controle proporcional: controle com ganho constante 
m Controle proporcional-integral: melhora o erro em regime permanente 
m Controle proporcional-derivativo: melhora as respostas transitórias 


m Controle proporcional-integral-derivativo: tenta combinar o melhor dos tipos de 


y(t) 
| 


controladores anteriores 


t) 


sinal de referencia 


tealimentacao 
controle em malha fechada 
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Ё. AÇÕES BÁSICAS DE CONTROLE 


% Controlador: 
m Dispositivo que realiza um conjunto de operações matemáticas sobre 


um sinal de erro e(t). 


» Objetivo: 
m Produzir um sinal de controle u(t) para ser aplicado no processo; 
» Tenta atender determinados objetivos do sistema operando em 


malha fechada. 
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Е. _ AÇÕES BÁSICAS DE CONTROLE 


7 AcAo “BANG-BANG” 


% O controlador pode assumir apenas 2 valores, conforme o erro seja 


positivo ou negativo: 


U, se e(t) > 0 


nct) = 
U, se e(t) <D 


4 OBS.: Caso um dos sinais (U, ou U5) sejam nulos, a ação “bang-bang” é 


chamada de ação liga-desliga. 
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E AÇÕES BÁSICAS DE CONTROLE 


mb АСАО “BANG-BANG” 


% Para evitar que ruídos em torno do erro e(t) possam levar a 


chaveamentos desnecessários, pode-se utilizar um controle “bang-bang' 


com histerese: u(t) 


El e(t) 


% Se u(t) = U,, é necessário que o valor de e(t) desça abaixo de - E,(t) 
para que haja um chaveamento para U.; 
4$ Se u(t) = U}, é necessário que o valor de e(t) ultrapasse o valor de E, (t) 


para que haja um chaveamento para U}. 
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- AÇÕES BÁSICAS DE CONTROLE 


END AcAo “BANG-BANG” 
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E AÇÕES BÁSICAS DE CONTROLE 
E AcAo “BANG-BANG” 


*• Pode ser considerada a ação mais simples e económica existente; 


* Possui muitas limitações no que diz respeito ao comportamento do 


sistema; 


% Aplicação restrita a sistemas onde não é necessário precisão ou um bom 


comportamento dinâmico. 
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E AÇÕES BÁSICAS DE CONTROLE 


RED AcAo “BANG-BANG” 


* Exemplo: Termostato de geladeira e ar-condicionado, Controle de 


nível usando boias 


— 


| 

! 

| , Flutuador 
| I 


⁄ 


R 
Sistema de controle de nível de líquido com controlador liga-desliga. 
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- AÇÕES BÁSICAS DE CONTROLE 
AÇÃO PROPORCIONAL 


% O sinal de controle aplicado à planta é proporcional à amplitude 


do erro. 
r(t) + e(t) 
u(t) = K¿elt) o 
y(t) 


y(t) < r(t) > e(t) > 0= u(t) > 0 
y(t) > r(t) > e(t) «x0 u(t) < 0 
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E AÇÕES BÁSICAS DE CONTROLE 


END ACAO PROPORCIONAL 


“Tensão de 
- Saida 
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E AÇÕES BÁSICAS DE CONTROLE 


EP ACAO PROPORCIONAL 


“ Para uma planta do tipo O (sem integradores na malha direta): 
* Quanto maior K. > menor o erro em regime. 
% O aumento da ação proporcional melhora a precisão do sistema em 


malha fechada, embora nào consiga zerar esse erro. 


* Exemplo 1: calculando a função de transferência do sistema em malha 


20(s + 5) 
(s + 2)(s + 10)(s + 15) 
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fechada. 


r(t) + | 


- AÇÕES BÁSICAS DE CONTROLE 


БЕ» AÇÃO PROPORCIONAL 


* Exemplo 1: continuação 


4$» Supondo R(s) uma entrada degrau e aplicando o teorema do valor 


final: 
y(oo) = lim sY (s) = limsT(s)R(s) 
s>0 s>0 


| 20K,(s + 5) 1 100K, 
y (00) = lim | s———— —— r skJ— | = ———— 
sol (s + 2)(з + 10)(s + 15) + 20K,(s +5)s) 300 + 100K, 
T 100K, 
оо) = 1 — ————— 
Е 300 + 100K, 


К. > о = e(0)50 f 
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E AÇÕES BÁSICAS DE CONTROLE 


ER ACAO PROPORCIONAL 


* Exemplo 1: continuação (Matlab) 
4 n=20*[1 5]; 
* а = conv(conv([1 2],[1 10]),[1 15]); 
4$ G = tf(n,d) 
* Yo Kc = 1 10 100 1000 (escolha de 4 valores de Kc) 
« Gmf = feedback(Kc*G,1) 
$* step(Gmf) 


** Avalie o transitório e o erro de regime permanente. 
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- AÇÕES BÁSICAS DE CONTROLE 


END ACAO PROPORCIONAL 


% Conclusões: 
* Ao aumentar-se o ganho proporcional, diminui-se o erro em regime 


permanente, porém sem torná-lo nulo; 
« Um aumento indiscriminado no valor do ganho proporcional quase 


sempre torna o sistema muito oscilatório ou pode até mesmo levar o 


sistema à instabilidade. 
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E AÇÕES BÁSICAS DE CONTROLE 


END ACAO INTEGRAL 


% Consiste em aplicar um sinal de controle u(t) proporcional à integral do 


sinal de erro e(t): 


u(t) = E | e(t)dt 


% T; é chamado de tempo integral ou “reset time"; 
“ А ação integral tem uma função de armazenamento de energia; 


« Quando o erro se anula (e(t) = 0), o sinal de controle u(t) será mantido 
em um valor constante, que é proporcional à energia armazenada até o 
instante t. 
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E AÇÕES BÁSICAS DE CONTROLE 


END ACAO INTEGRAL 
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Е. _ AÇÕES BÁSICAS DE CONTROLE 


ED ACAO INTEGRAL 


* O controlador integral pode ser representado por: 


“ Vantagens: 
“ Aumenta em 1 o tipo do sistema, permitindo obter erro nulo para uma 
entrada constante; 


* Melhora a precisão do sistema de regime permanente. 


“ Desvantagens: 
* O polo na origem tende a tornar a resposta do sistema mais lenta; 


+ Pode vir a tornar o sistema instável. 
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E AÇÕES BÁSICAS DE CONTROLE 


END ACAO DERIVATIVA 


€ Corresponde a aplicação de um sinal de controle proporcional à 


derivada do sinal de erro: 
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E AÇÕES BÁSICAS DE CONTROLE 
AÇÃO DERIVATIVA 


A 


mu nm 
ra | É | 


? a Het vo г 
зжалі I 


dVi Ps 
/o (1) C di 
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- AÇÕES BÁSICAS DE CONTROLE 


ED ACAO DERIVATIVA 


% A implementação física de um derivador puro пао é possível na 


prática; 


4$ Para ganhos elevados o sistema seria extremamente sensível а 


ruídos de alta frequéncia; 


% Na prática inclui-se um polo juntamente com a ação derivativa 


para limitar o ganho em alta frequência. 
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E AÇÕES BÁSICAS DE CONTROLE 


NED ACAO DERIVATIVA 


% Ação derivativa com polo extra 


Dilerenciador puro Dilerenciador com polo 
40 40 
35 
зо 
зо 
20 
25 
3 10 3 20 
[v] © 
15 
o 
10 
10 
5 
' o i z o o i 2 3 
107 10 10 10 10 10 10 10 
wírad/s) w(rad's) 


limitacáo do ganho em alta freqüëncia 
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- AÇÕES BÁSICAS DE CONTROLE 


22 AÇÃO DERIVATIVA 


Ф 


% Aplicar um sinal de controle derivativo é o equivalente a aplicar uma 


ação de controle baseada na tendência de evolução do erro; 


LÁ 


* A ação derivativa é antecipatória ou preditiva, pois baseia-se na 


«t. 


tendência de evolução do erro; 


«t. 


* Faz com que o sistema reaja mais rápido, melhorando a resposta 


transitória do sistema; 


«t. 


* Quando o sinal de erro for constante, a ação derivativa será nula, ou 


seja, nào tem qualquer efeito no regime permanente. 
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E CONTROLADORES INDUSTRIAIS 


ШЕЕ CONTROLADOR PROPORCIONAL (P) 


* é essencialmente um amplificador com ganho ajustável; 


* a relação entre o sinal de saída do controlador u(t) e o sinal de erro 


atuante e(t) é dada pela relação: 


an 


Diagrama de blocos do controlador proporcional 
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E CONTROLADORES INDUSTRIAIS 


= CONTROLADOR PROPORCIONAL (P) 


A Ay 
Eai О saq; 
R 
HN, | Re || 
R ie » | + | 
ev( + | 0 
| | Ms | 
B Tensão de 
mha ЭР ! IR ору 

( : | | [S v( Saida 
da aia o 


Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


E. CONTROLADORES INDUSTRIAIS 


NED CONTROLADOR PROPORCIONAL (P) 


% Efeito da variação do ganho proporcional 


Step Response 


Amplitude 


Time (seconds) 
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- CONTROLADORES INDUSTRIAIS 


EDD CONTROLADOR PROPORCIONAL-INTEGRAL (PI) 


* А ação de controle de um controlador proporcional-integral é definida 


por: 


Ї 


u(t) = K e(t) + А | e(t)dt 


г 0 


% O principal objetivo da ação integral é eliminar о erro em regime 


permanente para uma entrada (referência) constante. 


+ É necessário ter cuidado para não tornar o sistema muito oscilatório ou 


até mesmo instável ao se ajustar a ação integral. 
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E CONTROLADORES INDUSTRIAIS 


EPA CONTROLADOR PROPORCIONAL-INTEGRAL (PI) 


% Efeito da variação do tempo integral. 


Step Response 


о 
со 


0.6 


Amplitude 
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E CONTROLADORES INDUSTRIAIS 


RE» CONTROLADOR PROPORCIONAL-INTEGRAL (PI) 
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E CONTROLADORES INDUSTRIAIS 


ED CONTROLADOR PROPORCIONAL-DERIVATIVO (PD) 


% O efeito da ação de controle proporcional-derivativa pode ser expresso 


pela relação: 


de(t) 


u(t) = K e(t)+ K T, 


« A ação derivativa tem a função de antecipar a ação de controle para que 
O processo reaja mais rápido, evitando o aparecimento de transitórios 


com característica oscilatória. 


4$ Pode-se dizer que o sinal de controle a ser aplicado é proporcional а 
uma predição da saída do processo. 
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E CONTROLADORES INDUSTRIAIS 


RE» CONTROLADOR PROPORCIONAL-DERIVATIVO (PD) 
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- CONTROLADORES INDUSTRIAIS 


NED CONTROLADOR PROPORCIONAL-DERIVATIVO (PD) 


% Efeito da variação do tempo derivativo. 


Step Response 


Amplitude 


Time (seconds) 
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- CONTROLADORES INDUSTRIAIS 


RE» CONTROLADOR PROPORCIONAL-INTEGRAL-DERIVATIVO (PID) 


4 Combina as vantagens dos controladores proporcional, proporcional- 
integral e proporcional-derivativo, possuindo as vantagens de cada uma 


das trés acóes individuais: 
* a ação proporcional está ligada à rapidez da resposta do sistema; 


% a ação integral está relacionada à precisão do sistema (garante erro nulo 


em regime permanente para uma referéncia constante); 


* a ação derivativa tende a aumentar a estabilidade relativa do sistema 
(torna a resposta do sistema mais rápida devido ao seu efeito 
antecipatório). 
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- CONTROLADORES INDUSTRIAIS 


BEES CONTROLADOR PROPORCIONAL-INTEGRAL-DERIVATIVO (PID) 
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- CONTROLADORES INDUSTRIAIS 


ШЕ; CONTROLADOR PROPORCIONAL-INTEGRAL-DERIVATIVO (PID) 


€ A equação de um controlador PID é dada por: 


u(t) = K e(t) + T | edt +K.L, di 


i 0 
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- CONTROLADORES INDUSTRIAIS 


ET? CONTROLADOR PROPORCIONAL-INTEGRAL-DERIVATIVO (PID) 


% О controlador PID é a estratégia de controle mais utilizada na 
indústria: 
* consegue afetar a rapidez do sistema, além de seus períodos 


transitórios e em regime permanente. 


“ Popularidade de controladores PID: 
% Desempenho robusto em uma ampla gama de condições de 
funcionamento; 
% Simplicidade funcional, que permite uma operação de forma simples 


e direta (basta escolher os parámetros K,, T; e Ta). 
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E. CONFIGURACAO E VARIACAO DO PID 


NEP — ] . | 


% A estrutura ideal (na equação matemática clássica presente em 
documentos técnicos) do controlador PID tem pouca utilização em 


aplicações práticas; 


* As estruturas de implementação do PID podem ser classificadas: 
o pelo tipo de interação entre as ações P, I, D; 


o pela distribuição das ações P e D na malha de controle; 


< PID ideal 
UCs) 
E(s) 


=G.(s)=K, Es 
T.s 


l 
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Е. _ CONFIGURAÇÃO E VARIAÇÃO DO PID 


[О > PID IDEAL 


proporcional integral 


e(t) u(t) 


derivativo 


- l de(t) 
«oe Ken | edt «T, = | 


Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


= _ CONFIGURACAO E VARIACAO DO PID 
PID SERIE 


derivativo 
u(t) 


proporcional 


integral 


O cde aes Hs 
E(s) Ts 
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= _ CONFIGURAÇÃO E VARIAÇÃO DO PID 
mix PID PARALELO 


proporcional U(s) 


=G (s) =K + m + К |5 
Е( 5) 5 


inteeral 


u (t) 


derivativo 
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= _ CONFIGURACAO E VARIACAO DO PID 


ШЕЕ PI+D 


r(t) e(t) 
id integral proporcional 


u(t) f yq) 
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derivativo 
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i CONFIGURAÇÃO E VARIAÇÃO DO PID 


ЖҮЛ» I+PD 


proporcional 


u(t) р 


inteeral 


r(t) © e(t) @ 


derivativo 


Eng. de Controle e Automação – Rejane de Barros Araújo/André Ferreira 


E CONFIGURAÇÃO E VARIAÇÃO DO PID 


EE» OBSERVAÇÕES 


% Os parámetros de projeto dos controladores são diferentes e as funções 


de transferência não são sempre equivalentes; 


% No PID série, os zeros do controlador somente podem ser reais, 
enquanto que no PID paralelo e ideal, podem ser reais ou complexos 


conjugados; 


% A característica de polos complexos conjugados pode ser útil quando o 


sistema a ser controlado tem polos complexos conjugados dominantes; 
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E CONFIGURAÇÃO E VARIAÇÃO DO PID 


> OBSERVAÇÕES 


4 O PID série é denominado de interativo, pois o ajuste de parcela 


derivativa afeta a parcela integral; 


% A estrutura paralela é a mais real, porém seus parâmetros não têm 


interpretação física; 
“ Quando se trabalha com conceito de polos e zeros, a estrutura série 


tem a vantagem de que cada zero do PID está relacionado com um 


parámetro do controlador (T; ou T ¿); 
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E CONFIGURAÇÃO E VARIAÇÃO DO PID 


NED OBSERVAÇÕES 


% Um controlador PID ideal pode sempre ser representado como um 


controlador PID série, cujos parâmetros são calculados como: 


ETST A. TT 
KeK a Т, =T+T, T, =" 
T; 1 +, 


% Caso se implemente o controlador somente com as parcelas P, Pl ou 


PD, entáo as formas ideal e série sáo equivalentes; 


% Na estrutura I+PD, o “kick” da parcela proporcional é subtraindo; 
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E SINTONIA DE CONTROLADORES PID 


% A sintonia do controlador PID envolve a seleção de valores adequados 


dos ganhos deste controlador; 


% A sintonia "ideal" do PID é um procedimento subjetivo e depende da 


planta a ser controlada; 


% Diferentes técnicas vêm sendo proposta pela literatura de controle ao 
longo dos anos e os resultados indicam que: 
o 30% dos controladores instalados operam em modo manual; 
o 30% das malhas são inconsistentes (projetos mal elaborados); 
o 25% das malhas utilizam o ganho default; 


o 30% das malhas apresentam problemas nos equipamentos. 
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E SINTONIA DE CONTROLADORES PID 


“ Controle servo: рага uma mudança degrau unitário na referência. 
Depende da dinâmica da planta e da sintonia do controlador. 
Rastreamento de referência (reference tracking): 


o Overshoot, tempo de subida, tempo de acomodação. 


4$ Controle regulatório: na presença de uma perturbação de carga na 
variável de saída, tem-se uma típica resposta desejada de malha 


fechada (disturbance rejection); 
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E SINTONIA DE CONTROLADORES PID 


“ Métodos empíricos: utilizam resultados de simulação ou prático do 
processo; 
o Ziegler-Nichols: curva de reação e curva de oscilação; 
o Cohen-Coon; 


o Chien, Hrones and Reswick; Hàgglund-Astróm; outros 


* Métodos analíticos: utilizam a função de transferência do modelo do 
processo e/ou índices de desempenho de malha fechada; 
o Alocação de polos; 
o Lugar geométrico das raízes; 
o IMC; GPC; GMV; outros 
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- SINTONIA DE CONTROLADORES PID 


EE AJUSTE DE PARÂMETROS 


« Maior desafio na utilização do controle PID; 


% Problema ainda maior quando não se tem muito conhecimento do 


comportamento dinâmico do sistema a ser controlado; 
“ Sintonia é geralmente feita utilizando métodos empíricos: 


* Ex.: Método de Ziegler e Nichols (ganho crítico e curva de reação); 


% Servem como uma primeira tentativa de sintonia do PID. 
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E SINTONIA DE CONTROLADORES PID 


EDD MÉTODOS DE ZIEGLER-NICHOLS 


* a) Método do Ganho Crítico: 


% Utilizado em malha fechada; 
4$ Aplica-se um controle puramente proporcional e aumenta-se о 
ganho do sistema até que o mesmo fique no limiar da estabilidade 


(oscilação sustentada); 


% Nesta situação utiliza-se o valor do ganho crítico (К„) e do período 


crítico (P,) para ajustar os parámetros do PID (K., T; e Ta); 
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E SINTONIA DE CONTROLADORES PID 


ED MÉTODOS DE ZIEGLER-NICHOLS 


* a) Método do Ganho Crítico: 


e 
CUT Processo E 


Esquema de controle para utilização do método do ganho crítico 


y(t) Per 


+ 


тоѕ Araújo/André Ferreira 


E SINTONIA DE CONTROLADORES PID 


NED METODOS DE ZIEGLER-NICHOLS 


4$» a) Método do Ganho Crítico: 


e [ws [= | 9 


0.45K,, Pa 
1.2 
0.6K P p 
o fo | mp | mg 
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E SINTONIA DE CONTROLADORES PID 


ED MÉTODOS DE ZIEGLER-NICHOLS 


% b) Método da Curva de Reação: 


« Utilizado em malha aberta; 
* Aplica-se um sinal do tipo degrau observando-se se a resposta do 
sistema possui forma de “S” (característica típica de muitos sistemas 


industriais); 


4 Caso essa resposta seja observada, utilizam-se algumas 


características da curva em “S” para ajustar os parâmetros do PID. 
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E SINTONIA DE CONTROLADORES PID 


ED MÉTODOS DE ZIEGLER-NICHOLS 


% b) Método da Curva de Reação: 


y(t) 


YE (----------------------=еенннененанне) 


Reta tangente no ponto 


de máxima inflexáo 


yo | 7 


A to t2 
ZN 
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E SINTONIA DE CONTROLADORES PID 


ED MÉTODOS DE ZIEGLER-NICHOLS 


% b) Método da Curva de Reação: 


rn 
ЕИ NUN 


PI 0.97 
A po 
PID 1.27 20 0 
m [Wap = | % 


* Exemplo 2: qe 
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E SINTONIA DE CONTROLADORES PID 


“ Os métodos de Ziegler-Nichols podem ser usados como uma primeira 


tentativa de sintonia do controlador PID. 


* Outros métodos de ajuste podem ser encontrados na literatura: 
*• Cohen-Coon 
* Tyreus-Luyben 
* Chien-Hrones-Reswich 


+» Astrom 
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E SINTONIA DE CONTROLADORES PID 


% Regra de Cohen-Coon 
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E SINTONIA DE CONTROLADORES PID 


“ Para fazer um “ajuste fino” nos parámetros do controlador PID deve-se 
levar em conta que, geralmente, o seguinte comportamento do sistema 


em malha fechada para uma entrada degrau é observado: 
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- EXERCÍCIOS 


4$ Projete um controlador proporcional para o processo de 2º ordem, que 
atenda as seguintes especificações: ovsershoot de 15% e tempo de 
acomodação de 0.4 segundos. 


1 


Gls) = (os + DDGS + 5) 


*• Projeto um controlador Pl para o processo de 1º ordem, que atenda as 


seguintes especificações: overshoot de 10% e tempo de pico de 1 


ў 
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segundo. 
| 


ls) = 642) 


- EXERCÍCIOS 


% Utilizando-se o simulink, simule as seguintes estruturas do controlador 
PID. Observe os sinais de saída e controle. Utilize dois osciloscópios. 

* PID ideal 

* PID série 

* PID pararelo 

* PI+D 


& 1+PD 
2,25 


G,(s) = ——  _ 
р(5) 52 +25 + 1 
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- EXERCÍCIOS 


% Utilizando-se o simulink, sintonize o controlador PID pelos seguintes 
métodos. 
% Ziegler-Nichols (Método do ganho crítico) 
% Ziegler-Nichols (Método da curva de reação) 


4 Cohen-Coon 


apis) = ç +s 
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- EXERCÍCIOS 


€ Desenvolva em código no Matlab, o projeto do controlador PID ideal, 
sintonizando pelos seguintes métodos: 
4$ Ziegler-Nichols (Método do ganho crítico) 
% Ziegler-Nichols (Método da curva de reação) 


+ Cohen-Coon 


apis) = ç + Ds 
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F. SATURACAO E ANTIWINDUP 


% Os processos industriais estão submetidos a restrições; 


% Um controlador opera num intervalo limitado de 1-5 V ou 4-20 mA, uma 
válvula não pode abrir mais que 100% e menos que 0%, um motor 


funcionando como atuador tem um limite de velocidade, etc; 


4 A entrada real da planta é temporariamente diferente da saída do 
controlador. Quando isto acontece, se o controlador é inicialmente projetado 
para operar em uma região linear, o comportamento dinâmico do sistema de 


malha fechada deteriora-se em relacáo ao desempenho linear esperado; 


4 Esta inadequação de desempenho (saturação na malha direta) é 


denominada windup. 
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yr(t) e(t) u (t) шї t 
| E 


Wmax(t) , Up(t) > Umax(t) 
ut) =+ up(t) , Uminlt) S up (t) S umax(t) 
Wmin(t) , t Ut) < Uminkt) 
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F. SATURACAO E ANTIWINDUP 


% Exemplo 3: considere a planta contínua de 2º ordem representada pelo 
seguinte modelo: 


1 
(s + 1)(5.175 + 1) 


Gp (s) = 

% sendo controlada por um controlador Pl+D sem e com a presença do 
limitador no sinal de controle (umin = 0 V е Umax = 5 V). Considerando 

a planta com uma dinámica de 1° ordem, K, = 1, т = 5.17 560 = 15, е 


aplicando a 1? regra de Ziegler-Nichols, obtém-se a seguinte sintonia: 


К. = 6.204, Т; = 2 e T4 = 0.5. d 
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Е... SATURAÇÃO E ANTIWINDUP 


* Exemplo 3: continuação 


вака гей га, controle 
E № 
iA LS ue 


sem limitador com limitador 


Simular no simulink 
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- TÉCNICA ANTIWINDUP 


% Quando a variável de controle atinge o limite de operação máximo (ou 


mínimo) do atuador, ocorre a saturação do sinal de controle; 


% A saturação faz com que a malha de realimentação seja de certa forma 
quebrada, pois o atuador permanece no seu limite máximo (ou mínimo) 


independentemente da saída do processo; 


* A consequência disto é que a resposta transitória do sistema torna-se 
lenta e oscilatória, características extremamente indesejáveis em um 


processo industrial; 
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- TÉCNICA ANTIWINDUP 


Se um controlador com ação integral é utilizado, o erro continua a ser 
integrado e o termo integral torna-se grande, ou seja, tende a carregar-se 


demasiadamente (windup); 


Para que o controlador volte a funcionar na região linear é necessário que 


o termo integral se "descarregue"; 


Deve-se esperar que o sinal de erro mude de sinal e, por um período de 


tempo aplicar na entrada do controlador um sinal de erro oposto; 


A técnica back-calculation impedi que o integrador continue a se 


carregar quando a saturação no sinal de controle ocorre. 
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- TÉCNICA ANTIWINDUP 


e()=y,() — уб) | ut) 


Cawlt) 
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- TÉCNICA ANTIWINDUP 


% A técnica back-calculation funciona da seguinte forma: quando a saída 
do atuador satura, o termo integral é novamente calculado de forma que 


seu valor permaneca no limite do atuador; 


% É vantajoso fazer esta correção não instantaneamente, mas 


dinamicamente com uma constante de tempo Taw; 


* Deve-se ter cuidado na seleção de Taw, especialmente em sistemas com 


ação derivativa; 


% Na prática deve-se ter Tw > Ta e menor que T;. Regras sugeridas na 


literatura de controle de processos são: Taw = 4/ TiTa ou Taw = Тү; 


É. TÉCNICA ANTIWINDUP 


“ Exemplo 4: mesmo exemplo anterior com a técnica antiwindup. 


sión (Wanqa contre 
no м [I 
"n № М „ 


sem antiwindub com antiwindup 


Simular no simulink 
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F. LUGAR GEOMÉTRICO DAS RAÍZES (LGR) 


* O LGR (Root Locus) é uma representação 
gráfica dos polos de malha fechada, quando 


um parâmetro do um sistema varia. 
e Proposto por Evans (1948; 1950); 


e Control System Synthesis by Root Locus 
Method in AIEE Transactions, v. 69, 1950. 
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LUGAR GEOMÉTRICO DAS RAÍZES (LGR) 


* Artigo sobre o LGR 
publicado na AlEE 


Transactions em 1950. 


Control System Synthesis by Root 
Locus Method 


WALTER R. EVANS 


MEMBER л 


Synopsis: The root locus method deter- 
mines al] of the roots of the differential 
equation оё a control system by a graphical 
plot which readüy permits synthesis for 
desired transient response ог frequency 
response. The hase points for this plot 
оа the complex plane are the zeros and poles 
of the open loop transfer function, which are 
readily available. The locus of roots is a 
plot of the values of e which make this 
transfer function equal to — 1 as loop gain is 
increased from zero to infinity. The plot 
can De established in approximate foem by 
inspection amd the significant parts of the 
locas caleulated accurately and quickly by 
use of a simple device, For maltipie loop 
systems, one solves the innermost Joop 
first, which them permits the next loop to be 
svad by asother root locus plot. The 
resultant plot gives a complete picture of 
the system, which is particularly valuable 
for mnusmal systems or those which have 
wide variations im parameters, 


HE root locus method is the result of 

an effort to determine the roots of the 
differential equation of a control system 
bv using the concepts now associated 
with frequency response methods, The 
roots are desired, of course, because they 
describe tbe natural response of the sys- 
tem. The simplifying feature of the 
control system problem is that the open 
loop transfer function is known as a 
product of terms. Each term, such as 
1/(14T5), сап be easily treated in the 
same manner as ап admittance such as 
1/(RHjx). It is treated as a vector in the 
sense used by electrical engineers in 
solving a-c circuits. The phase shift and 
attenuation of a signal of the form e" 
being transmitted is represented by 
1/(14-T4) in which s in general is a com- 
px number. The key idea in the root 
locus method is that the values of s 
which make transfer function around the 
loop equal to —1 are roots of the differ- 
ential equation of tbe system. 
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The opening section im this paper, 
Background Theory, outlines the overall 
pattern od analysis. The following sec- 
tion on Root Locus Plot points out the 
great usefulness of knowing factors of the 
open loop transfer function in finding the 
TOOLS. 

The graphical nature of the method 
requires that specific examples be 
used to demonstrate the method itself 
under the topics: Single Loop Example, 
Multiple Loop System, and Corrective 
Networks The topic Correlation with 
Other Methods suggesta methods by 
which experience in frequency methods 
сап be extended to this method. The 
topic Other Applications includes the 
classic problem of solving an mth degree 
polynomial Finally, tbe section on 
Graphical Calculations describes the key 
features of a plastic devie called а 
“Spirale”, which permits calculations to 
be made from direct measurement on the 
plot. 


Background Theory 


The over-all pattern of analysis can be 
outlined before explaining the technique 
of sketching a root locus plot. Thus 
consider the general single loop system 
shown in Figure 1 

Note that each transfer function is of 
the form KG(a) in which K is а static gain 
constant and G(s) is a function of the 
complex number. In general, G(s) has 
both numerator and denominator known 
infactored form. The values of a which 
make the function zero or infinite can 
therefore be seen by inspection and are 
called zeros and poles respectively. The 
closed loop transfer function can be ex- 
pressed directly from Figure 1 as given in 
equation 1 


Esrans— Control System Synthesis 


a OK) 
петка) Касаба) а 


The problem of finding the roots of the 
differentia] equation bere appears in the 
form of finding values of s which make 
the denominator zero. After these values 
аге determined by tbe root locus method, 
the denominator can be expressed in 
factored form, The zeros of the function 
&,/8, сал be seen from equation 1 to be the 
zeros of Gale) and the poles of Gale). 
The function can now be expressed as 
shown in — 2 


The constant K, and the exponent y 
depend upon the specific system hut for 
control systems y is often zero апа K, is 
often 1, 

The full power of the Laplace Trans- 
form? or an equivalent method now cam 
be used. The transient response of the 
output for a unit step input, for example, 
is given by equation 3 


=. 
Rih == Y Aja" (3) 
it 
The amplitude A, is given by equation 
4 


a-[; Т! 14) 
5 =" 


The closed loop frequency response, on 
the other hand, can be obtained bv sub- 
stituting s-jw into equation 2. For- 
tunately, the calculation in finding A, or 
бј) involves the same problem of 
multiplying vectors that arises in making 
a root locus plot, and can be calculated 
quickly from the resaltant root locus 
plot. 
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AIEE TRANSACTIONS 


F. LUGAR GEOMÉTRICO DAS RAÍZES (LGR) 


“É um método “poderoso” para análise e projeto do 
comportamento do sistema com relação a sua estabilidade e 


características de resposta transitória. 


* Traz informações relevantes sobre o comportamento do 
sistema (que podem ser complicadas de analisar para 


sistemas de ordem superior olhando simplesmente para sua 


função de transferência). 
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* О LGR pode ser usado para descrever qualitativamente o 
desempenho de um sistema quando determinados parámetros 
variam: 

* Ex: efeito da variação de um ganho sobre máximo percentual 


de sobressinal, tempo de acomodação e tempo de pico. 


% Mostra também claramente as faixas de estabilidade e 


instabilidade do sistema, além dos valores de ganho que levam o 


sistema a um limiar de estabilidade. 
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Forward 
. transfer 
` “| o š 
Actuating incio 
signal 


Input Output 


R(s) KG(s) Cís) 
1+ KG(s)H(s) 


Feedback 
transfer 
function 


(a) (b) 


+ KG(s)H(s)> FT de malha aberta; 


* (a) Sistema de malha fechada; 


* (b) FT equivalente de malha fechada. 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


мттитонр 


* Polos de malha aberta podem ser facilmente encontrados e 
geralmente não dependem dos ganhos do sistema. 
e Polos de KG(s)H(s). 


“ Polos de malha fechada são mais difíceis de serem calculados: 
• É necessário obter o denominador da FT de МЕ; 


e Sofrem influência do valor do ganho do sistema; 


KG 
° Polos de S ree oHG 
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` O PROBLEMA DO SISTEMA DE CONTROLE 


Forward 
А transfer 
Actuating fanclion 
Input signal 
Ej s) 


Output 
C(s) 


R(s) KG(s) Cís) 
1 + KG(s)H(s) 


Feedback 
transfer 
function 


(a) (b) 


_ Ме (5) 
Юс (s) 


G(s) 


ro = ——K NEH) 
Рс(5)рһ (s) + KNç(s)Nu(s) 

Мң (5) 

Оң (5) 
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H(s) = 


мттитонр 


río = — КМоб)0н (S) 


D¿(s)Dy(s) + KNç(s)Nr:(s) 


* Zeros de T(s): 


e Zeros de G(s) e Polos de H(s). Podem ser obtidos facilmente. 


** Polos de T(s): 
e Responsáveis pela resposta transitória e estabilidade do sistema; 


e Dependem do valor de K. 


% LGR: 
e Mostra graficamente a posição dos polos de T(s) quando o “O 
KG(s)varia. 
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f; REPRESENTAÇÃO VETORIAL DE NÚMEROS 
COMPLEXOS 


% Números complexos: 
e Forma cartesiana: o + Jo 


e Forma polar: MZ0 


“ Qualquer número complexo pode ser descrito graficamente como um 


vetor. 
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COMPLEXOS 


+ Se um número complexo é substituído em uma função complexa F(s), 


o resultado também é um número complexo: 


F(s) =s+a ——»  F(s)= (c +a)+ jo 


$ = o + JO 


4 2 possíveis representações: 
jo 
s-plane 
jo 
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COMPLEXOS 


4 Conclusão: 


e (s + a)é um número complexo que pode ser representado por um 


vetor desenhado do zero da função até o ponto s. 


“ Exemplo 1: (s + 7); 5+; 


• É um número complexo desenhado do zero da função (—7) até о 


ponto s = 5 + j2. 
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COMPLEXOS 


~ Generalizando para uma função F(s) qualquer: 


F(s) = [ils +2) _ [I[fatfres c2mplexBs d8 numerad?r 
С MC + pj) E | fatBires celimplexPls dEl denPiminad Er 


* Na forma polar: 


Р módulo dos 2егоѕ 1 IS Zi 
e Módulo: M =P mocuo dos zeros = оин A 

П módulo dos polos П? , |(5+ру)| 
e Fase: 0 = У, дпеишВ dBls хегЁБ — > дпоиВ dis pel O 


= Ziza) da “(s + р;) 
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E; REPRESENTAÇÃO VETORIAL DE NÚMEROS 
COMPLEXOS 


* Exemplo 2: 


(s+1) 


e Para F(s) = Ыы? 


calcule seu valor para s = —3 + j4 


« Contribuição dos polos е zeros: 
Zero em -1: V20.116.62 N 
Polo em 0: 5z7126.9* 
Polo em -2: /17 7104.0? 


j s-plane 


“ Resultante na forma polar: 


20 
М 20 = л 16.6°—126.9°—104.0°=0.217—114.3° 


5117 
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f; REPRESENTAÇÃO VETORIAL DE NÚMEROS 


COMPLEXOS 
+ Exercício 
e Рага F(s) = т calcular о seu valor no ponto. 
s = —7 + J9 
e Resposta: 


—0.0339 — j0.0899 = 0.096 — 110.7º 
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- DEFININDO O LGR 


~ Exemplo 3: 
e Câmera de segurança com rastreamento automático (usada para 


seguir objetos em movimento): 


Subjects Motor Camera 
position Sensors Amplifier and camera position 
RS) + C(s) 


e к» 
K e. AN 


R(s) K C(s) 
52+ 105+ K 


where K = K, K> 
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m DEFININDO O LGR 


K 
52+ 105+ K 


where К = K, K> 


* Posição dos polos com a variação do ganho K: 


K Pole 1 Pole 2 


0 -10 0 

5 —9.47 —0.53 
10 —8.87 -113 
15 -8.16 -1.84 
20 -7.24 -2.76 
25 28 -5 
30 —5 + j2.24 5:32 24 
35 —5 + j3.16 -S- j3.16 
40 —5 + j3.87 —5 – 3.87 
45 —5 + j4.47 -5 — j4.47 


50 -545 -S-js 
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DEFININDO O LGR 


R(s) K C(s) 
52+ 105+ K 


where К = K, K> 


* Posição dos polos com a variação do ganho K: 
jo 


s-plane 


K=0 5 1015 20 20 15 105 


-10 -9 8 -7 6 -5 -4 3 -2 -I 
-ji 
30X 32 
35 X -j3 
40X a 
45X A 
K=50X 
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E DEFININDO O LGR 


% LGR do sistema: 
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[Р DEFININDO O ГОВ 


s% Análise: 


e Fornece a posição dos polos рага K > 0. 


e Polos são reais e distintos para K < 25: 


o Sistema superamortecido. 


e Polos reais e iguais para K = 25: 


o Sistema criticamente amortecido. 


e Polos complexos conjugados рага K > 25: 
o Sistema subamortecido. 
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Е. PROPRIEDADES ро LGR 


% Partindo da FT de MF de um sistema de controle: 


T(s) = KG(s) 


1 + KG(s)H(s) 


o Os polos de Т (5) podem ser obtidos a partir de: 


Кб (5)Н (5) = —1 = 1Z(2k + 1)180º к= 0,+1,+2, +3,... 


° Uma forma alternativa de representar esses polos é dada por: 


IKG(s)H(s)| = 1 ZKG(s)H(s) = (2k + 1)180° | 
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E. PROPRIEDADES DO LGR 


* Para determinado polo que atenda a restrição de ângulo (múltiplo 


impar de 180º), o seu ganho associado pode ser obtido usando: 


1 


“COMO 
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F. PROPRIEDADES DO ГОВ 


* Exemplo 4: 


R(s) 


— K(s +3)(5 + 4) n 
KG(s)H(s) — (Gt06Gt2) salis 
— К (5 + 3)(s + 4) 


(1 + K)s^ + (3 + 7K)s + (2 + 12K) 
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Ё. PROPRIEDADES ро LGR 


* Exemplo 4: continuação 
e O ponto —2 + j3 é um polo de T(s)? 


o Fazer o teste do ângulo! 


0, + 0; — Өз — 0, = 56.31? + 71.57? — 90º — 108.43? -Cross?) 


jo 


o Conclusão: 


—2 + j3 não é um polo de MF de T(s). à рше 12 


Não é um ponto sobre o LGR 
Não é um polo de MF para 


algum valor de K 
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Não é múltiplo fmpar de 180º 


E PROPRIEDADES DO LGR 


R()) + O K(s + 3) (s +4) C(s) 
x . ГЕ QS (з + 1) (5+2) 
* Exemplo 4: continuação ш Шай 


° О ponto —2 + j (2/5) é um polo de T(s)? 
O, + 6, — 04 — 0, = 19.4712? + 35.2644º — 90? — 144.7356? = — 180º 


e (Conclusão: 


—2 + j (7) é um polo de MF de Т (s). — ^u 
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Ё. PROPRIEDADES ро LGR 


* Exemplo 4: continuação 


e Qual o valor do ganho que fornece o polo de MF em —2 + j (7/4) ? 
K 1 _ 1  [[módulB dEs pele 
—|G(s)IH(sS) M  [[módull dfs zergs 
s-plane 
l -2 + j 0,7071 
NP) 


Lab, | 5 (1.22) ú 
“LL (212)(122) ` 
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E PROPRIEDADES DO LGR 


“ Exercício 1: Dado um sistema com realimentação unitária que possui 
a seguinte função de transferência: 


K(s + 2) 

G(s) = ——— 
(52 + 4s + 13) 

e a) Calcule o ângulo de G(s) no ponto (—3 + j0) determinando а 

soma algébrica dos ângulos dos vetores desenhados a partir dos 


zeros e dos polos de G(s) para o ponto dado; 


• b) Determine se o ponto especificado em (a) está sobre o LGR; 


e C) Se o ponto especificado em (a) estiver sobre o LGR, 0 


o ganho K, utilizando os comprimentos dos vetores. 
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F- REGRAS PARA ESBOÇAR O LGR 


* 1) Número de ramos: 
e Cada polo de malha fechada se move quando o ganho do sistema 


varia. 


e Definindo ramo como o caminho que um polo percorre, então 
existirá um ramo para cada polo de malha fechada. 
o O número de ramos do LGR é igual ao número de polos de 


malha fechada do sistema. 
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- REGRAS PARA ESBOCAR O LGR 


K 
«^ 1 " 
* 1) Numero de ramos: 


Kz0 5 10 15 20 
-10 9 -8 -7 6 -5 43 -2 -l 


F- REGRAS PARA ESBOÇAR O LGR 


* 2) Simetria: 
e Se os polos complexos de malha fechada não são complexos 
conjugados, isso resultaria em um polinômio com coeficientes 


complexos. 


e Sistemas físicos não possuem coeficientes complexos em suas 
funções de transferência. 


o O LGR é simétrico com relação ao eixo real. 
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F- REGRAS PARA ESBOÇAR O LGR 


Ris) K C(s) 
5% | ia` 
* 2) Simetria: 


K=0 5 10 15 20 
-10 -9 -$ -7 4 
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Ж REGRAS PARA ESBOÇAR O LGR 


* 3) Segmentos no Eixo Real: 
e À contribuição angular dos polos e zeros complexos de malha 


aberta sempre irá se anular. 


e À contribuição angular dos polos e zeros reais de malha aberta 


vale 0º ou 180º. 


je 
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Ж REGRAS PARA ESBOÇAR O LGR 


* 3) Segmentos no Eixo Real: 
e Para K > 0, qualquer ponto do eixo real que ficar à esquerda de 
um número impar de singularidades (polos ou zeros) localizadas 


também no eixo real é um ponto do LGR. 


Jo 
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F- REGRAS PARA ESBOÇAR O LGR 


* 3) Segmentos no Eixo Real: 
e Exemplo 5: 


R(s) 
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F- REGRAS PARA ESBOÇAR O LGR 


* 4) Pontos Iniciais e Finais: 


e Onde o LGR começa (ganho nulo) e termina (ganho infinito)? 


e Analisa-se a FT de malha fechada para as duas situações: 


KNç(s)Du(s) 


TG) CODO + KNç (SN, (5) 
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Ж REGRAS PARA ESBOÇAR O LGR 


* 4) Pontos Iniciais e Finais: 
e Para K > 0: 


K Nc (5) Оң (5) 
De (s)Du(s) + К№(5) № (5) 


| 


_ КЇїє(5)Юн (5) 
TG) = p (урн (5) + € 


Т (5) = 


e O LGR começa nos polos de G(s)H(s) (FT de malha et 
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Ж REGRAS PARA ESBOCAR O LGR 


* 4) Pontos Iniciais e Finais: 


e Para K oo: 


K Nc (s)Dgu(s) 
Dg; (s)Du(s) + К№(5) № (5) 


| 


T(s) = K Nc (s)Du(s) 
€ + К№(5)№ (5) 


Т (5) = 


e O LGR termina nos zeros de G(s)H(s) (FT de malha bet 
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F- REGRAS PARA ESBOÇAR O LGR 


* 4) Pontos Iniciais e Finais: 
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F- REGRAS PARA ESBOÇAR O LGR 


* 5) Comportamento no infinito: 
e Aregra anterior diz que o LGR começa nos polos de malha aberta 


e termina nos zeros de malha aberta. 


e О que fazer quando o número de polos de malha aberta não é 
igual ao número de zeros de malha aberta”? 
o Uma FT G(s) pode ter polos ou zeros infinitos; 
o Polo infinito: 


G(s) — со quando s > оо. 


o Zero infinito: 
G(s) > 0 quando s > oo. © 
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F- REGRAS PARA ESBOÇAR O LGR 


* 5) Comportamento no infinito: 
e Exemplo 6: 
оа) G(s) = ѕ 


Zero em 0 e Polo no infinito. 


ob) G(s) = t/s 


Polo em 0 e Zero no infinito. 


* Todo sistema tem um número igual de polos e zeros, caso sejam 


considerados todos os polos e zeros finitos e infinitos. Q 
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F- REGRAS PARA ESBOÇAR O LGR 


<+ 5) Comportamento no infinito: 


e Suponha que uma função de transferência possui um número de 


polos finitos maior do que o de zeros finitos. 


e Como o LGR começa nos polos de malha aberta e termina nos 
zeros de malha aberta, é necessário conhecer a localização de 


todos os polos e zeros do sistema. 
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F- REGRAS PARA ESBOÇAR O LGR 


<+ 5) Comportamento no infinito: 
e Nesta caso, onde estão localizados os zeros infinitos”? 


e O LGR se aproxima de assíntotas quando o lugar geométrico se 


aproxima de infinito. 


e А equação dessas assíntotas é definida pelo ponto onde 


interceptam o eixo real (o4) e pelo ángulo (0, ). 
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F- REGRAS PARA ESBOÇAR O LGR 


<+ 5) Comportamento no infinito: 


> polos finitos — >, zeros finitos 


"  N? de polos finitos — N? de zeros finitos 


o (2k + 1)n 


^  N? de polos finitos — Nº de zeros finitos 


К = 0,-+1,+7,... 
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F- REGRAS PARA ESBOÇAR O LGR 


* 5) Comportamento no infinito: 
° Exemplo 7: 


o Esbocar o LGR para o sistema abaixo: 


Kis + 3) 
a s(s + 108 + 2)(s + 4) Ш 


о Polos finitos: 4 


о Zeros finitos: 1 


о Zeros infinitos: З 
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F- REGRAS PARA ESBOÇAR O LGR 


* 5) Comportamento no infinito: 


K(s + 3) 
sS ж 
T sis + 1s + 2)(5 + 4) im 


o Definição das assíntotas: 


e Exemplo 7: 


Ris) + 


ds (2k + 1)x 
' Жыен MED REST EEE A 
#finite poles — #finite zeros 
a = 122-49 _4 =x/3  fork=0 
a 4-1 3 БЫ E 


= g fork=1 


= Sm/3 for k = 2 
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F- REGRAS PARA ESBOÇAR O LGR 


% 5) Comportamento no infinito: 


Kis +3) 
sis + 1Ks + 2)(s + 4) 


e Exemplo 7: 


R(s) + 


jo 


s-plane 


_C1-2- 4 - (3) _ 
.Ci-2-9-C3. 


Ga 


4 
3 


(2k + 1) a Asymptote 
Oa aT Fw oh 
#finite poles — #finite zeros 

= л/3 for k = 0 

= g fork = 1 

= 5л /3 for k = 2 


Asymptote 
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F. LGR 


* Exercício 2: Esboce о LGR e suas assíntotas para o sistema com 


realimentação unitária que possui a seguinte função de transferência: 


K 


G(s)= Үс у(х + y + 6) 
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F. REFINANDO O Е$ВОСО DO LGR 


* As 5 regras mostradas anteriormente permitem se obter um esboço do 


LGR rapidamente; 


* Caso seja necessário um detalhamento maior, deve-se obter pontos 


específicos do LGR e seus ganhos associados. 
* Para maior precisão no LGR deve-se obter as seguintes informações: 
. Pontos no eixo real onde o LGR entra ou sai do plano complexo; 
e Pontos de cruzamento no eixo jæ; 


. Ângulos de partida e chegada dos polos e zeros complexos. Ө 
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* Considere o esboco do LGH feito com as regras mostradas 


anteriormente: 


K(s — 3)(s — 5) 


AS TIED 


e Qual o valor dos pontos o; e oz? 
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"PONTOS DE ENTRADA E SAÍDA DO EIXO REAL 


* Método usando diferenciação: 


e Deve-se minimizar ou maximizar o ganho K. 


e Para qualquer ponto no LGR vale a relação: 
o 1 
... G(s)H(s) 


• Para pontos sobre o eixo real, onde os pontos de entrada e saída 


podem existir, tem-se s — o: 


1 
m © G(c)H(o) Ө 
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“PONTOS DE ENTRADA E SAÍDA DO EIXO REAL | 


* Método usando diferenciação: 


Е 1 
mE G(o)H(o) 


-3 -2 -9g,-1 0 1 2 305 4 5 
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"PONTOS DE ENTRADA E SAÍDA DO EIXO REAL 
e Partindo dos polos e zeros de malha aberta: 
_ K(s—3)(s—5) К(5? – 85 + 15) 
оса (5 + 1)(5 +2) (52+ 35 + 2) 
e Para todos os pontos до LGR, KG(s)H(s) = —1 e para os 


segmentos по eixo real, s = o. Então: 


К(о? — 80 + 15 (52 
2 o m ve (c^ + Зс + 2) 
(c^ +30 + 2) (c? — 8с + 15) 


e Diferenciando K com relação a се igualando a zero: 


ак (11e? — 266 — 61 B Ее: 
E nn dao as LPs © 


do (02 – 80 + 15)? 
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"PONTOS DE ENTRADA E SAÍDA DO EIXO REAL 


* Método da transição: 


e Não precisa utilizar a diferenciação; 


e Os pontos de entrada e saída satisfazem a relação: 


n 
1 > 1 
ctz Z +р 


o onde z; e p; são, respectivamente, o negativo dos zeros e polos 


de G(s)H(s). 


As 
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SENE SEE! — K(s - 3)(5 — 5) 
esq ain AA (s + 1)(s + 2) 


—1.45 
110? – 260 — 61= 0 —> ped e 
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PONTOS DE CRUZAMENTO NO EIXO JO 


5 É o ponto do LGR que separa a operação estável do sistema da 


operação instável. 


* O valor do ganho do LGR no ponto de cruzamento pode ser calculado 
utilizando o critério de Routh-Hurwitz: 
e Forçar uma linha de zeros na tabela de Routh para obter o ganho 


no ponto de cruzamento; 


e Resolver o polinômio na linha acima da linha cheia de zeros para 


achar o ponto de cruzamento. O 
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PONTOS DE CRUZAMENTO NO EIXO JO 


Kis + 3) 
s(s + IMs + 2)(s + 4) 


* Exemplo 8: 


jo 


s-plane 
Asymptote 


jl 


Asymptote 


Asymptote 
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PONTOS DE CRUZAMENTO NO EIXO JO 


EN . R(s) + < — N33 _ _ A) 
» Exemplo 8: жЕ стщ 


о Obter a ЕТ de MF do sistema: 


K(s + 3) 


TO = DTL 
(s) s^ + 7s? + 14s? + (8 + K)s + ЗК 


o Montar a tabela de Routh: 


gi 1 14 3K 


p 7 8+K 
s? 90 — K 21K 
E —K? — 65K + 720 

90 — K 
ҳо 21К 
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PONTOS DE CRUZAMENTO NO EIXO JO 


* Exemplo 8: 


gs 1 14 3K 
$5 7 8+K 
s 90 – K 21K 
5! К? — 65K + 720 
90 — K 
s? 21K 


e Para valores positivos de K, a linha s! pode ser formada por zeros: 


—K^—65K +720 = 0 — P X>, K=9.65 © 
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PONTOS DE CRUZAMENTO NO EIXO JO 


* Exemplo 8: 


1 14 3K 


sá 
s 7 8+K 
52 90 — К 21К 
ғ -K? — 65K + 720 

90— K 


= 


21K 


с^ 


e Formando o polinômio par da linha s? рага K = 9.65: 
(90 — K)s? + 21K = 80.355 + 202.7 20 —»  s=+j1.59 


e Conclusão: 
o Sistema estável para O < K < 9.65 © 
o Ponto de cruzamento do eixo jæ: s = +j1.59 
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t: ÂNGULOS DE PARTIDA E DE CHEGADA DOS 
POLOS E ZEROS COMPLEXOS 


* ОТОН começa nos polos de MA e termina nos zeros de MA; 


+ Para um LGR mais preciso é necessário calcular o ángulo de partida 


dos polos complexos e o ângulo de chegada dos zeros complexos. 


* Assume-se um ponto e no LGR próximo ao polo complexo que se 


deseja calcular o ángulo de partida: 
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t: ÂNGULOS DE PARTIDA E DE CHEGADA DOS 
POLOS E ZEROS COMPLEXOS 


* А soma dos ángulos de todos os polos e zeros finitos com relação ao 
ponto e deve ser um múltiplo ímpar de 180º. 


e 


s-plane 


—04 +05 +03 —0, —05+06= (2k + 1)180 


с 
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t: ÂNGULOS DE PARTIDA E DE CHEGADA DOS 
POLOS E ZEROS COMPLEXOS 


+ Da mesma forma, pode-se calcular o ángulo de chegada de um zero 


complexo: 


jo 


—04 +05 +03 —0, —05+06= (2k + 1)180 
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t: ÂNGULOS DE PARTIDA E DE CHEGADA DOS 
POLOS E ZEROS COMPLEXOS 


+ Exemplo 9: K(s * 2) 


S 
1 (536732512) = 


• Polos complexos: —1 + j1 


s-plane 


—f — 85 + 84 — 84 = 180º 


1 
—8, — 90° + tan”! 9 — tan” (5 — 180* 


| 


0, = —251.6° = 108.4° rlocus(G) 
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а 


e EXEMPLOS DO LGR 
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- EXEMPLOS DO LGR 
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F. EXEMPLO COMPLETO 


% Considere o sistema mostrado abaixo: 


K(s + 2) 


= (s? — 45 + 13) 


• a) Desenhe o esboço do LGR 


e b) Ache o cruzamento com o eixo jo 


) 
e C) Calcule o ganho K no cruzamento com o eixo jo 
e d) Ache o ponto de entrada no eixo real 

) 


e е) Calcule o ângulo de partida dos polos complexos O 
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É EXEMPLO COMPLETO 
A) DESENHE O ESBOÇO DO LGR 


K(s + 2) 


E (s2 — 4s + 13) 


e Polos de malha aberta: 2 + j3 


e Zero de malha aberta: —2 


jar 
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: EXEMPLO COMPLETO 


мттитонр 


B) АСНЕ O CRUZAMENTO COM О Е!ХО JO 


K(s + 2) 


s (s2 — 4s + 13) 


e Obtenha a FT em MF e depois utilize o critério de Routh-Hurwitz: 


G)- 1x6) Fs + ОК 13) 
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E. EXEMPLO COMPLETO 
B) ACHE O CRUZAMENTO COM O EIXO JO 


e Obtém-se uma linha cheia de zeros para K = 4. 
e Da linha s? com K = 4, tem-se: 


e Conclusão: 
o Ponto de cruzamento com o eixo jæ é 


$2+21=0 —> 5= +01 
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- EXEMPLO COMPLETO 
C) CALCULE O GANHO K NO CRUZAMENTO COM O EIXO JO 


e Do item (b), sabe-se que K = 4 no cruzamento com o eixo jæ. 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


ie EXEMPLO COMPLETO 
D) ACHE O PONTO DE ENTRADA NO EIXO REAL 


í — 1 " 1 
0-2 0q-2-;j3 0-2+j3 


0? — 40 + 13 = (0 + 2)(20 — 4) 


g^ — Ат + 13 = 262: —B 


2 = Е - F4 
а +40 21=0>0=( 5 
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u EXEMPLO COMPLETO 
D) ACHE O PONTO DE ENTRADA NO EIXO REAL 
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в EXEMPLO COMPLETO 


мттитонр 


- D) ACHE O PONTO DE ENTRADA NO EIXO REAL 


_ 1 |. K(s+2) 
la С (s)H (s) EDU (s? — 4s + 13) 
(s? — 4s + 13) 
Е Ш 
% Fazendo s = —7: 
"-— (C7)2 —4(—7) +13) + _ (49 + 28 + 13) 


((77) + 2) E (—5) 


K — 18 
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É EXEMPLO COMPLETO 
E) CALCULE O ÂNGULO DE PARTIDA DOS POLOS COMPLEXOS 


jo 


* Desenhe os vetores para um ponto є perto do polo 


complexo: 


2. ângulo dos zeros — > ângulo dos polos = —180° 


3 
tan`! (5) — 90° — 0 = —180° 2 


0 = 126.9° 
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F. EXEMPLO COMPLETO 
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E. LGR USANDO O MATLAB 


K(s + 2) 
>> num=[1 2] G(s) — (s? — 4s + 13) 
num = 
1 2 
>> den=[1 -4 13] 
den = 
1 -4 13 


>> g-tf(num,den) 


$^2-45+13 


Continuous-time transfer function. 
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- LGR USANDO O MATLAB 


>> rlocus(g) : Root Locus 

I T 

4 - > E -| 

/ мА 
£ | x 

/ | 

p^ 2- j | | 
кеј / | 
c / | 
[e] / l 
© | | 
Ф | I 
o, | 

Об Oee сс RITE NS A 
& | | 
D | | 
[э] | 
E | 

g 2. | E 
E | 

-4 E | _ 

-6 | | | | | | | | | 
-16 -14 -12 -10 -8 -6 -4 -2 0 2 4 


Real Axis (seconds!) 
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E LGR USANDO O MATLAB 


Root Locus 
6 | | 
^^ P I b N 
| / EN 
System: g z є 
Саіп: 18 ү, | 
= 2 Pole: -7 / | 
W E " Ў l = 
о Damping: 1 / | 
9 Overshoot (96): 0 I i 
ЕЯ Frequency (rad/s): 7 | | 
2 0 E 2 = 
< | | 
> | | 
Е \ | 
© -2- | = 
- | 
4. | = 
-6 [ | | | | | | | | 
6. 4 «12 -1O -8 -6 -4 -2 0 2 4 


Real Axis (seconds!) 
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- LGR USANDO O MATLAB 


Root Locus 
6 | | 
Ыы System: g ' м A 
Gain: 4.01 | N 
Pole: -0.00647 + 4.58i | X 
= / Damping: 0.00141 | 
wr / Overshoot(%):996 - 
5 .. Frequency (rad/s): 4.58 | 
3 x 
- | 
о 0 | Sci E A E ln A 
E | | 
> | | 
E \ 
© -2 + Ñ | | 
E ; 
-4 x | 
-6 [ | | | | | | | | 
46- i da dO — 38 -6 -4 2 0 2 4 


Real Axis (seconds `!) 
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F. ASSUNTOS COMPLEMENTARES 


4 LGR Generalizado: 


e Utilizado para verificar como os polos em MF do sistema variam 


em funcáo de um parámetro do sistema diferente do ganho. 


e Ex.: 
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F. ASSUNTOS COMPLEMENTARES 


* LGR para sistemas com realimentação positiva: 


TO KG(s) 
ISE KOHO 


e Como a função de transferência de malha fechada muda, as regras 


para traçar o LGR também se alteram. O 
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UNIDADE IV: PROJETO DE COMPENSADORES 
USANDO O LGR - PARTE 2 


o Engenharia de Controle e Automacáo 
1 o Prof. André Ferreira/ Prof? Rejane de Barros Araújo 


PROJETO DE COMPENSADORES USANDO O LGR 


+O LGR permite escolher o ganho adequado para 
atender — determinadas ` especificações ае 
desempenho transitório, desde que essas 


especificações façam parte do LGR. 
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“PROJETO DE COMPENSADORES USANDO O LGR 


Ho 


> Ponto А pode ser atendido 
com ajuste de ganho, 


diferente do ponto B. 


X= Closed-loop pode 


X= А | ENT 
к + Apenas com a variação do 
(a) ganho é possível atender 
etn Polsat В mA apenas a especificação de 


máximo  sobressinal, mas 
não o que se deseja com 
relação ao tempo O 
acomodação. 

(h) 
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PROJETO DE COMPENSADORES USANDO O LGR 


+ Solução: 


e Compensar o sistema com polos e zeros adicionais; 


e O sistema compensado deve possuir ит LGR que passe pelos 
polos dominantes que atendam às especificações de 


desempenho desejadas; 


e Deve-se garantir que os novos polos e zeros não afetem 


significativamente o comportamento dos polos dominantes do 


sistema. o 
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PROJETO DE COMPENSADORES USANDO O LGR 


+ Compensadores: 
e Usados para melhorar, de forma independente, as respostas 


transitória e em regime permanente do sistema. 


* Ajuste de “ganho puro”: 
e TK le, TM(%) 
e LK М) T e 


“+ Compensadores Dinâmicos: 
e Não apresentam resposta instantânea (não são ganhos e sim FTs); 
e Conseguem atender a ambas as especificações "0 


regime permanente) simultaneamente. 
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E COMPENSADORES DINÂMICOS 
CONFIGURAÇÕES 


х a) Cascata: 


Cascade Original 
compensator controller Plant 


* b) Em realimentação: 


Original 
controller Plant 


Ri ‚ 


Feedback 
compensator 
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- COMPENSADORES DINÁMICOS 


+ Ambas as configurações mudam o posicionamento dos polos e 
zeros de MA, criando um novo LGR que deve passar pela 


localizacáo dos polos de MF desejados. 


+ Tipos: 
e Lag 
• Pl 
e Lead 
e PD 
e Lead-Lag 
e PID 
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Ё. COMPENSADORES IDEAIS 


+ Integrador puro > melhora erro de regime permanente; 


+ Derivador puro Ə melhora o comportamento transitório, porém é 


“ruidoso”. 


> Observação: Esses compensadores são mais caros, pois só podem 
ser implementados com redes ativas (usando amplificadores 


operacionais e fontes de potência externas). O 
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- COMPENSADORES МАО IDEAIS 


+ Podem ser implementados com componentes passivos 
(resistores e capacitores), porém não são integradores e 


derivadores puros; 


“+ São mais baratos e não necessitam de fontes de potência 


externas; 


“+ Não conseguem zerar o его de regime nos casos em que о 


compensador ideal consegue. @ 
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. COMPENSADORES IDEAIS X COMPENSADORES МАО 
IDEAIS 


“ À escolha entre um ou outra solução envolve custo, 
peso, desempenho desejado, função de transferência 
e a interface entre o compensador e o restante do 


sistema. 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ а) Compensador Integral Ideal (РІ): 


e Sistema Original Gain Plant 
ele 
jo 


—04 ES 0- E Өз — (2k + 1)180º 


s-plane 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ а) Compensador Integral Ideal (РІ): 


e Integrador puro (afeta o transitório do sistema original) 


Compensator Plant 


R(s) + C(s) 


x G(s) 


jo 


—01 — 0, — 04 — 0, = (2k + 1)180? 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ а) Compensador Integral Ideal (РІ): 
e Integrador puro (sem afetar o transitório do sistema original) 
Compensator Plant 


R а 


jo 


$-plane 


8c 


SER 


—0, — 0; — 03 — Ope + Ө = (2k + 1)180? 
A 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


* а) Compensador Integral Ideal (РІ): 
e Sistema sem compensador atendia especificações de desempenho 


transitório, mas com erro de regime permanente; 


e Inclusão de um compensador */< atendia ás especificações de 


regime permanente, mas náo a de transitório; 


e O compensador K(s+)/. atende ambas as especificações, quando 
a está perto do polo na origem, recebendo o nome de 
COMPENSADOR INTEGRAL IDEAL. o 
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USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 1: Para o sistema abaixo, deseja-se um amortecimento de 


6,470,174. 


Gain Plant 
Els) \ 
mm (s+1)(s+2)(s+10) a 


*** Proj) LGR exl.m *** 


Eng. de Controle e Automação – André Ferreira e Rejane de Barros Araújo 


f; MELHORANDO O ERRO DE REGIME PERMANENTE 
| USANDO СОМРЕМ$АСАО EM CASCATA 


+ Exemplo 1: continuação. 


Root Locus 
T 


Imaginary Axis 


| | | | | | 
-35 -30 -25 -20 -15 -10 -5 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 1: continuação. 


+ Dar zoom e achar polos para $420,174. 


Root Locus 
T T T T 
Sr | d 
System: g1 
Gain: 161 
2- Pole: -0.686 + 3.88i 
o Damping: 0.174 
& Overshoot (%): 57.4 
D 0 š " Frequency (rad/sec): 3.94 
E 
D 
А] 
£ 
= 
-4L 
6 ы 
| | | | 
-2 -1.5 -1 -0.5 0 


Real Axis 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 1: continuação. 
+ Para К=161: 
* polos dominantes: 
o —0,686 + j 3,88 
o 6= 0,174 
o 0, = 3,94 rad/s 
о М(%) = 57,4% 


* 3º polo: 
o —11,6 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
| USANDO СОМРЕМ$АСАО EM CASCATA 


+ Exemplo 1: continuação. 


161 


+ ЕТ de МА: K6G)HG)=——5G+2G+10 


161 161 


GON GE 


8,05 


1 
PIE 


p 


е(20) = = 0,1105 
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USANDO COMPENSAÇÃO EM CASCATA 


Exemplo 2: Deseja-se novamente um amortecimento de $420,174, 


porém agora o compensador não é mais um ganho puro. 


Compensator Plant 


*** Proj LGR ex2.m *** (mostrar sisotool) © 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
| USANDO СОМРЕМ$АСАО EM CASCATA 


Exemplo 2: continuação. 


Root Locus 
T 


Imaginary Axis 


l l | | | | 
-35 -30 -25 -20 -15 -10 -5 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


* Exemplo 2: continuação. 


“Dar zoom e achar polos para $470,174. 


Root Locus System: g2 
T 


Gain: 155 LL 
Pole: -0.67 + 3.79i LA 
6r Damping: 0.174 J 
Overshoot (%): 57.4 
Frequency (rad/sec): 3.85 
4L = 
2L -| 
2 
< 
> 0 
o 
E 
[9)] 
3 
= DE 
-4- 
-6L 
l | | [Í | 


Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


* Exemplo 2: continuação. 


Step Response 
T 


I 
Sistema não-compensado 
Sistema com integrador ideal 


Amplitude 


| | | 


| 
0 5 10 15 20 25 30 35 
Time (sec) 
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E MÉTODO ALTERNATIVO 


+ Controlador РІ: 
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USANDO COMPENSAÇÃO EM CASCATA 


“+ b) Compensação em atraso (Compensador Lag): 


e A compensação usando um integrador ideal requer o uso de uma 


rede ativa. 


e Para o projeto de um compensador usando uma rede passiva, 
pode-se colocar o polo e o zero perto da origem, conforme a 


figura abaixo: 


jo 
IZel > [рє 


(S + Ze) 


8 5 ажр) 


“le Pe 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ b) Compensação em atraso (Compensador Lag): 
e Diminui o ganho do sistema em altas frequências, sem reduzir o 
ganho em baixas frequências (melhora a precisão em regime 


estacionário); 
e Reduz a largura de faixa; 
e O sistema fica mais lento; 


e Por causa da diminuicáo do ganho em altas-frequéncias o ganho 
total do sistema pode ser aumentado, melhorando a precisáo O 
regime permanente 
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USANDO COMPENSAÇÃO EM CASCATA 


+ b) Compensação em atraso (Compensador Lag): 


e Sistema original: 


Gain Plant 


(s zs + 22) + + + 
K ¡q íÍ _ __ —— —— ELCESXEOEnnÓ£rrrren"qánibOCeseà] AEG 
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USANDO COMPENSAÇÃO EM CASCATA 


+ b) Compensação em atraso (Compensador Lag): 


e Sistema com compensador em atraso: 


Compensator Plant 
K(s + z.) 
(s + p.) 


KZZ ...\ 7 
Kag -( = JE 
P1P2 --. J Dc 
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USANDO COMPENSAÇÃO EM CASCATA 


+ b) Compensação em atraso (Compensador Lag): 
e Se o polo e o zero do compensador lag estão perto, a 
contribuição angular do compensador para o ponto P é 


aproximadamente zero: 


o Não causa um efeito significativo na resposta transitória do sistema original. 


jo 


s-plane 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ b) Compensação em atraso (Compensador Lag): 
e Com essa configuração o ganho para o polo dominante no ponto 
P é aproximadamente o mesmo para o sistema compensado e 


para o sistema não compensado: 


Zc 


Кос = Kyo — > Kyo 
Dc 


e Pode-se melhorar o ganho estático do sistema não compensado 


através da razão entre o polo e o zero do compensador lag. O 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


* b) Compensação em atraso (Compensador Lag): 
e Deve-se lembrar que, para não afetar muito o desempenho 


transitório, o polo e o zero do compensador devem estar próximos. 


° Exemplo: 


—0,01 ME 
р, = 0,000) Ке = Ко 0,001 ^ “Pe = 10 Ko 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ b) Compensação em atraso (Compensador Lag): 
e Embora um integrador ideal elimine o erro em regime permanente 
para uma entrada degrau, um compensador lag (com um polo que 
não está na origem) irá melhorar a constante de erro estático de 


um fator “/,.. 


e Escolhendo z e р. perto da origem também garante uma 


alteração mínima na resposta transitória do sistema compensado. 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 3: Melhorar o erro de regime permanente de um fator 10, 
usando um compensador lag para o sistema mostrado no exemplo 


anterior. 


Compensator Plam 


K(s +2.) l 
(s + 1)(5 + 2)(s + 10) 


(s t p.) 


*** proj LGR ex3.m *** e 
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MELHORANDO O ERRO DE REGIME PERMANENTE 


USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 3: continuação 


+ Do sistema original sabe-se que: 


| 161 
E= 161 > K, = lim 


5—0 (s + D(s + 2)(5+ 10) — 


e, (20) = = 0,1105 


ItK 


po 


+ Deseja-se para o sistema compensado: 


eo (2) 
10 


€, (00) = = 0,01105 
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8,05 


.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 3: continuação 


+ Assim, tem-se: 


1 
е9) = Ko. e. (a) 


1+ Kpc = 90,497 > Kpc = 89,497 


Ze Ky 89497 .. a 
e и. HE б 


po 


+ Selecionando arbitrariamente p.: 


р. = 0,01 =z, = 11,2 х 0,01 


z, & 0,111 
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USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 3: continuação 
+ Obter o LGR do sistema compensado e achar os polos para 
E450,174. 


Compensator Plant 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
| USANDO СОМРЕМ$АСАО EM CASCATA 


+ Exemplo 3: continuação 


Root Locus 
T [| 


Imaginary Axis 


| | | | | | | 
-35 -30 -25 -20 -15 -10 -5 0 Б 10 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 3: continuação 


Root Li 
System: 03 | 
Саіп: 156 NC d 
Pole: -0.674 + 3.82i 
6r Damping: 0.174 - 
Overshoot (%): 57.4 
Frequency (rad/sec): 3.87 
4- m - 
2- А 
E 
< 
E 
S 0 x 
D 
[4] 
Е 
э -| 
-AL -| 
-6L -| 
| | | | | 
-2 -1.5 -1 -0.5 0 0.5 
Real Axis 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


Y 


+ Exemplo 3: continuacáo 


* Para К=156: 


Y 


œ polos dominantes: 
o —0,674 +j 3,82 
o 6= 0,174 
o 0, = 3,87 rad/s 
o M,(%) = 57,4% 


+ Outros polos: 
o —11,6 
o —0,101 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 3: continuação 


(s + 1)(5 + 2)(s + 10) (s + 1)(s + 2)(s + 10)(s + 0,01) 
A / 


805 89,497 
eo) _ 0,1105 | 000 | 


Polos 
—0,686 + j 3,88 —0,674 tj 3,82 
dominantes 
р -0,101 


0,111 


рс 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 3: continuação 


Step Response 
T 


Sistema não-compensado 
Sistema com compensador LAG 


1L 
\ Дм 
V 
o 08 
© 
2 
E 
5 0.6 - 
0.4 - 
0.2 
0 | | | | 
0 5 10 15 20 25 


Time (sec) 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 3: continuação 


Step Response 
1 4 T T T T T 


Sistema com integrador ideal 
Sistema com compensador LAG 


0.8 H 4 


Amplitude 


0.6 4 


0.4 | Е 


0.2 Е 


| | | 


| 
0 5 10 15 20 25 30 35 
Тіте (ѕес) 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 
+ Exemplo 3: continuação 
o Compensador lag não é capaz de eliminar o erro em regime 
permanente para uma entrada degrau da mesma forma que o 
controlador Pl (embora possa ter um desempenho bastante 


aproximado ao desejado). 


o Compensador lag é formado por uma rede passiva (não usa 


amplificadores operacionais e fontes externas). 
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USANDO COMPENSAÇÃO EM CASCATA 


“ a) Compensador Derivativo Ideal (PD): 
o Caso o LGR do sistema original não passe pela localização dos 
polos que atendam às especificações de desempenho transitório, 
a inclusão de um derivador ideal na malha direta pode fazer com 
que o LGR do sistema compensado apresente o comportamento 


desejado. 


Compensator Plant 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO СОМРЕМ$АСАО EM CASCATA 


+ Exemplo 4: Sistema não compensado (£,=0,4): K 
(s + 1)(s + 2)(s + 5) 
Root Locus 
10 I 
8L 
6L 
4L 
o 2p 
2 
© 0 
4L 
-6L 
-8L 
-10 | | | | | | l 
-16 -14 -12 -10 -8 -6 -4 -2 4 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: Sistema não compensado (£,=0,4): K 
(s + 1)(s + 2)(s + 5) 


Root Locus 
| T | T System: g1 | 
4L Gain: 23.7 — 
Pole: -0.938 + 2.15i m ud 
Damping: 0.4 
3r Overshoot (96): 25.4 
Frequency (rad/sec): 2.35 


Imaginary Axis 
o N 
T T 
N 


| | | | | | | | 
-1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0 
Real Axis 
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MELHORANDO O ERRO DE REGIME PERMANENTE 


USANDO COMPENSAÇÃO EM CASCATA 
+ Exemplo 4: Sistema não compensado (6,=0,4): K 
(s + 1)(s + 2)(s + 5) 
System: mf1 
Peak amplitude: 0.866 
RM E Sep Responso 
Р 234 
08 | DEA NÉ 
0.7 E 
| System: mf1 1 
0.6 - | Settling Time (sec): 3.75 = Е. ET 
| el) OE 0,2967 
Ф 0.5 F} A 
Ë 0.4 L E 


0.3- 


0.2 | 


Time (sec) 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: Compensação 1 (64=0,4): K(s + 2) 


(s + 1) (5 + 2) (5 + 5) 


Root Locus 
10 
8 - 
6 + 
4L 
a 2r 
£ 
= 0 >x 
£ 
[9)] 
[t] 
E +25 
-4 
-6L 
-8- 
= 0 l l l [ 
-6 -5 -4 -3 -2 -1 0 1 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: Compensação 1 (64=0,4): K(s + 2) 
(s + 1)(s + 2)(5 + 5) 
Root Locus 
ER E 
System: 92 
6- Gain: 51.3 
Pole: -3 + 6.88i 
5. Damping: 0.4 
Overshoot (96): 25.4 
Frequency (rad/sec): 7.51 
zm 
Ë 
2L 
1- 
0 o 
q 


Real Axis 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO СОМРЕМ$АСАО EM CASCATA 
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+ Exemplo 4: Compensação 1 (64=0,4): 


Amplitude 


System: mf2 

Peak amplitude: 1.14 
Overshoot (%): 25.4 
Attime (sec): 0.459 


Step Response 
T 


System: mf2 
Settling Time (sec): 1.12 


Time (sec) 
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K(s + 2) 
(s + 1)(s + 2)(s + 5) 


p im "UEFI ijai 
e eree ` ' 


1 
е(оо) = = 0,0888 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: Compensação 2 (64=0,4): K(s + 3) 
(s + 1)(s + 2)(s + 5) 
Root Lo 
10 
8L 
6L 
4 - 
wm 2r 
ы 
s or © £ —— —x 4 
=Й = 
-6L 
-8 L 
-10 | | | | i 
-6 -5 -4 -3 -2 -1 0 1 
Real Axis 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: Compensação 2 (64=0,4): K(s + 3) 
(s + 1)(s + 2)(s + 5) 


Root Lo, System: g3 
Gain: 35 


10 СЕ 
Pole: -2.43 + 5.55i 
aL Damping: 0.4 
Overshoot (96): 25.3 
a Frequency (rad/sec): 6.06 
| E 
AL 
o 2Һ 
E 
Ë 0 Б x O- — AAA 
£ 
о 
б 
Е -2Ь 
-4 |- 
-6L 
-8 L 
-10 l l | | | | 
- -5 -4 -3 -2 -1 0 1 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: Compensação 2 (64=0,4): K(s + 3) 
(s + 1)(s + 2)(s + 5) 


System: mf3 


Step Response 
I T 


Peak amplitude: 1.16 
Overshoot (96): 27.3 
At time (sec): 0.566 K. = lim 35(s + 3) = 105 
=. 1 P s0(s+1)(s+ 2)M(s+ 5) Р 
4L System: mf3 | 
с a ЖИИ. — Settling Time (sec): 1.38 | — à 
dE po po ——, (20) a 0,087 
= el 20) = — Ë 
E 1+ K, 
E 
E 
« 4 
| | 
1.5 2 2.5 


Time (sec) 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 
+ Exemplo 4: Compensação З (64=0,4): K(s + 4) 
(s + 1)(s + 2)(s + 5) 


Root Lo 
10 
8L 
6L 
4L 
a 2p 
£ 
5  — —— H— A V EE 
£ 
D 
[t] 
E -2- 
-4L 
-6L 
-8L 
-10 l 
-6 5 4 -3 2 1 0 1 
Real A 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: Compensação З (64=0,4): K(s + 4) 
(s + 1)(s + 2)(s + 5) 
Root Locus 
10 | System: 94 EM 
Gain: 20.7 
8- Pole: -1.87 * 4.27i 
Damping: 0.4 
6- Overshoot (96): 25.4 
Frequency (rad/sec): 4.66 
4 m 
o 2- 
£ 
© 0 SO — *— AS — 
-4L 
-6L 
-8- 
-10 І 
-6 5 4 -3 2 1 0 1 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE - 
USANDO COMPENSAÇÃO EM CASCATA | 


+ Exemplo 4: Compensação З (64=0,4): K(s + 4) 
(s + 1)(s + 2)(s + 5) 


Step Response 
T T 


System: mf4 
Peak amplitude: 1.13 
Overshoot (%): 27.1 20,7(s + 4) 
; é | K, = lim —————noO—-——V = 8,28 
Attime (sec): 0.737 р == (s 4 ET 4+ 2)(s 4 5) 
System: mf4 -| 
Settling Time (sec): 1.8 
= = 0,1078 


Amplitude 


| | 
15 2 2;5 3 3.5 
Time (sec) 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: 


Sistema nào Е _ l 
Compensação 1 | Compensação 2 | Compensação 3 
dominantes | Polos —0,938 +j 2,15 —3 +j 6,88 -2,43+]5,55 -1,87 +j 4,27 
ominantes 


M,(%) 

desej ado 
mr umma masma mama mar [mm 
o ÉS na 
[ S*plo |  -&122 | sss 
Zero | x | x | 


Comentários | Aproximacáo de à == de | Aproximação de 
2º ordem OK 22 ordem OK 22 ordem OK 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 4: resposta ao degrau 


Step Response 
T T 


T 
Sistema náo compensado 
——— Compensação 1 
ur ——— Compensação 2 


————— Compensação 3 
1 L 
i N HESSE 


Ф 0.8 | la 
o 
= 
Ë 
< 06 
0.4 
0.2 
0 | | | | | | 
0 1 2 3 4 5 6 7 
Time (sec) 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


мтттонр 
má 


USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 5: Projetar um compensador derivativo ideal para que o 
sistema abaixo apresente um sobressinal máximo de 16% com 


redução de 1/3 no tempo de acomodação original. 


ЖЕ 


М„(%) = 16% = 0,16= e * 


—1,8326 = — 


1-2 


03403(1— 4)= č 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
| USANDO СОМРЕМ$АСАО EM CASCATA 


+ Exemplo 5: continuação. 


Imaginary Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 5: continuação. 


Root Locus 
T T 


System: 91 

Gain: 43.2 LA 
Pole: -1.2 + 2.05i NT | 

4 Damping: 0.504 + 

Overshoot (%): 16 

Frequency (rad/sec): 2.38 

2. = 4 


Ітадіпагу Ахіѕ 
о 
| 


-1.8 -1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0 
Real Axis 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 5: continuação. 
+ Рага K = 43,2: 
+ polos dominantes: 
o —1,2 +j 2,05 
o É = 0,504 
o 0, = 2,38 rad/s 
o М(%) = 16% 


+ 3º polo: 
o —7,59 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 5: continuação. 


Step Response 


1.4 
System: mf1 
Settling Time (sec): 3.47 
CS SS SP SS SS ESSES EE TT 
Ф $ Д 
кеј 
= 
E 
| | | i | | 
2 2.5 3 3.5 


Time (sec) 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 5: continuação. 


t, 3,335 "ER e ul 
sdes = 3 E EN = 1,1117s š m. O ES 

4 4 Со, = 3,5982 E 
Odes = Eo = LAT 3,5982 E 0,504 > OD = 7,14 rad/s 


Q4 = On [ — 2? = 6,1668 rad/s 


| 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 5: continuação. 
polos desejados = — 3,5982 + 6,1668 


+ Contribuição angular dos polos de MA: 


ZS|s=-3,5982+j6,1668 = 120,2627º 
> 8, =2752553° 


“(5 + 6)ls=-3,5982+j6,1668 = 68,7204° @ — » Ө, = —180° 


В..=95:2553° 
2(5 + 4)ls=-3,5982+j6,1668 = 86,2721º 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 5: continuação. 
jo 


ё = 0.504 
jl 
Desired 
-3.613 + /6.193 compensated 6 
š J 
dominant pole 

j5 
j4  s-plane 
ГА, 


Uncompensated f 
dominant pole —1-205 + /2.064 


ZE 4 <j 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 5: continuação. 


6,1668 


95,2553º 
| le 


3,5982 — z, 


6,1668 
3,5982 — Z4 = 108719 


tg0— tg(180? — 95,2553?) = 10,8719 


Z4 — 3,5982 — 0,5672 
6,1668 _ 6,1668 


HE. = I pa. 
gU T TUE sis Za tg 


z, = 3,03 
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MELHORANDO O ERRO DE REGIME PERMANENTE 


USANDO COMPENSAÇÃO EM CASCATA 
+ Exemplo 5: continuação. K(s + 3,03) 
| 4)(s + 6 
+ Sistema compensado TTE 


Overshoot (96): 16 
Frequency (rad/sec): 7.15 


Imaginary Axis 
о сл 
| 
| 
| 
| 
x 
| [E 
d 
l 


Real Axis 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO СОМРЕМ$АСАО EM CASCATA 


мтттонр 
má 


* Exemplo 5: continuação. 
% Sistema compensado 


Step Response 


I 
System: mf2 System: m 
; А Peak amplitude: 1.15 
Peak amplitude: 1.12 ; 
: Overshoot (96): 15 
Overshoot (96): 12.1 At time (sec): 1.67 
12r At time (sec): 0.529 id 
TT) System: mf1 
| __ Settling Time (sec): 3.47 
1E === = = =s 


| System: m2 
| Settling Time (sec): 1.15 


Amplitude 


Sistema náo-compensado 
Sistema com compensação derivativa ideal 


H H H l l 
0 0.5 1 1.5 2 2.5 3 3.5 
Time (sec) 
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4 


4.5 


K(s + 3,03) 
s(s + 4)(s + 6) 


E MÉTODO ALTERNATIVO 


+ Controlador PD: 


K 
К+ E zero do compensador 


K2 


2 


G¿(s) = K-s + Ki = К» S + 


К, > ganho do compensador em МА 


° Problemas: 
o Circuito ativo 
o Ruído no sinal devido à diferenciação 
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USANDO COMPENSAÇÃO EM CASCATA 


“+ р) Compensação em Avanço (Compensador Lead): 
o Vantagens: 

o Pode ser implementado com uma rede passiva e reduz o ruído 

de diferenciação. 
o Desvantagens: 

o À presença de um polo extra não reduz o número de ramos que 
tendem a cruzar o eixo jæ como no caso da inclusão de apenas 
um novo zero. 

[Zel < Incl 
(s + Zc) 
(s + pc) 
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Gols) = 


.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Para o sistema abaixo, reduza o tempo de acomodação 


pela metade, mantendo o sobressinal máximo em 30%. 


Compensator 


ss + 4)(5 + 6) 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
| USANDO СОМРЕМ$АСАО EM CASCATA 


+ Exemplo 6: continuação 


М,(%) = 30% > ¿= 0,358 


Root Locus 


Imaginary Axis 


-15 | | | 
-20 -15 -10 -5 


Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: continuação 


М,(%) = 30% > ¿= 0,358 


Root Locus 
15 i 
10 System: g1 
| Gain: 63.2 7 
Pole: -1.01 + 2.63i 
Damping: 0.357 
5L Overshoot (96): 30 | 
Frequency (rad/sec): 2.81 
E^ ш 
: Й 
= 0 x x 
£ 
[9)] 
[v] 
E 
-5L 4 
-10- aj 
-15 | | l I | | l l | | 


Real Axis 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: continuação 
* Para K=63,2: 
<+ polos dominantes: 
o –1,01 +j 2,63 
o € = 0,358 
o 0, = 2,81 rad/s 
o M.(%) = 30% 


+ 3º polo: 
o —7,99 
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MELHORANDO O ERRO DE REGIME PERMANENTE 


USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: continuação 


t. 3,976 
Usdes = 2 Es 2 sdes — 1,9885 


4 4 
todos 1,988 


Odes = 


= Oje, = 2,012 


to, —2,0122 q, = 


0, — polos desejados > — 2,012 + j 5,248 


Og = On 1-4 = 5,248 rad/s o 
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MELHORANDO O ERRO DE REGIME PERMANENTE 


мтттонр 
má 


USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: continuação 


io 


¿= 0,358 


polos desejados > — 2,012 + 5,248 


j4 s-plane 


polos sist. não compesado > — 1,01 + j 2,62 


s-plane 


{| 
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USANDO COMPENSACAO EM CASCATA 
* Exemplo 6: continuação 


« Projeto do compensador LEAD: 
o Arbitrar a localização do zero e calcular posteriormente o polo que 


faz o LGR passar pelo ponto desejado. 


o Contribuição angular dos polos em MA: 


ZSls=-2.012+j5,248 = 110,976" 


As + 4)]5- 2124/5248 = 69,253? 2 8, = 232,997° 


“(5 + 6)1;-_2,012+;5,248 = 52,768° Q 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 1 (continuação) 


s+5 
S + Dci 


Zei = 5 > Ga (S) = 
Ms + 5)ls=-2,012+j5,248 = 60,345º 


Ozc1 — Oper — k A @, = — 180? 


Opa = 60,345? — 232,997? + 180? 


Opc, = 7,348? 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 1 (continuação) 


Өг = 7,348? 


‚24 
= ——— 50,129 — 2,012) = 5,248 
tan( pci) Pei — 2,012 (Der ) 


Per = 42,71 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 
+ Exemplo 6: Compensador 1 (continuação) 


s+5 
Ge (S) = 


s + 42,71 


Root Locus 
T 


Imaginary Axis 
o 
& 
x 


| | | | | | | | 
140 -120 -100 -80 -60 -40 -20 0 20 40 60 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 1 (continuação) 
* Definir a função de transferência com o compensador LEAD 


s+5 
s + 42,71 


Gc1 (5) = 


>> gz-tfi[1 5],convi[1 42.71],сопу{[1 O],convi[i 3],[1 6])))) 


Transfer function: 
B + 5 


A+ ado 8" 3 + 451.1 a2 + X825 a 
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USANDO COMPENSACAO EM CASCATA 
+ Exemplo 6: Compensador 1 (continuação) 
+ Usar a função “rlocfind” para achar o ganho LGR no ponto 


desejado 


>> [k2,p2]-rlocfind(g2,-2.0124j*5.248] 


s+5 


Es É Em a Jem 
k2 — 165) = 4127 


1.4127е+003 


р2 = 


3.3919 
-2.0120 + 5.24801 
-2.0120 - 5.24801 
—5.1341 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 1 (continuação) 
+ Verificar se o amortecimento desejado é obtido. 


s+5 
s+42,71 


Gc1 (5) zx 


>> dampipz! 

Eigenvalue Damping Freq. (rad/s) 
—4.36e-001 1.00e+000 4.36e+001 
-2.О1е+000 + 5.25e+000i 3.58e-001 5.62e+000 
-2.О1е+000 – 5.25e+000i 3.58e-001 5.62=+ü00 
—5.13e+000 1.00e+000 5.13e+000 
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USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 1 (continuação) 
+ Para K=1412,7: 
e polos dominantes: 


o—2,012 + j 5,248 


e 3º e 4º polos: 
o —5,1341; -43,5519 


o Zero: 
o—5 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 1 (continuação) 


UR S 1.3 Sistema 
Overshoot (%): 30.5 Pa Ш s sponse 
At time (sec): 0.608 Tm At time (sec): 1.37 | | : - | 
ess Sistema não-compensado 
———— Compensador Lead 1 
di System: mf2 System: mf1 | 


Settling Time (sec): 1.96 Settling Time (sec): 3.97 


Amplitude 


Time (sec) 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
| USANDO COMPENSAÇÃO EM CASCATA 
+ Exemplo 6: Compensador 2 


s+4 
S + Dc 


LR D- 
“(5 + 4)ls=-2,012+j5,248 = 69,253? 


0,2 = б,» EI ba Ө, LM 180º 


Өр. = 69,253° — 232,997? + 180° 


Ө... = 16,256° 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 2 


б» = 16,256º 


5 
tan( 0.5) = —— — p 0,2916(p,, — 2,012) = 5,248 
( €5c2) p.; — 2,012 (Pez G. Gi s+4 


s + 20,01 


р. = 20,01 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 2 


s+4 
Eds — 
c2 
s + 20,01 
Root Locus 
60 
40 
20 
K^] 
£ 
a 0 Meroe a 
© 
S 
E 
-20 
-40 
-60 | | l l l l 
-70 -60 -50 -40 -30 -20 -10 
Real Axis 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 2 (continuação) 
* Definir a função de transferência com o compensador LEAD 


s+4 
s + 20,01 


Go (s) E 


>> g3-tfi[1 4],comei[1 ad0.01],comveir[i1 0],convr[i 4],[1 67111) 
E locus (дз! 


Transfer function: 
а + 4 


3^4 + 30.01 =” + 224.1 502 + 480.2 3 
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USANDO COMPENSACAO EM CASCATA 
+ Exemplo 6: Compensador 2 (continuação) 
+ Usar a função “rlocfind' para achar o ganho LGR no ponto 


desejado 


>> [k3,p3]7"rlocfind[g3,-Z.0124j*5.248) 


694.5274 


p3 = 


-21.9861 

-2.0120 + 5.2480i 
-2.0120 - 5.24801 
-4.0000 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 2 (continuação) 
+ Verificar se o amortecimento desejado é obtido. 


>> damp [p3] 

Eigenvalue Damping Freq. (rad/s) 
-2.20e-001 1.00e+000 2 .20e+001 
-2. 01e+000 + 5.25e4-0001 3.5Be-001 5.62e+000 
-2.01e4-000 - 5.2¿5e+0001 2.5Be-001 5. 62e+000 
–4 . 00e+000 1.00e+000 4.00e+000 
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USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 2 (continuação) 


+ Para К=694,527: 
e polos dominantes: 
o -2,012 +) 5,248 
e 3? e 4? polos: 
o —4; 21,9861 
° Zero: 


o-4 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 2 


System: mf3 System: mf1 
Peak amplitude: 1.29 Peak amplitude: 1.28 sponse 
Overshoot (%): 28.7 о Overshoot (%): 27.8 | | | 
id d uid ЛОШИ (За МЗГ Sistema não-compensado 
PR a —— Compensador Lead 2 
1.2 + | 
System: mf3 System: mf1 


Settling Time (sec): 1.98 Settling Time (sec): 3.97 


Amplitude 


0.6 | | | | _ 


Time (sec) 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
| USANDO СОМРЕМ$АСАО EM CASCATA 


+ Exemplo 6: Compensador 3 


s+2 
S + р.з 


Ze3 — 2 = С.з (5) = 
As + 2)1;-_2,012+5,248 = 90,131? 


бз — без — У Op = — 180° 


бз = 90,131° — 232,997º + 180° 


Өз = 37,134° 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 3 


бз = 37,134? 


= 0,757(p,4 — 2,012) = 5,248 
сз 2,012 5+ 2 


s + 8,943 


52 
tan( 05,3) = Dos — 2,012 
бз (s) == 


p,, = 8,943 
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f; MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 3 


s+2 
Ооз (5) = ———— 
3 
€ s + 8,943 
Root L 

10 

8- 

6- 

4L 
a 2p 
& 
S o ° 
Ё 
° 
Ë -2L 

-4L 

-6L 

SB 

-10 | 

-20 15 -10 5 0 5 

Real A 
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USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 3 (continuação) 
* Definir a função de transferência com o compensador LEAD 


s+2 
s + 8,943 


Ga(s) = 


>> gá-tfi[i1 2] ,conví[1 8.943], convili 0] ,convi[1 4],[1 6])))) 


Transfer function: 


3^4 + 18.54 3^3 + 113.4 Ho o4 + 214.6 8 
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USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 3 (continuação) 
+ Usar a função “rlocfind” para achar o ganho LGR no ponto 
desejado 


>> [k4,p4]-rlocfind(g4,-2.0124j*5.248) 


ка = 


344.4025 


-13.2766 
-2.0121 + 5.24801 
-2.0121 - 5.24801 
-1.6423 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 3 (continuação) 
+ Verificar se o amortecimento desejado é obtido. 


s+2 
s + 8,943 


Goa (s) EE 


>> damp (pá) 

Eigenvalue Damping Freq. (rad/3) 
-1.33e+001 1.00e+000 1.33e+001 
-2.01е+000 + 5.25e+0001 3.58e-001 5.62e+000 
-2.01e+000 - 5.25e+0001 3.58e-001 5. 62e+000 
-1.64e+000 1.00e+000 1.64e+000 
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t MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


Y 


+ Exemplo 6: Compensador 3 (continuação) 
+ Para K=344,4025: 


Y 


e polos dominantes: 


o—2,012 + j 5,248 


e 3? e 4? polos: 
o —1,6423; -13,2766 


e Zero: 


o —2 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 3 


System: mf 1 
Peak amplitude: 1.28 sponse 
System: mf4 Overshoot (96): 27.8 | | | 
: : At ti 11.37 
or T ү did as d Sistema náo-compensado 
vershoot (96): 14.2 => 
At time (sec): 0.727 | ш | 
с | | System: mf1 
| | System mia Settling Time (sec): 3.97 


l Settling Time (sec): 1.64 a 


Amplitude 


| | il | 


2 3 4 5 6 
Time (sec) 
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.. MELHORANDO O ERRO DE REGIME PERMANENTE 
USANDO COMPENSAÇÃO EM CASCATA 


+ Exemplo 6: Compensador 3 
+ Comparação entre os compensadores projetados 


Step Response 
1 4 T Ї T T T 
Sistema náo-compensado 
————— Compensador Lead 1 
12r — ———- Compensador Lead 2 
— ———- Compensador Lead 3 


Amplitude 


Time (sec) 
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MELHORANDO O ERRO DE REGIME PERMANENTE 
E A RESPOSTA TRANSITÓRIA 


NiTTUTOFID 
má 


“+ Combinando-se as técnicas mostradas anteriormente, deseja- 
se melhorar o desempenho transitório e o erro de regime 
permanente do sistema compensado simultaneamente: 


e 1) Melhora-se o transitório e depois ajusta-se o regime 


permanente. 


o Desvantagem: geralmente a velocidade da resposta 


diminui quando se ajusta o regime permanente. 
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"m MELHORANDO O ERRO DE REGIME PERMANENTE 
E A RESPOSTA TRANSITÓRIA 


+ Combinando-se as técnicas mostradas anteriormente, deseja-se 
melhorar o desempenho transitório e o erro de regime permanente 
do sistema compensado simultaneamente: 

e 2) Melhora-se primeiro o regime permanente e depois a resposta 


transitória. 


o Desvantagem: o ajuste da resposta transitória modifica as 
condições de regime permanente podendo piorar (o que não 
atenderia às especificações desejadas) ou melhorar a resposta 


em regime (geralmente não apresenta maiores problemas, o 


pode causar um aumento no custo do processo). 
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Е A RESPOSTA TRANSITORIA 


+ Método (1) será utilizado. 


„> a) Controlador PID 


K K 
2 1 2 
K Kıs + K, + Kas? Кз (s +з +) 
сд) = к ++ кув = Misto OT TR + К) 


Eng. de Controle е Ашотасао — André Ferreira е Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Etapas de projeto: 
• 1) Avalie o desempenho do sistema não compensado para 


determinar quanto deve-se melhorar a resposta transitória. 


° 2) Projete o controlador PD para atender às especificações de 
resposta transitória (localização do zero e ganho de malha 


aberta). 


e 3) Simule o sistema para garantir que as especificações foram 


atendidas. 


* 4) Caso as especificações de resposta transitória não sejám) 


satisfeitas, refaça o projeto. 
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MELHORANDO O ERRO DE REGIME PERMANENTE E 
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“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Etapas de projeto: 
° 5) Projete o controlador Pl que leva o sistema a atender às 
especificações em regime permanente. 


• 6) Determine Ki, K; e K3. 


e 7) Simule o sistema para garantir que todos os requerimentos 


(resposta transitória e regime permanente) foram atendidos. 


° 8) Reprojete os controladores caso as especificações obti 


° 


estejam fora do que é solicitado. 
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MELHORANDO O ERRO DE REGIME PERMANENTE E 
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A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Projete um controlador PID para que o sistema abaixo 
passe a operar com um tempo de pico igual a 2/3 do tempo de pico 


do sistema não compensado, com um sobressinal máximo de 20% e 


com erro de regime nulo para uma entrada degrau. 


| K(s + 8) 
(s + Hs + 6X5 + 10) 
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E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: LGR do sistema não compensado 


Root Locus 
15 | 
10 + 
5L 
£ 
> 
q 0 © 
© 
о 
[o] 
E 
-5 
-10 L 
-15 l | | | | 
-12 -10 -8 -6 -4 -2 0 2 
Real Axis 
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| MELHORANDO O ERRO DE REGIME PERMANENTE E 
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A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: LGR do sistema náo compensado 


Root Locus 
15 
10 m 
System: g1 
Gain: 122 
Pole: -5.41 + 10.6i 
5L Damping: 0.456 
Overshoot (96): 20 
E Frequency (rad/sec): 11.9 
> ob " FEN " | w 
E 
[9)] 
o 
E 
-5L 
-10- 
-15 


| | | l [| 
-10 -8 -6 -4 -2 0 
Real Axis 
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E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Usando a função "rlocfind 
>> [ki,pi]2rlocfind(gi1,-5.414j*10.6] 


ЕІ = 


122.1669 


pl = 


-5.4160 +10.59959i 
-5.4160 -10.5999i 
-8.1680 
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MELHORANDO O ERRO DE REGIME PERMANENTE E 
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A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 
+ Exemplo 7: continuação 
<+ Obtendo o coeficiente de amortecimento a partir do sobressinal: 


>> HMpdes-ü0.z; 


ecsi=- log (Mpdes) ¿sart (pitz+log (Mpdes)*2)] 
М,(%) = 20% > ¿= 0,456 


+ Usando a função "damp" para verificar o amortecimento: 


>> damp (pl) 

Eigenvalue Damping Freq. (rad/s) 
-5.42e+000 + 1.06e+001i 4.55e-001 1.19e+001 
-5.42e+000 - 1.06e+001i 4.55e-001 1.19e+001 
-8.17e+000 1.00e+000 8.17e+000 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 
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t MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: continuação 
+ Para K=122,17: 
° polos dominantes: 
o—5,416 + j 10,6 
об = 0,456 
o 0, = 11,9 rad/s 
o M,(%) = 20% 


• 3º polo: 
o —8,168 
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| MELHORANDO O ERRO DE REGIME PERMANENTE E 


NET TUTOFED 
mú 


A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Resposta do sistema não compensado 


Step Response 
T 


System: mf1 

1.2L Peak amplitude: 1.02 
Overshoot (%): 20.8 
At time (sec): 0.293 


Amplitude 


| | | | | | 
0 0.2 0.4 0.6 0.8 1 1:2 
Time (sec) 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


œ Exemplo 7: continuação 


« Tempo de pico do sistema não compensado: 


e Calculado: 


л 


t, = — 
p Oj 


• À partir do sistema simulado (lembrar que o sistema não é um 
sistema de 2º ordem): 


t, = 0,2935 


М„(%) = 20,8% © 


Eng. de Controle е Ашотасао — André Ferreira е Rejane de Barros Araújo 


| MELHORANDO O ERRO DE REGIME PERMANENTE E 


NET TUTOFID 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


* Exemplo 7: Especificações de regime transitório: 


2 
t E х 0,293 = 0,1953 s 


pdes 


T 
> @asaəs = 16,08 rad/s 


pdes 


Dades = 


16,08 


Oades — Ondes 1 — Č S Ondes = === 
н ü ; des JA — (0,456)? 


Ün = LOndes = 0,456 x 18,07 > Ojos = 8,24 


polos desejados > — 8,24 + 16,08 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


= 18,07 rad/s 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
|. .A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: continuação 


polos desejados > — 8,24 + 16,08 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Contribuição angular dos polos 


“(5 + 3)ls=-8,24+516,08 = 108,05? 
zs + 6)ls=-8,24+716,08 = 9/,93º 
ds + 10)].— 224.516,08 — 83,75? 


$ Ө, = 289,73° 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Contribuição angular do zero 


“(5 + 8)ls=-s,24+716,08 = 90,85º 
» Ө, = 90,85? 


+ Contribuição angular do zero do PD 


£ 0. + ë, -X 8, = —180° 


90,85? + O. — 289,73? = —180º 


z. = 18,88° 


Eng. de Controle е Automação – André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
| A RESPOSTA TRANSITÓRIA USANDO РО+РІ (PID) 


* Exemplo 7: Cálculo do zero do PD 


16,08 
tan( 6.) = — 74 = 0,342(2, — 8,24) = 16,08 


G.pp(s) = (s+ 55,26) О 


E 5526 


Eng. de Controle e Automação – André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


* Exemplo 7: 


<+ Definir a FT com o compensador PD: 


>> gz-tficonvi[1 8],[i 55.256]],convi[1 3],convi[1 6],[1 10]1!1) 


Transfer function: 
B^2 + 63.26 3 + 442.1 


33 + 19 5^2 + 108 s + 180 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: LGR do sistema com compensador PD 


Root Locus 


Imaginary Axis 


-40 L 


-60 | | | | 
-200 -150 -100 -50 0 50 


Real Axis 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


| MELHORANDO O ERRO DE REGIME PERMANENTE E 


NET TUTOFED 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: LGR do sistema com compensador PD 


Root Locus 
60 Т | T T | I | 


T T 
System: 02 
Gain: 5.56 
ЛО Pole: -8.25 + 16.1! 
Damping: 0.456 
Overshoot (%): 20 
20 L Frequency (rad/sec): 18.1 


Imaginary Axis 


l | l | l l | l l 
-60 -50 -40 -30 -20 -10 0 10 20 30 40 
Real Axis 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: 
« Usando a função “rlocfind”: 


>> [k2,p2]-rlocfind(gz,-8.244j*16.08) 


p2 = 
-8.2402 +16.0799i 
-8.2402 —16.07991 
-B.0773 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


* Exemplo 7: 
+ Usando a função "damp" para verificar o amortecimento: 


>> dampipz] 

Eigenvalue Damping Freg. (rad/s) 
-8.24е+000 + 1.61e+001i 4.56e-001 1.81e+001 
-8.24е+000 - 1.61e+001i 4.56e-001 1.81e+001 
-8.08e-000 1.00e4+000 8. 08e+000 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: 
+ Para K=5,56: 
° polos dominantes: 


o-8,2402 + j 16,0799 


° 3º polo: 
o —8,0773 


* Zeros 
o —8; —55,26 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Simulação do sistema com compensador PD: 


Step Response 
T 


T T 
Sistema não-compensado 
Sistema com compensador PD 


Amplitude 


| | | 
0.6 0.8 1 1.2 


Time (sec) 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Simulação do sistema com compensador PD: 


Step Response 


System: mf2 T T | | | 
Peak amplitude: 1.13 pM Sistema náo-compensado 
Overshoot (76): 21.7 усеш Ѕіѕіета сот сотрепѕааог РО 
At time (sec): 0.171 Peak amplitude: 1.02 mp | 
— | Overshoot (96): 20.8 
At time (sec): 0.293 
Тир кык OO _ 
o 09r | | | 
Ko) 1 
2 
E | | 
< 06+ | | 3 
0.4 |- a 
0.2 | | | - 
0 i | | | | | | 
0 0.2 0.4 0.6 0.8 1 1:2 


Time (sec) 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


* Exemplo 7: 
“+ Colocando o zero do РІ perto da origem (-0,5): 


s + 0,5 
Gp (s) == " 


+ Definir a FT com os compensadores PD e РІ: 


> g3strf(vronv(rt1-0:5],.comve([ft^8]5[1:55.261)] ,conv([r-01];convt[l1-31;conv([p1-^6];T1:40])3))j 


Transfer function: 
3^3 63. 16 3^2 4.473.730 221 


874 P159 scr IUB 3^7 F' 180.8 o 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: LGR do sistema com compensadores PD e Pl. 


Root Locus 


Imaginary Axis 


-40 L 


-60 | | | | 
-200 -150 -100 -50 0 50 


Real Axis 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


* Exemplo 7: 


* Usando a função “rlocfind”: 


>> [k3,p3]=rlocfind(g3,-8.24+3*16.08) 


k3 = 
5.6270 
p3 = 
-8.0363 +16.15191 
-0 . U30 -10. di 
-8.0815 
-0.4729 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


* Exemplo 7: 


+ Usando a função "damp" para verificar o amortecimento: 


>> dampip3)] 


Eigenvalue Damping Freq. (rad/s) 


-8 .04e+000 + 1.62e+0011 4.45e-001 1.80e+001 
-8 .04e+000 – 1.62e+0011 4.45e-iI 1.80e+001 
-8.08er000 1. 00e+000 8. 08e+000 
= . 132-001 1.00e+000 4.73e-üU0U01 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


* Exemplo 7: 
+ Para K=5,62: 
° polos dominantes: 


о-8,0363 + j 16,1519 


° 3º e 4º polos: 
o —8,0815 


° zeros 
o —0,5; -8; -55,26 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 
má 


A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Simulação do sistema com compensadores PD e Pl: 


Step Response 
T 


1.4 | | 
Sistema náo-compensado 
Sistema com compensador РО 
ъё Sistema com controlador PID 
TN 
М 
o 08 
o 
2 
E 
< 0.6 
0.4 
0.2 
0 l l l l l 
0 1 2 3 4 5 6 


Time (sec) 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


MELHORANDO O ERRO DE REGIME PERMANENTE E 


мтттонр 


“A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Simulação do sistema com compensadores PD e Pl: 


Step Response 


I T T 


44, System: mf3 


Peak amplitude: 1.17 | = 
Overshoot (%): 16.5 Sistema não-compensado 


At time (sec): 0.174 Sistema com compensador PD 
Sistema com controlador PID 


stem: mf2 
ak amplitude: 1.135298 
'ershoot (96): 21.7.5 
time (sec): 0.171 |j | Peak amplitude: 1.02 

0.8L | | Overshoot (%): 20.8 
i | At time (sec): 0.293 


Amplitude 


| | | | | 


0 1 2 З 4 5 6 
Time (sec) 


Eng. de Controle e Automação – André Ferreira e Rejane de Barros Araújo 


E. MELHORANDO O ERRO DE REGIME PERMANENTE E 
A RESPOSTA TRANSITÓRIA USANDO PD+PI (PID) 


+ Exemplo 7: Parâmetros do PID 


К, ‚К 
2 1 2 
К Kıs + К, + Kas? к, (5 +з +É) 
G.(s] = К, +— +K s=—— Ç : =—— 


(s + 55,26)(5 + 0,5) _ ¿ ¿> (s? + 55,26s + 27,63) 


Gpip(s) = К 5 
К, = 5,62 


K 
F3 = 55,7762 K, = 313,37 


3 


K 
= = 27,63> K, = 155,28 


3 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO О ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Etapas de projeto: 


e 1) Avaliar o desempenho do sistema não compensado. 


• 2) Projetar o compensador LEAD para atender às especificações 
de resposta transitória (localização do polo, localização do zero e 


ganho de malha). 


• 3) Simular o sistema para verificar se as especificações foram 


E 


atendidas. 


° 4) Reprojetar o compensador LEAD, caso o desempenho transi 


náo esteja satisfatório. 


Eng. de Controle e Automacáo — André Ferreira e Rejane de Barros Araújo 


fis MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 


TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


“ Etapas de projeto: 


5) Verificar o desempenho do sistema com compensação LEAD 
com relação ao regime permanente para determinar quanto o 


mesmo deve ser melhorado. 


6) Projetar o compensador LAG para atender às especificações de 


regime permanente. 


7) Simular o sistema com compensador LEAD-LAG para verificar 


se todas as especificações foram atendidas. 


8) Reprojetar os controladores caso as especificações estejt) 


fora do que foi especificado. 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Ë MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Exemplo 8: Projete um compensador LEAD-LAG рага que o sistema 
abaixo opere com máximo sobressinal de 20% e uma redução pela 
metade no seu tempo de acomodação do sistema sem compensação. 


Além disto, o sistema compensado deve apresentar um erro em 


regime permanente para uma entrada rampa 10 vezes menor que o 


A as 
s(s + 6)(s + 10) 


sistema original. 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


“+ Exemplo 8: LGR do sistema não compensado. 


Root Locus 
30 | 
20 - 
10 + 
n 
£ 
B o ë А 
© 
D 
[t] 
E 
-10 - 
-20 - 
-30 l | | 
-40 -30 -20 -10 20 
Real Axis 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Ë MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


<+ Exemplo 8: LGR do sistema não compensado. 


M, (96) = 20% > ¿= 0,456 


Root Locus 
T T T Т Т Т Т 
| 
10 + 1 
8L System: g1 | 
Gain: 192 
Pole: -1.79 + 3.5i 
6 + Damping: 0.456 


Overshoot (96): 20 | 
Frequency (rad/sec): 3.93 


Imaginary Axis 
o N 


l l | | | | | | | 
-3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 
Real Axis 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Usando a função “rlocfind” 
>> [kli,pi]95rlocfindigil,-1.7943*3.5) 
кі = 


192.0875 


pi = 
-12.4130 


—3.7935 T+ 3.501231 
-1.7935 - 3.5012i 


** Usando a função "damp" para verificar o amortecimento 


>> damp (pl) 

Eigenvalue Damping Freq. (rad/s) 
-1.24e+001 1.00e+000 1.24e+001 
-1.79e+000 + 3.50e+000i 4.56e-001 3.93e+000 
-1.79e+000 - 3.50e+000i 4.56e-001 3.93e+000 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


Exemplo 8: 
+ Para K=192,09: 


° polos dominantes: 


o -1,7935 + j 3,5012 


• 3º polo: 
o-12,413 


Eng. de Controle e Automação – André Ferreira e Rejane de Barros Araújo 


| 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Simulação do sistema não compensado 


System: mf1 Step Response 
T T 


1.4 — Peak amplitude: 1.19 — 
Overshoot (%): 18.8 
Attime (sec): 1 


System: mf1 
Settling Time (sec): 2.19 


Amplitude 


| | | | | | 
1 1.5 2 2.5 3 3.5 
Time (sec) 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
| TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: 
t, 2,19 


£ а 


= t, =10958 


+ 


TO. — — ÓÓ GÀ 
sdes roues 7595 — 0,456 x 1,095 6 


Ogdes = Ondes |1 m p = 7,13 rad/s 


—À -—- 
Odes = SOndes > Odes = 3,65 


polos desejados > — 3,65 +] 7,13 


| 
Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Ë MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: arbitrar a posição do zero 


jo 


polos desejados > — 3,65 + j 7,13 


Ip.| > 6 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Contribuição angular dos polos 


As l 10)ls=-3,65+;7,13 = 48,31 


“(5)|5=-з,в5+у7 4з = 117,11" 


Ӯ, 0, = 165,429 


Eng. de Controle е Automação – André Ferreira е Rejane de Barros Araújo 


Ë MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


% Exemplo 8: Contribuição angular do polo do compensador LEAD 


=: -X 0, = —180° 
Go. = 180° -X 8, 


Oo. = 14,580 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Cálculo do polo do compensador LEAD 


7,13 
tan( 9.) = — <£ > 0,2601(p. — 3,65) = 7,13 
ds s+6 


s + 31,06 


GreaD (s) = 


p. = 31,06 


Eng. de 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


“ Exemplo 8: Definindo a FT com o compensador LEAD 


s+6 
s+ 31,06 


б кдр (5) = 


>> gz-t£r[i 6],convi[i 31.06],come([1 0) cony([1 6], [1 1D])))] 
Transfer function: 


а + 6 


374 + 47.06 s”3 + 557 s“2 + 1864 8 


Eng. de Controle e Automação – André Ferreira e Rejane de Barros Araújo 


Ё. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: LGR do sistema com o compensador LEAD 


Root Locus 
80 


+ + 


Imaginary Axis 
о 
ERR 


l l l l l l 
-120 -100 -80 -60 -40 -20 0 20 40 
Real Axis 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


% Exemplo 8: Usando a função “rlocfind” 


>> [k2,p2]7rlocfind(g2,-3.6543*7.13) 


-33.7602 
-36499 FE 7.13001 
-3.6499 — 7.13001 
-6.0000 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Usando a função “damp” para verificar o amortecimento 


>> damp [pel 


Eigenvalue Damping Freg. (rad/s) 


-3.38e+001 1.00e+000 3. 38e+001 
—3. 65e+000 + 7.13e+0001 4.56e-001 8.01ec000 
—3.65eTr000 — 7.13e+000i 4. 56-001 8.01er000 
—6.00e+000 1.002+000 6. 002+000 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO О ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Simulação do sistema com o compensador LEAD 


Step Response 
T T 


T 
Sistema não-compensado 
Sistema com compensador LEAD 


Amplitude 


| | | 


l 
1 1.5 2 2.5 3 3.5 
Time (sec) 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Ё. MELHORANDO О ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Simulação do sistema com o compensador LEAD 


t, 219 
Les == = = Е = L09585 
2 2 
System: mf2 System: mf1 se 
Peak amplitude: 1.19 1 — 3: Peak amplitude: 1.19 = | 
Overshoot (%): 19.3 Overshoot (%): 18.8 
At time (sec): 0.465 At time (sec): 1 
1.2 = 
| | System: mf1 
M, (%) = 20% = с = 0,456 Р | A o —9« ` sal сЕ 00 
| | | System: mf2 | 
|| Settling Time (sec): 1.07 | 
0.8 E | E 
3 E | 
= E 
£ 06 E - 
04 I : 
0.2 E E _ 
E Sistema não-compensado 
E Sistema com compensador LEAD 
0 | | | | | | | 
1 1.5 2 29 3 3.5 


Time (sec) 


Eng. de Controle e Automação – André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


« Exemplo 8: Sistema não compensado 


192 


nO SG 6) 10 


( 192 _ 192 


K, = lim |s ( ———— ——— 
v So l As(s + 6)(s + 10) 


s=0 


1 
encrampa (S°) = E 0,3125 


р 


Cães (оо) = 0,03125 


Eng. de Controle е Automação – André Ferreira e Rejane de Barros Araújo 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


% Exemplo 8: Sistema com o compensador LEAD 


Gira (S)Gyc (s) = ur) 


2166 


Grean (S) Gyc(s) = s(s+ 10)(s + 31,06) 


К li | ( EM ) TEL ENS 
БАР sol \5(5 + 10)(5+ 31,06) 10 x 31,06 à; 


— 0,1434 


€LEADra mpa (со) кз 
vLEAD 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Exemplo 8: Sistema com compensador LEAD-LAG: 
• Como a especificação do problema pede uma melhora de 10 
vezes no erro em regime permanente para uma entrada rampa, o 
compensador LAG deve melhorar o erro em regime permanente 


de um fator: 


@LEADrampa (со) 0,1434 
edes (00) 0,03125 


° Desta forma: 


FA 


Inl 


— 4,59 
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Е. MELHORANDO О ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Exemplo 8: Compensador LAG: 
• Assumindo: 
р. = 0,01 z, = 0,0459 
• Desta forma: 


s + 0,0459 


= ору 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Exemplo 8: Definindo a FT com o compensador LEAD-LAG: 


s+6 s + 0,0459 


ME des G ritos. reli 
5 + 31,06 rac (5) - 001 


Скар (5) = 


>> g3i=tf£f(conv([1 6], [1 0.0459]),convií[1 0.01], соп ([1 31.06], соқ ( [1 0], сопу ([1 6], [1 10]))))) 


Transfer function: 
з^2 + 6.046 s + 0.2754 


95 + 47.07:9º4 + 557.4 s^3 + 23865 ^72 + 28.645 
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Ë MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Exemplo 8: ІСК do sistema com compensador LEAD-LAG: 


Root Locus 
T 


Imaginary Axis 
N 
o o 

x 
"m 


2 1 | | | | i l | 
140 -120 -100 -80 -60 -40 -20 0 20 40 60 
Real Axis 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE E A RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


“+ Exemplo 8: Usando a função “rlocfind”: 


>> [k3,p3]=rlocfind(g3,-3.65+3*7.13) 


-33.7624 


-3.6307 + 7.1254i 
-6.0000 
-0.0461 
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Е. MELHORANDO О ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


> Exemplo 8: Usando a função “damp” para verificar o amortecimento: 


>> damp (p3) 

Eigenvalue Damping Freq. (rad/sz) 
—3 . 38e+001 1.00e+000 3. 38e+001 
—3.63e+000 + 7.13e+0001 4. 54e-001 8.00e4000 
—3.63e-r000 - 7.13e+0001 4.54e-001 8. 00e+000 
—6. 00e+000 1.00e+000 6. 00e+000 
—4.61e-002 1. 00e+000 4.61e-002 
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Е. MELHORANDO O ERRO DE REGIME PERMANENTE E А RESPOSTA 


TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Exemplo 8: Simulação do sistema com compensador LEAD-LAG 


Step Response 
T T 


T 
Sistema não-compensado 
Sistema com compensador LEAD 
Sistema com controlador LEAD-LAG 


Amplitude 


1.5 2 
Time (sec) 
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2.5 


3.5 


| 


Е. MELHORANDO О ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Exemplo 8: Simulação do sistema com compensador LEAD-LAG 


Step Response 


1 4 T T T T T T 
t 2,19 
1 2 E Ü.  ________ ч tsdes = 2 Е 2 БЛ tsaes а 1,095 5 
| System mf2 System: mf1 
| Settling Time (sec): 1.07 Settling Time (sec): 2.19 


li System: mf3 | 
|| Settling Time (sec): 1.04 | 
0.8 l | 


Amplitude 


0.6 


0.4 


Sistema não-compensado _| 
Sistema com compensador LEAD 


| | | Sistema com compensador LEAD-LAG 
0 i i | I | | I 
0 0.5 1 1.5 2 2.5 3 3.5 © 
Time (sec) 
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0.2 


Е. MELHORANDO О ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


+ Exemplo 8: Simulação do sistema com compensador LEAD-LAG 


System: mf3 System: mf 1 se 


Peak amplitude: 1.2 — — — Peak amplitude: 1.19 ~ 


Overshoot (96): 20.1 Overshoot (96): 18.8 
At time (sec): 0.467 At time (sec): 1 
1 2 YI ns 


System: mf2 | 
Peak amplitude: 1.19 | 


Overshoot (%): 19.3 ы шсш s= Mae 


At time (sec): 0.465 | 


0.8 


Amplitude 


Sistema não-compensado 


0.2 


Sistema com compensador LEAD 


Sistema com compensador LEAD-LAG 
| | | 


| ! | 
0 0.5 1 1.5 2 2.5 
Time (sec) 
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3 


3.5 


М„(%) = 20% > ¿= 0,456 


Е. MELHORANDO О ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Simulação da resposta do sistema a uma entrada rampa 


Step Response 


1600 

Rampa H 
1400 - Sistema náo-compensado 

Sistema com compensador LEAD 

Sistema com compensador LEAD-LAG 
1200 - L 


1000 - 


Ф 
Е 
E 800 
o 
£ 
< 

600 

400 E 

200 

0 | | 
0 500 1000 1500 


Time (sec) 
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| 


Е. MELHORANDO O ERRO DE REGIME PERMANENTE E А RESPOSTA 
TRANSITÓRIA USANDO COMPENSADOR LEAD-LAG 


* Exemplo 8: Simulação da resposta do sistema a uma entrada rampa 


Step Response 


15 T | | A 
—— Катра H 
14.5 Sistema não-compensado 
Sistema com compensador LEAD 
Wr Sistema com compensador LEAD-LAG 
13.5 
13r 
Ф 
Ë 
E 125 
о 
Е 
< 
12 Ł 
11.5 
11 + 
10.5 
10 l l l l l | | | 
10 10.5 11 11.5 12 12:5 13 13.5 14 14.5 15 


Тте (ѕес) 
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CARACTERÍSTICAS DOS COMPENSADORES 
EM CASCATA 


+ Pl e LAG: 


Function Compensator Transfer function Characteristics 
$ + Z; 
Improve steady-state error PI K Ea 1. Increases system type. 
2. Error becomes zero. 
3. Zero at —z, is small and negative. 
4. Active circuits are required to implement. 
Improve steady-state error Lag > is 1, Error is improved but not driven to zero. 
c 


2. Pole at —p, is small and negative. 

3. Zero at —z, is close to, and to the left of, the 
pole at —p.. 

4. Active circuits are not required to implement. 
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CARACTERÍSTICAS DOS COMPENSADORES 
EM CASCATA 


мтттонр 
má 


+ PDe LEAD: 


- 


Function Compensator Transfer function Characteristics 


Improve transient response PD K(s + ze) 1. Zero at —z, is selected to put design point on 
root locus. 


2. Active circuits are required to implement. 
3. Can cause noise and saturation; implement 
with rate feedback or with a pole (lead). 


А $ + Ze 
Improve transient response Lead K r. 1. Zero at —z, and pole at —p, are selected to put 
+ Pe design point on root locus. 
2. Pole at —p, is more negative than zero at —z,. 


3. Active circuits are not required to implement. 
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CARACTERÍSTICAS DOS COMPENSADORES 
EM CASCATA 


мтттонр 
má 


+ PID e LEAD-LAG: 


k SG + гем) 


Improve steady-state error and PID š 


transient response 


1. Lag zero at —zi,, and pole at origin improve 
steady-state error. 

2. Lead zero at —Zieag improves transient 
response. 

3. Lag zero at —zi4 is close to. and to the left of, 
the origin. 

4. Lead zero at —Zieag is selected to put design 
point on root locus. 


л 


‚ Active circuits required to implement. 


Can cause noise and saturation; implement 
with rate feedback or with an additional pole. 


EN 


js + гы: 
Improve steady-state error and Lag-lead (5 + zugs zie) 1. Lag pole at -Pay and lag zero at — Zag are used 


transient responsc (s + Pus)(s + Рем!) to improve steady-state error. 


2. Lead pole at —p,44 and lead zero at —Zjead are 
used to improve transient response. 


3. Lag pole at —p,,, is small and negative. 
4. Lag zero at —zjag is close to, and to the left of, 
lag pole at —py,. 


5. Lead zero at —Zieag and lead pole at —p,.,4 аге 
selected to put design point on root locus. 


6. Lead pole at —py¿,y is more negative than lead 
Zero at —Zicad- 


7. Activc circuits are not required to implement. 
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INSTITUTO FEDERAL DE EDUCAÇÃO, CIÊNCIA E TECNOLOGIA DO PARÁ Ë 
CURSO DE ENGENHARIA DE CONTROLE E AUTOMAÇÃO мумән 


Unidade VI: SISTEMAS DE CONTROLE NO ESPAÇO DE ESTADOS 


Disciplina: Engenharia de Controle | 


Professora: Prof.º Dr.? Rejane de Barros Araújo 


Introduçao E. 
Teoria de Controle i 


>” Clássico. 


% Moderno (a partir da década de 1960). 


CONTROLE CLÁSSICO CONTROLE MODERNO 


Sistemas monovariáveis Sistemas multivariáveis 
(SISO: single-input single-output) (MIMO: multiple-input multiple-output) 


Sistemas lineares e invariantes no tempo Sistemas lineares ou não-lineares, variantes ou 
não variantes no tempo 

Enfoque no domínio da frequência complexa Abordagem centrada no domínio do tempo 

Função de Transferência (modelo externo de | Espaço de Estados (modelo interno) 

entrada e saída) 
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Definições importantes - 


4 Estado: é o menor conjunto de variáveis (chamadas de variáveis de estado) com 
valores conhecidos em t = t,, que juntamente com o conhecimento dos valores dos sinais 
de entrada para t > tj, permitem determinar o comportamento do sistema para qualquer 


instante t > tọ. 


% Variáveis de Estado: são as grandezas cujo conjunto de valores determina o estado de 
um sistema dinámico. Se forem necessárias pelo menos n variáveis x4,x5,..,x, para 
descrever completamente o comportamento de um sistema dinámico, diz-se que essas n 


variáveis sáo um conjunto de variáveis de estado. 
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Definicoes importantes - 


4 As variáveis de estado não precisam ser necessariamente grandezas físicas 
mensuráveis ou observáveis, o que representa uma vantagem no método do espaco de 
estados. Na prática, porém, é conveniente escolher grandezas facilmente mensuráveis 
para serem as variáveis de estado do sistema, pois as leis de controle ótimo requerem a 


realimentação de todas as variáveis de estado com a ponderação adequada. 


** Vetor de Estados: é um vetor x constituído das n variáveis de estado necessárias para 


descrever completamente o comportamento do sistema dinámico. 


X1 

X2 
A = (=; 

Xn 
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Definicoes importantes - 


4$ Espaço de Estados: é o espaço n-dimensional cujos eixos coordenados consistem nos 


eixos x1,x2,...,Xn. Qualquer estado pode ser representado por um ponto no espaço de 


estados. 


4$ Equações no Espaço de Estados: modelagem de sistemas dinâmicos no espaço de 
estados 

“ variáveis de entrada 

$* variáveis de saída 


< variáveis de estado 
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Definicoes importantes - 


% A representação de um sistema no espaço de estados não é única, porém sabe-se 


que o número de variáveis de estado é o mesmo em qualquer representação. 


% Os integradores atuam nos sistemas dinâmicos contínuos no tempo сото 
“dispositivos de memória”, de tal forma que os sinais de saída desses integradores 
podem ser considerados como valores que definem o estado interno de um sistema 


dinâmico. 


% O número de integradores em um sistema dinâmico é também o número de variáveis 


de estado necessárias para definir completamente a dinâmica desse sistema. 
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Definiçoes importantes - 


+ т sinais de entrada: и, (t), и, (t), ..., u, (t). 
% m sinais de saída: у; (t), у (t), ..., уњ (t). 


% n variáveis de saída dos integradores (variáveis de estado): x, (t), x; (t), ..., x, (t). 


(0) = (66 t) > equação de estado 


y(t) = gx, u, t) > equação de saída 


* Linearizando em torno do ponto de operação 


x(t) = A(t)x(t) + B(t)u(t) > equação linear de estado 


y (t) = C(t)x(t) + D(t)u(t) > equação linear de saída 
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Definições importantes 


% Onde: 
A(t) matriz de estado 
B(t)> matriz de entrada 
C(t) matriz de saída 


D(t)> matriz de transmissão direta 
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Definicoes importantes - 


+ Supondo o sistema linear e invariante no tempo: 


x(t) = Ax(t) + Bu(t) 
y(t) = Cx(t) + Du(t) 


* Exemplo 1: circuito RLC série 


u(t) (0) с uc (t) 
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Definicoes importantes - 


* Exemplo 1: circuito RLC série (continuação). 


i(t) = c ee (1) u(t) = Ri(t) + 10 + uc(t) (1) 


duc(t) d*uc(t) 


(I) em (II) > u(t) = RC ДЕ + LC di + uc(t) 
2 
1с“ чс) + RC Tuc + uc(t) = u(t) сш j- : e puc culo) 


d^uc (t) Е R duc (t) 


ze uc(t) treu) 
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Definicoes importantes - 


* Exemplo 1: circuito RLC série (continuação). 


4 Fazendo 


xi (t) = uc(t) xa (1) = 200 
4 Tem-se: 
dx, ON dxo(t) | B il 


йй) O- Lx (D) + ul) 


dt dt Ic 


0 
1 
o + = u(t) 


КЕ E Е 0 


| (t) 
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Definicoes importantes Е. _ 
Ехегсісіоѕ I 


“+ 1) Para o circuito RLC série do exemplo anterior, monte a equação de estados supondo 


que x, (t) = u¿(t) e x>(t) = i(t). 


** 2) Monte as equações de saída para i(t) eu (t) do sistema anterior. 


** 3) Dado o circuito abaixo, ache sua representação no espaço de estados considerando 


x4 (t) = Volt), x2 (t) = ij (t) e y(t) = ig (t). 
V, (t) 


IO 
vo (*) V) = Velt) 
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Funções de transferência x Equações no Espaço de Estados e. 


= = Ax(t) + Bu(t) _ G(s) = Y(s) 
y(t) = Cx(t) + Du(t) U(s) 


* Aplicando a transformada de Laplace 


sX (s) — x(0) = AX (s) + BU(s) 
Y (s) = CX(s) + DU(s) 


% Para obter a função de transferência С (s), tem-se x(0) = 0: 


sX(s) = AX(s) + BU(s) 

(sI — A)X(s) = BU(s) 

(sI — A) 1 (sl — A)X(s) = (sI — A) 1BpU(s) 
X(s) = (sI — A) *BU(s) 
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Funções de transferência x Equações no Espaço de Estados e 


*& Substituindo X(s) em Y (s): 
Y(s) = C(sI — A) ! BU(s) + DU(s) —  Y(s)- [C(sl — A) !B + D]U(s) 


m = G(s) = [С(51 — A) 1 B + D] 


+ Definindo a matriz de transição de estados Ф = (sI — A) 1, tem-se: 


G(s) = [COB + D] G(s) = C(sI — AJ 1B * D 


polinômio(s) 


GS) = A 
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Funções de transferência x Equações no Espaço de Estados a. 


“ Observação 1: Inversa de uma matriz М. 


4 M é a matriz adjunta (matriz dos co-fatores transposta) 
“* Co-fator de a; jz>(-1)t det(sub — matriz retirando a linha i e a coluna j) 
* Para uma matriz de 2º ordem: 


T1141 ES = 1152 => | 


= М = 
—1П21 T1144 


M = | 
M21 ni» 
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Funções de transferência x Equações no Espaço de Estados B 


% Observação 2: O polinômio característico де G(s) é o determinante de (sI — A). Desta 


forma, os auto-valores de A correspondem aos polos de G(s). 


4 Auto-valores e auto-vetores de uma matriz А 


* Av= v,onde Zé о auto-valor e v é o auto-vetor de A. 


Àu — Av = 0 
(М — A)v = 0 


Para v = 0>det(AI—- А) = 0 
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Funções de transferência x Equações no Espaço de Estados - 


* Exemplo 2: Converter o sistema representado no espaço de estados para a forma de 


funcáo de transferéncia. 
«o -[; o + gua а= J в 


y(t) = [1,5 0,625]x(t) 


С = [1,5 0,625] D=0 
G(s) = [C(sI — 4у-18 + D] 21-4) = [6 [o с] 4 t 


(sI — A) 1 = Adj(sl — A) det(sl — А) = (s 4)5 +6 = s? + 45 + 6 


1 
det(sI—A) 
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Funções de transferência x Equações no Espaço de Estados - 


* Exemplo 2: continuação 


| fs —1,5 
AdjGI — А) =|) s+4 
* Desta forma 

-1,5 

=4 _ À 14 $+41 
Sa 52 +45 +6 

_ _ дү-1р — 
G(s) = CI AB == [1,5 0,625]|; - . ||| 
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Conversão de uma Função de Transferência para = _ 
Е$расо de Estados i 


* Como a escolha das variáveis de estado de um sistema dinámico não é única, existem 
inúmeras possibilidades de representação desse sistema na forma de equação de estados. 
Entretanto, independente da representação no espaço de estados, há uma única função de 


transferência que representa esse sistema. 


* Apesar da infinidade de representações possíveis para um sistema dinâmico no espaço 
de estados, algumas representações chamadas de FORMAS CANÔNICAS podem ser 
obtidas. 
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Conversão de uma Função de Transferência para E 
Espaço de Estados Е 


< Caso 1: 
Y(s) | 1 


H(s) = Ra 
(s) U(s) s? + as? +aĢzs + as 


(s? + a4,s? + azs + аз)Ү (з) = U (s) 
s3Y(s) + a4s?Y (s) + assY(s) + asY(s) = U (s) 


% Aplicando a anti-transformada de Laplace 


а?у(0) mp d’ y(t) Fa dy(t) 
dt3 1 dt? 2 dt 


+ azy(t) = u(t) 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


* Caso 1: colocando na forma de diagrama de blocos 


(X2) (хз) 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


% Caso 1: fazendo 
y 3 
X3=X2 
X2=X1 


X1- = 04X — 43X, = азХз tu 


« Arrumando as variáveis na equação de estados 


X1 — 04 — (02 — 3 X4 1 X4 
TEE 0 a y=[0 o 1il|xe 
x| lo 1 ollsl lo хз 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


* Caso 1: generalizando para um sistema de ordem k: 


Y (s) 1 
H(s) = — Vk k-1 k-2 
U(s) s*+ays“1+ а,5%? +. + Gk 
X= — а1Х — 09X2 — *** — Gk Xk + и 
X2=X1 
Хк-1=Хк-2 
Хк=Хк-1 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


* Caso 1: Então, tem-se 


x(t) = Ax(t) + Bu(t) f 
y(t) = Cx(t) + Du(t)' | 


— 42 єє. —Gk-1 — (Uk 1 
0 0 0 0 
1 0 0 в} 
: : 0 
0 1 0 - kx1 
0 1lixk pag 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


% Caso 2: 


Y(s) 5? + bzs + bs 


H(s) = o A cR M 
(s) U(s) s? + as? + а5 + as 


% Incluindo a variável intermediária Z(s), pode-se escrever: 


Y(s) Y(s)Z(s) 
U(s) Z(s)U(s) 


< De tal forma que: 


Y). diu" 
"im bjs“ + bs + bs 


Z(s) . 1 
U(s) sš+as2+a2s+as 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


< Caso 2: 


U(s) Z(s) Y (s) 


1 b,s? + bas + bs 
s? + as? + ass + as 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


“* Caso 2: expandindo o diagrama 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


4 Caso 2: eliminando os derivadores 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


«+ Caso 2: fazendo 
у= b4x4 T box; + b3X3 


X3=X2 
X2=X1 
Ki=— Xi — doX» = d3Xs FU 


* Arrumando as variáveis na equação de estados 


X —a, —az —@3][Х1 1 X1 
X2 = | 1 0 0 х2|+ JO[u у= [bi b; b3] E 
x| lo 31 oll lo хз 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


* Caso 2: Generalizando para um sistema de ordem k: 


H(s) Y(s) bis 1 + b,s*7? +... + by 
5) = —— = AOASKÁKÁ 

U(s) st+tasti+tast2+.-.+ay 
Xi- — d4X4 = бай» — s — Qi TU 

X22X4 
y = bixi + Бх +- + рех 
Xk-1-Xk-2 
Хк=Хк-1 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


“ Caso 2: Forma de equação de estados 


x(t) = A¿x(t) + B-u(t) | 
m = Cox(t) + D,u(t) PM 


241 -а: Gk—1 ак 
1 0 0 0 1 
0 1 0 0 0 
: Ў В, = 
0 0 1 0 kxk 0 Кх1 


Ce = [b1 b, .. bk|ixk 


essa forma particular dá-se o nome de forma canônica do controlador. 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


< Outras formas canônicas 


% Forma canônica do observador 


Ao = Ac 
B = M 
Co = Bo 
D, = 0 


4 Forma canónica de Jordan: é obtida fazendo-se a expansão em frações parciais da 
função de transferência do sistema, com a matriz de estados sendo uma matriz 


diagonal com os polos desta função de transferência. 
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Conversão de uma Função de Transferência para E 
Espaço de Estados i 


* Exercício: Dado o sistema G(s), obtenha a sua representação no espaço de estados 


na forma canônica de controlador. 


s+2 


G(s) = 216545 
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Controlabilidade e Observabilidade Е. _ 


* Considere a equação de estados com n estados, r entradas e m saídas: 


x(t) = Ax(t) + Bu(t) 
y(t) = Cx(t) + Du(t) 


* com A ERP В є В" e C e Rmzn. 
* O раг (А,В) é dito ser controlável se, para qualquer estado inicial xy e para qualquer 
estado final x,, existir uma entrada u(t) que transfere o estado de x, para x, em tempo 


finito. 
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Controlabilidade e Observabilidade - 


“ A definição requer apenas que se possa mover qualquer estado inicial no espaço de 
estados para qualquer estado final em tempo finito, não havendo restrição quanto à 
trajetória a ser seguida ou à magnitude do sinal de entrada. É importante notar que a 
equação de saída não tem influencia sobre a controlabilidade do sistema. 

** A matriz de controlabilidade de um sistema é dada por: 


TO =[В AB А2В .. A"-1g] 


“ Diz-se que um sistema é controlável se a matriz GO tiver rank (número de linhas ou 


colunas linearmente independentes) completo de linhas. 
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Controlabilidade e Observabilidade E 


** O par (A,C) é dito ser observável se para qualquer estado inicial xy existir um tempo t; 
tal que o conhecimento da entrada u e da saída y no intervalo [0,t,] seja suficiente para 
se determinar de maneira Única xo. 


* A matriz de observabilidade de um sistema é dada por: 
G 
oB-| (4 
САт-1 


“ Diz-se que um sistema é observável se a matriz OB tiver rank completo de colunas. 
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Controlabilidade e Observabilidade Be 


* No MATLAB, os comandos “ctrb” е “obsv” retornam, respectivamente, as matrizes de 


controlabilidade C e de observabilidade O. Utilizando o comando “rank” é possível 


determinar se um sistema é controlável e/ou observável. 


“ Exemplo 3: Considere o seguinte sistema dinâmico mostrado abaixo. Calcule as 
matrizes de controlabilidade e observabilidade e determine se o mesmo é ou não 


controlável e observável: 


2 0 0 1 
x(t) = [ 1 ] x(t) + [0| u(t) 
4 0 3 0 


у@)=[2 —1 Olx(t) 
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Sistemas Equivalentes 


4. Considere o sistema 


x(t) = Ax(t) + Bu(t) 
y(t) = Cx(t) + Du(t) 


* Seja х = Рх, сот P não singular. Então, pode-se escrever o seguinte sistema 


equivalente: 
X(t) = Ax(t) + Bu(t) 
y(t) = Cx(t) + Du(t) 


< comA-PAP,B-PB,C-CP 'eD-D. 
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Sistemas Equivalentes E 


* Pode-se mostrar que as propriedades de estabilidade, controlabilidade е 
observabilidade são preservadas pela transformação de equivalência e que as matrizes de 


controlabilidade e de observabilidade se relacionam da seguinte forma: 
CO=PTO 


OB = OB P-1 
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Realimentação de estados Е. _ 


% A ideia básica da realimentação de estados consiste em alocar os polos de malha 


fechada (autovalores da matriz dinâmica), modificando, assim, a dinâmica do sistema. 


“ Dada uma representação em variáveis de estado de um sistema 


x(t) = Ax(t) + Bu(t) 
y(t) = Cx(t) + Du(t) 


ração — André Ferreira e Rejane de Barros Araújo 


Realimentação de estados Е. _ 


+ Usando realimentacáo de estado, cada variável de estado é multiplicada por um ganho e 


realimentada para o terminal de entrada, ou seja: 

u(t) = —Kx(t) + r(t) 
4 onde К = [k, k, .. kn] é o vetor de ganhos de realimentação. 
+ Assim, temos: 


x(t) = Ax(t) + Bu(t) 
u(t) = Kx(t) + r(t) 


x(t) = (A — BK)x(t)* Br(t) 
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— x(t) = Ax(t) + B(-Kx(t) + r(t)) 


Realimentação de estados Е. _ 


“+ Devemos ter acesso a todos os estados do sistema. 


* Teorema: Se (A,B,C,d) for controlável, usando u(t) = -Kx( + r(t podemos escolher 


bitrariamente os autovalores de (A - ВК). 
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Realimentação de estados Be 


* Método da substituição direta para a Determinação da Matriz de Ganhos K. 


% Se o sistema for de ordem reduzida (n < 3), a substituição direta da matriz K no 
polinômio característico desejado poderá ser mais simples; 
% Por exemplo, se п = 3, então escreve-se a matriz de ganho K de realimentação de 
estados como: 
K=[k k2 k3] 
% Substitua a matriz K no polinômio característico desejado |SI — A + BK| e iguala a 


(s + 14)Cs + и»)(5 + из), isto é: 


IST А + BK| = (s + m)(s + uo)(s + из) 


Eng. de Controle e Automação — André Ferreira e Rejane de Barros Araújo 


Realimentacao de estados ie 


% Exemplo 4.1: Dado o sistema em espaço de estados, usando u(t) = —Kx(t) + r(t). 
Determine K para que os polos desejados de malha fechada estejam localizados em s = 


-2 +j4es = –10. 


0 1 0 0 
х= |0 0 1|x-4l0l[u 
—1 -5 -6 1 


** simulacáo matlab 

“А = [0 1 0;0 0 1;-1 -5 -6]; 
$ В = [00 1]; 

4$ J = [-24j*4 -2-j*4 -10]; 
* K = place(A,B,J) 
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Realimentação de estados Е. _ 


* Fórmula de Ackermann para Determinação da Matriz de Ganhos К. 
* 1) Formar A(s) = s" +as””!+...+ta,.4s + a, com polos desejados. 


% 2) Calcular K da seguinte forma: 


K=[0 0 +... 1]U tada) 


О = [В AB А2В .. pn1p] 
q. (A) = А" + аА"! +... + anl 
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Realimentacao de estados - 


4 Exemplo 4.2: Dado o sistema em espaço de estados, usando u(t) = —Kx(t) + r(t). 


Determine K para que os autovalores do sistema sejam -1 e -2. 


+=, alere 


y=[1 1]x 


simulação matlab 

“А = [0 1 0;0 0 1;-1 -5 -6]; 
$ В = [00 1]; 

4$ J = [-24j*4 -2-j*4 -10]; 
% K = acker(A,B,J) 
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Observadores de estados E 


* O observador de estados consiste em um mecanismo (algoritmo) para estimação dos 


estados da planta. É uma solução útil quando os estados reais da planta não estão 


accessíveis, situação muito comum na prática . 


% Seja 


x(t) 2 Ax(t) + Bu(t) 
y(t) = Cx(t) 


* conhecendo A, B e C e as medidas de y(t) e u(t), constrói-se o estimador: 


£(t) = ARCE) + L(y(t) – $(0) + Bu(t) 
$(t) = CX(t) 
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% onde L = [L 


l2 


Observadores de estados Е. _ 
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Observadores de estados Е. _ 
Erro de estimação i 


** O erro entre x e 2, conhecido como erro de estimação (ou erro de observação), é dado 
por; X(t) = x(t) — X(t), derivando-se, temos: 

X(t) = x(t) — X(t) 

X(t) = [Ax(t) + Bu(t)] — [A2(t) + LOC) — Ў(@)) + Bu(t)] 


X(t) = A(x(t) — 2(t)) — LC(x(t) — £(t)) 
X(t) = (A — LO)X(t) 


* Para que lim Z(t) = 0 é necessário que os autovalores de (A – LC) sejam estáveis. 
—00 


* Teorema: Se (A,B,C,D) for observável, então um estimador de estado assintótico com 
isquer autovalores pode ser construído. 
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Observadores de estados Е. _ 


* Fórmula de Ackermann рага Determinação da Matriz de Ganhos do observador L. 
* 1) Formar A(s) = s" + a,s* ^! +- + а„_15 + a, com polos desejados. 


4 2) Calcular L da seguinte forma: 
Е= W o 0 « 4 


onde 


V=[C CA C4? .. cart 
q,(4) = А" + аА" +... + anl 
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Observadores de estados 


com autovalores -3 e -3. 


+=, zerie 


y=[1 1]x 
% simulação matlab 


c(t) 
de E > A WM E + А u(t) 


«c mo 


+ Exemplo 5: Dado o sistema em espaço de estados, projetar um observador de estados 
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Realimentacao de estados observados - 


“ Para a realimentação de estados é necessário que todos os estados reais da planta, x(t), 
sejam mensuráveis. Quando isto nào ocorre, há a necessidade de construir um observador 


de estados. Neste caso a realimentação é feita a partir dos estados estimados, X(t): 


u(t) X(t) y(t) 
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Realimentacao de estados observados - 


x(t) = Ax(t) + Bu(t) 
y(t) = Cx(t) 


% Se X(t) é uma estimativa de x(t) , então, na realimentação de estados utiliza-se: 


u(t) = KX(t) + r(t) 


“ Porém, se a realimentação é feita a partir dos estados estimados, a dinámica do 


estimador precisa ser considerada. Desta forma, tem-se: 
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Realimentação de estados observados Е. _ 


10 = 
ES E lic (А— BK — ic) "d “ B uit) 


* Exemplo 6: Utilizar a realimentação de estados em conjunto com o observador de 


estados projetados anteriormente. 
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Seguidores de referëncia ou servosistemas - 


“ Para que um sistema descrito por variáveis de estado possa, além de possuir a dinâmica 
desejada (garantida pela alocação de polos por realimentação de estado), seguir uma 


determinada entrada, com erro zero, usamos o princípio do modelo interno. 


+ Considere entradas de referência descritas por equações do tipo: 
Xy = Ap Xy 
т = Cyxy 
* com condições iniciais desconhecidas. Um modelo equivalente para entradas de 


referência é: 
ro) apar D+ + at + оог = 0 


g onde 7) é a n-ésima derivada de r(t). 
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Seguidores de referência ou servosistemas E. 


* Exemplo 7: 
* Degrau unitário: r(t) 21; t20; r=0: 


x, = 0x 
* Escolhendo x, = r, temos: ра - i 
== T 


* Rampa ипіїагіа: r(t) =, t> 0; E : 
EXE | 
XH = E 
ku Escolhendo | "1 '. temos: 4 2 a 
Arg = T 


r-[ 0] > 
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Seguidores de referëncia ou servosistemas 
Princípio do modelo interno para referência do tipo degrau 


NETITUTOFED 


+ Considere o sistema 


x(t) = Ax(t) + Bu(t) 
y(t) = Cx(t) 


*& Definimos o erro de rastreamento como: e(t) = y(t) — nt). 
* Das características do sinal de referência, temos: é(t) = y(t) — r(t)= y(t) = Cx(t). 


4. Definimos novas variáveis de estado como: 


e = cr 


goa x . 
| =" entao temos: l "T 


w = 
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Seguidores de referëncia ou servosistemas 
Princípio do modelo interno para referëncia do tipo degrau 


NETITUTOFED 


* Na forma matricial 
(= [р alEl+ isl” 


* Se o sistema for controlável, então, existe uma lei de controle por realimentação de 
estado, da forma; w=k,e+k,z, tal que os polos do sistema aumentado podem ser 


posicionados arbitrariamente. 


* Desde que os polos do sistema aumentado sejam alocados na regiáo de estabilidade, o 
erro de rastreamento será estável. Assim, o objetivo de rastreamento assintótico com erro em 


regime nulo será alcançado. Ou seja, a resposta do sistema abaixo é assintoticamente 
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Seguidores de referëncia ou servosistemas 
Princípio do modelo interno para referëncia do tipo degrau 


El E E A KM И 


NETITUTOFED 


** O sinal de controle u(t) é obtido pela expressão: 


t 


t 
ui) = | м(т)ат = ks | е(т)ат + k;x(t) 
0 0 


t 


X = Ax + Bu завк | «(dt + kax) 
0 
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Seguidores de referência ou servosistemas 
Princípio do modelo interno para referência do tipo degrau 


+ Definindo: x,,4(t) = f е(т)ат, temos: X, 1 (t) = e(t) = r(t) – y(t) = r(t) — Cx(t). 


** Logo: 


Је lelt i 
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Seguidores de referëncia ou servosistemas E. 


* Exemplo 8: Considere o sistema em espaço de estados. Projetar um controlador para 


que o sistema tenha erro nulo para uma entrada em degrau. 


йи № + [5] 


y=[1  0]x 
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ENGENHARIA DE CONTROLE | 


UNIDADE V: RESPOSTA EM FREQUÊNCIA DE 
SISTEMAS LINEARES 


Engenharia de Controle e Automação 
Prof. André Ferreira/Rejane de Barros Araújo 


ш SISTEMAS DE CONTROLE 


NEED. |  À]— 


« Histórico: 
*$ Nyquist (1932) 
* Critério de Estabilidade de Nyquist 


Harry Nyquist 1.5 


Harry Nyquist (1889-1976) 


February 7, 1889 (s) = 


Stora Kil, Nilsby, Vármland, a 4x41 
Sweden -1.5 , 


April 4, 1976 (aged 87) -0.4 -0.2 00 02 04 06 08 10 12 


Harlingen, Texas, USA Бе 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


a SISTEMAS DE CONTROLE 


% Histórico: 
* Bode (1945) 
*• Diagramas de Bode 


Hendrik Wade Bode M t —High Pass 
= Бы 
2 
o — Bode Pole 
oO 
2 
= 
c 
ob 
m 
= 1 10 100 1000 10000 
Frequency 
T " А —-High Pass 
Ф 
Ф 
юе 
— воде Pole 
- б 45 
Hendrik Wade Bode = 
Ф 
Вогп 24 December 1905 E: 
Madison, Wisconsin a 5 
Died 21 June 1982 (aged 76) 1 10 100 1000 10000 


Cambridge, Massachusetts 


Frequency 
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SISTEMAS DE CONTROLE 


< Histórico: 


“ Evans (1948) 
% LGR 


Walter R. Evans 


Root Locus of (5+2)/(59+352+5з+1) 


Br т T ша; Т q -— - — 


Imaginary Axis 


Born January 15, 1920 
Died July 10, 1999 (aged 79) 
Residence United States 
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E RESPOSTA EM FREQUÊNCIA DE SISTEMAS LINEARES 


Шы». — 


x(1) 


Output 
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8 RESPOSTA EM FREQUÊNCIA DE SISTEMAS LINEARES 


BED — . 


M,(@)Zó,(e)= М (о)м (@)Z]ó (0)+ plo)! 


Resposta da magnitude no 
domínio da frequência 


$(0)- ф,(®)—ф(®) —— dominio da trequência ” 
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8 RESPOSTA EM FREQUÊNCIA DE SISTEMAS LINEARES 


G. — 1] I ] |.) 


* Resposta da magnitude no domínio da frequência: 
“ Relação entre as magnitudes das senoides de saída e 


de entrada. 
“ Resposta da fase no domínio da frequência: 


« Diferença entre os ángulos de fase das senoides de 


saída e de entrada. 
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8 RESPOSTA EM FREQUÊNCIA DE SISTEMAS LINEARES 


Do .  h | | . 


4 Essas respostas são função da frequência e valem apenas para a 


resposta em regime estacionário das senoides aplicadas ao sistema. 


% Pode-se mostrar que a resposta em frequência de um sistema G(s) é 


dada por: 
As+ Bo 
c(s) = G(s) 
(5? + a?) 
— — | C=C) 
RG)- As+ Bo 
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F- RESPOSTA GRÁFICA DA RESPOSTA EM FREQUÊNCIA 


BED 


* A FT senoidal é uma função complexa da frequência w, sendo 


caracterizada pelo seu módulo e seu ángulo de fase. 
% Existem duas maneiras para representar graficamente a resposta em 
frequência de G(jw) = М. (jo)zdg (jo): 
* Diagrama de Bode: com gráficos de magnitude e fase separados 


(gráfico logarítmico; dB = 20 log M). 


*• Diagrama de Nyquist: gráfico polar (representação fasorial). 
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= _ RESPOSTA GRÁFICA DA RESPOSTA EM FREQUÉNCIA 


NS EXEMPLO 


% Magnitude: 


Е 1 — 1 1 [(2-joY 2-jo 
= ттт 2+ 2—jo)] w?-4 
% Fase: 
IG(jw)| = M(o) : 20log M(w) = 201 | > J 
w)| = M(w) = = O w) = og | ——— 
Vw? + 4 i x Vw? + 4 


ZG(jo) = $(jo) = —atan (5) 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


F- RESPOSTA GRÁFICA DA RESPOSTA EM FREQUÊNCIA 


BED . > | . | 5)  h. 


% Magnitude e Fase em gráficos separados: 


Bode Diagram 
T 
I 


20log(M) (dB) (dB) 


Phase (degrees) (deg) 


Frequency (rad/s) 
jarros Araújo 


F- RESPOSTA GRÁFICA DA RESPOSTA EM FREQUÊNCIA 


EE í | 


% Gráfico Polar (magnitude e fase no mesmo gráfico): 


Re G 


Ф = 0 = 
O [0 005 Ol 015 02 025 03 035 0. Aa a анн 
0.1 r/s 
-0.05 0.2 r/s 
0.3 r/s 
-0.10 0.4 r/s 
0.5 r/s 


Im G 


-0.15 


-0.20 


-0.25 
Note: r/s = rad/s 
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F- RESPOSTA GRÁFICA DA RESPOSTA EM FREQUÊNCIA 


BED =  h3 h hA /. |. 


% Esses gráficos podem ser obtidos com o auxílio do MATLAB. 


4 Definindo inicialmente a FT: 


>> g-tfil.polyií-z)]] 
Transfer function: 
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4$ Diagrama de Bode 


* >> bode(g) 


Bode Diagram 


* >> grid 


Magnitude (dB) 


Phase (deg) 


Frequency (rad/sec) 
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« Diagrama de Nyquist 


$ >> nyquist(g) 


Nyquist Diagram 


* >> grid 


Imaginary Axis 


Real Axis 
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E APROXIMAÇÕES ASSINTÓTICAS: 


Dom» DIAGRAMA DE BODE 


*• Os esboços dos diagramas de Bode podem ser simplificados através da 


sua aproximação por uma sequência de linhas retas. 


О pl p) (ss p.) 


Kt zz) mul 
dl р 6+ p.)- [65 p. 


© 
- 
Ы 


s>j0 
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E APROXIMAÇÕES ASSINTÓTICAS: 


һи 


ENTE: DIAGRAMA DE BODE 


“ Utilizando a representação em dB (dB = 20 log M), tem-se: 


20 log (G(¡w) = 20 log K + 20 log |(5 + z, ) + 20 log ((s + z, ) 
Te—20 log |s” |- 20 log (s 4 p=": 


s— jo 


€ Fazendo-se o mesmo procedimento para a fase, somam-se os ângulos 
de fase dos termos relativos aos zeros e subtraem-se os ángulos de 


fase dos termos referentes aos polos. 
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— DIAGRAMA DE BODE PARA G(s) = 5 + а 


G(s)=s+a > G(o) = jo +a = a (1 +1) 


€ Ет baixas frequências: 


о — 0 


э 20106 M = 20loga 
G(jo) = a > G(jw) = az0? | 


ф = 0° 
* Em altas frequências: 


п зей > 20 log М = 20log w 
j ф = 900 
G(jo) = a (12) = G(jo) = jw = 090% 
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E APROXIMAÇÕES ASSINTÓTICAS: 


EE DIAGRAMA DE BODE 


% Plotando o gráfico de dB x log w, tem-se: 


y = 20x 
“ onde: 
y = 20logM 
x = logo 


4$ Desta forma, pode-se dizer que para cada década (intervalo entre duas 
frequências, onde a última é dez vezes maior que a primeira) a 


magnitude sofre um aumento de 20 dB. 
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Ё. _ APROXIMAÇÕES ASSINTÓTICAS: 


ми 


EE DIAGRAMA DE BODE 


% Na frequência de quebra о = a, tem-se: 


G(jo)=a+ja >  G(jo) = Mz459 


« Assim sendo, o diagrama de Bode de G(s) = s + a pode ser desenhado 


utilizando as aproximacóes assintóticas: 
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Е. _ APROXIMAÇÕES ASSINTÓTICAS: 


MER DIAGRAMA DE BODE 


42 + 20 loga 
36 + 20 loga 
30 + 20 loga 
> 24 + 20 loga 


& 18+20loga 
12 + 20 log a 
v d EUM 
wga AU A Billi! 
ИШ PA 1 Hie 8081111 


0.01а 0.1а 


(а) 


90 , 

s ШШ ШШШ LEAL LL ELIT 
# c | ү LA ud |р 
Pg К ШШ E EHI 
3 Л ll eew LII 
E ШШ ИШ ү ИШ 


A ИШ ИШ LllHH 


fenis (rad/s) 
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= _ APROXIMAÇÕES ASSINTÓTICAS: 


EE DIAGRAMA DE BODE 


% Comparação entre a aproximação assintótica e resposta em frequência 


normalizada para G(s) =s +a 


Frequency 20 log z (dB) Phase (degrees) 

(rad/s) Asymptotic Actual Asymptotic Actual 
0.01 0 0.00 0.00 0.57 
0.02 0 0.00 0.00 1.15 
0.04 0 0.01 0.00 2.29 
0.06 0 0.02 0.00 3.43 
0.08 0 0.03 0.00 4.57 
0.1 0 0.04 0.00 5.71 
0.2 0 0.17 13.55 11.31 
0.4 0 0.64 27.09 21.80 
0.6 0 1.34 35.02 30.96 
0.8 0 2.15 40.64 38.66 
1 0 3.01 45.00 45.00 
2 6 6.99 58.55 63.43 
4 12 12.30 72.09 75.96 
6 15.56 15.68 80.02 80.54 
8 18 18.13 85.64 82.87 
10 20 20.04 90.00 84.29 
20 26.02 26.03 90.00 87.14 
32.04 32.04 90.00 88.57 

60 35.56 35.56 90.00 89.05 
80 38.06 38.06 90.00 89.28 


100 40 40.00 90.00 89.43 
—— n n n F [n -[] —əəƏ ç I———— ira/Rejane de Barros Araújo 


= _ APROXIMAÇÕES ASSINTÓTICAS: 


EE DIAGRAMA DE BODE 


% Comparação entre a aproximação assintótica e resposta em frequência 


normalizada para G(s) =s +a 


| I И 
| jl || И 
|... ] LLLLUL И 
a j LAU 
_ | | JI м ү 
_ Л || 
|] A ilt 


20 log(M/a) 
© 


ИУ 


F TEqUENCY (rad/s) 


Asymptotic approximation I 
10 
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E APROXIMAÇÕES ASSINTÓTICAS: 


m DIAGRAMA DE BODE 


% Comparação entre a aproximação assintótica e resposta em frequência 


normalizada para G(s) = ѕ +a 


o PTA 
A IH UA! 
= ue e NI 


ГИШГИШ 
V ПШ 


неве (rad/s) 


EN 
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= _ APROXIMAÇÕES ASSINTÓTICAS: 


25 4 DIAGRAMA DE BODE 


% Utilizando o Matlab 


>> g-tfi[l 10],1) 


Transfer function: 


s + 10 


>> роде [q] 
>> grid 


Magnitude (dB) 


Phase (deg) 


Bode Diagram 


Frequency (rad/sec) 
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F- DIAGRAMA DE BODE PARA G(s) = 1/(s +a) 


MD = 3 1 O —— 


1 1 1 
(s) stu Е (o) Jó + a a (E +1) 


* Em baixas frequências: 


& э Ü = lo 
Peg =ana (+) - 20 log a 


E o Ü 
G(jo) = => бо) = 220 ф= 0 
* Em altas frequências: 
Q) > a 1 
ШЕ > 20logM = 201og— 
G (jo) = > G(jw) =— = —¿ — 909 ф = —909 


17 ü 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


E APROXIMAÇÕES ASSINTÓTICAS: 


ми 


ШЕ» DIAGRAMA DE BODE 


% Na frequência de quebra w = a, tem-se: 


G(jw) = >  G(jo)- Mz —45? 


a + ja 


« Assim sendo, o diagrama de Bode de G(s) = s + a pode ser desenhado 


utilizando as aproximações assintóticas (gráfico normalizado): 
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E APROXIMAÇÕES ASSINTÓTICAS: 


| 4 DIAGRAMA DE BODE 


l 

l —6 dB/octave = || 
EH — 220 dB/decade Hi 

l 

Ë 


_ СОЕ HINL 
40 LH L a ШЕШ 
0080800008808 A 08080141 


0.0] 0.1 l 10 100 
Frequency (rad/s) 
a 


TIENE T TTE T LIII 
G (s) HHN LO lil 
T ES LILHI JII 


e L HINL ЕШШ TT Mi 
ps TICA 245º/decade 
—HEHHBE T JlIHHƏIP 80111111 


0.01 0.1 | 10 100 
Frequency (rad/s) 
а 


20 logMa 
© 


90 
45 
0 


Phase (degrees) 
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= _ APROXIMAÇÕES ASSINTÓTICAS: 


|. 8)» DIAGRAMA DE BODE 


% Utilizando o Matlab 


Bode Diagram 


>> g-tfíl,[1 10]) 


@ 
Transfer function: ср 
- š 
з + 10 
>> роде {gi 
>> grid 


Phase (deg) 


Frequency (rad/sec) 
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E. APROXIMAÇÕES ASSINTÓTICAS: 


| 385 DIAGRAMA DE BODE PARA INTEGRADOR 


% Usando o MATLAB: 


Bode Diagram 


>> g-tfíl,[1 0]) 


s 
g 
Transfer function: 5 
i > 
3 
>> bode (q) 2 
Ko) 
>> grid 2 
& 


Frequency (rad/sec) 


% OBS.: Cruza 0 dB para o = 1 rad/s 
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= _ APROXIMAÇÕES ASSINTÓTICAS: 


ШЫП» DIAGRAMA DE BODE PARA DERIVADOR 


4$» Usando o MATLAB: 


Bode Diagram 


>> g=tf([1 0],1) 


T 
Transfer function: š 
3 
>> bode (g) 
>> grid 


Phase (deg) 


Frequency (rad/sec) 
% OBS.: Cruza 0 dB para w = 1 rad/s 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


E DIAGRAMA DE BODE PARA FATORES 


ME DE 1º ORDEM 


“ Exemplo: 


% Frequências de quebra: 1 rad/s, 2 rad/s e З rad/s. 
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- DIAGRAMA DE BODE PARA FATORES 


BENE DE 1º ORDEM 


к 1 
2 13 


Tn P 


* O gráfico deve começar pelo menos uma década abaixo da menor 


G(s) = 


frequéncia de quebra e terminar pelo menos uma década acima da 


maior frequência de quebra: 
0,1 < o < 100 rad/s 


* O efeito de K é mover a curva de magnitude para cima (TK) ou para 
baixo (LK) de um valor 20 log K, com o mesmo não apresentando efeito 


na curva de fase. 
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E DIAGRAMA DE BODE PARA FATORES 


|o m» DE 1º ORDEM 


* Em о = 0,1 rad/s (baixas frequências) aproximar (2+ 1) com s— 0 е 


рага о termo s по denominador usar 0,1. 


|&(jo)|o=03 = —B їз 


K=1 > 0=01 > |G(jo)|ap = 2010g15 = 23,52 dB 
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E DIAGRAMA DE BODE PARA FATORES 


к» DE 1º ORDEM 


Ж и 
2 \3 


ез P 


% Cálculo da magnitude utilizando as aproximações assintóticas: 


G(s) = 


Frequency (rad/s) 
0.1 (Start: 1 (Start: 2 (Start: 3 (Start: 
Description Pole at 0) Pole at —1) Pole at –2) Zero at —3) 
Pole at 0 —20 —20 —20 —20 
Pole at —1 0 —20 —20 —20 
Pole at —2 0 0 —20 —20 
Zero at —3 0 0 0 20 


Total slope (dB/dec) —20 —40 —60 —40 
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E DIAGRAMA DE BODE PARA FATORES 


NE» DE 1º ORDEM 


3 К | 
2 \3 


n P 


% Gráfico da magnitude utilizando as aproximações assintóticas (efeito 


G(s)- 


individual dos polos e zeros): 


40 
, -HE 
173 — = 1 
E NM E - 
Ii 
-20 dB/dec ($6 dB/oct). ALTRES 
do 
E 0) 


Frequency (rad/s) 


LIL 


| 


20 log M 
кә 
© 
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E DIAGRAMA DE BODE PARA FATORES 
EPA DE 12 ORDEM 


`k ы 
2 \3 


EN 


Gráfico da magnitude utilizando as aproximações assintóticas (resultado 
final): 


G(s)- 


% 


—20 dB/dec (—6 dB/oct) 


aeneo | | II 
<y ҢА Hit 
| 
— 40 dB/dec (-12dB/oct) 
— HL ES 


TAI 
A A SESS 
Lit Lit 1 1 TERRI 


100 


Frequency (rad/s) 
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E DIAGRAMA DE BODE PARA FATORES 


к. DE 1º ORDEM 


K ы 


% Cálculo da fase utilizando as aproximações assintóticas: 


Frequency (rad/s) 
0.1 (Start: 0.2 (Start: 0.3 (Start: 0 (End: 20 (End: 30 (End: 
Description Pole at —1) Pole at —2) Pole at —3) Pole at —1) Pole at —2) Zero at —3) 
Pole at —1 —45 —45 —45 0 
Pole at —2 —45 —45 —45 0 
Zero at —3 45 45 45 0 
Total slope (deg/dec) —45 —90 —45 0 45 0 
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E DIAGRAMA DE BODE PARA FATORES 


кз DE 1º ORDEM 


3 К | 
2 \3 


n P 


4 Gráfico da magnitude utilizando as aproximações assintóticas (efeito 


G(s)- 


individual dos polos e zeros): 


LL 
iii >a 
T L IILI 
JL 1 LII 
l 10 


Frequency (rad/s) 
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E DIAGRAMA DE BODE PARA FATORES 


BED DE 1º ORDEM 


`k ы 
2 \3 


EN 


“ Gráfico da magnitude utilizando as aproximações assintóticas (resultado 


G(s) = 


final): 


== ЕЕЕ MM] 
ҮШ _ | 


Phase (degrees) 


Frequency (rad/s) 
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E DIAGRAMA DE BODE PARA FATORES 


ME DE 1º ORDEM 


ы 


+ Utilizando o Matlab 
>> g-ti([1 3],poly([O -1 -2])) 
Transfer function: 
3 + 3 


3^3 + 3 а^ + 2 8 


>> bode (g) 
>> grid 
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йы Cas kao hamana cd ue 
-—--кюк---++---+---+----+----) 


_ L 


ORDEM 


Bode Diagram 


1 
- 
E 
EI 

1 
E 

t 
_ 

1 
E 

1 

1 
q 

I 

І 

I 
m 

I 

1 

1 

1 

1 

1 

1 


ОЕ 1 
G(s) 


EN 


Es ium n om m ab m c. 


DIAGRAMA DE BODE PARA FATORES 


1 

1 

1 

1 

1 

1 

I 
[Тө] о 
с со 
= E 

' 


(ap) əpnlluBelN (бәр) əseud 
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E DIAGRAMA DE BODE PARA FATORES 


аз 4 DE 2º ORDEM 


G(s) = 5^ + 2¿@, S + o,” 


High-frequency LA 
asymptote 
Low-frequency 
asymptote 17 Z 
BENE 0.5 e d 
ER: nt 


-20 dl 
0.) l 10 


avo, 
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20 log (Mf? 


E DIAGRAMA DE BODE PARA FATORES 


4 4 DE 2º ORDEM 


G(s) = 5^ + 2¿@, S + G, ° 


NM 
— NW 


Phase (degrees) 
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E DIAGRAMA DE BODE PARA FATORES 
DE 2º ORDEM 


1 


G(s) = —— 
CU rasta” 
20 
Low-frequency 
asymptote 
10 
0 Е „т 
“a 
$ -10 
2 
o -20 
сч 


Ll 


1 
Oln 
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E DIAGRAMA DE BODE PARA FATORES 


IEP DE 2º ORDEM 


1 
G(s) = ————————— 
(s) sz + 220, S + Ot 


Eu. 


N Nx 
SS 
ES 
_ 
y 
|| 


NI 


| | ПИШ 
[SS PH 
LA 


Phase (degrees) 


a. 
Mil — 
= ng mmm 
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F- EXERCÍCIO 


BED- => y 


% Utilizando o MATLAB, obter os diagramas de Bode para 


1 


GCS) = = 
(s) s? + 270,5 + 0,* 


* Fazendo w=1 rad/s e com os seguintes coeficientes de 


amortecimento: 


ç = 0.1; 0.2; 0.3; 0.5; 0.7; 1.0; 1.5 
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E... VANTAGENS DA RESPOSTA EM FREQUÊNCIA 


EE NA FORMA DE DIAGRAMAS DE BODE 


* 1) Diagramas de Bode de sistemas ligados em série podem ser somados 


4 


diretamente. 


* 2) O comportamento do sistema desde baixas até altas frequências pode 


ser mostrado em um único gráfico. 


~ 3) Diagramas de Bode podem ser obtidos experimentalmente. 


e 


* 4) Compensadores dinámicos podem ser projetados utilizando apenas 


Diagramas de Bode. 
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CRITÉRIO DE NYQUIST 


4$ Relaciona a estabilidade de um sistema em malha fechada a resposta 


em frequência em malha aberta e à localização dos polos em malha 
aberta. 


G(s) = N 
G 
N : G(s)H(s) = №№ 
H()=>E Dg Dy, 
Dy 
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CRITÉRIO DE NYQUIST 


G(s)H (s) = он 
D; Dy, 


DD, + 
ЕСШ зу nec = СЕЕ 
D.D р.р, 


G H 


Е G(s) 
214 G()H (5) 


T(s) 
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- CRITÉRIO DE NYQUIST 


N N 
G(s)H(s) = С-Н ОРЕШЕЦ DD * No Ng 
D; Dy DD, D,D, 


% Conclusões: 
% 1) Os polos de 1 + G(s)H(s) são os mesmos de G(s)H(s), o sistema 
em MA. 


% 2) Os zeros de 1 + G(s)H(s) são os idênticos aos polos de T(s), o 


sistema em MF. 
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- CRITÉRIO DE NYQUIST 


% Um contorno (conjunto de pontos) A pode ser mapeado no contorno В 
substituindo-se cada ponto do contorno A na função Fís) e 


representando-se graficamente os nümeros complexos resultantes. 


je Q Contour À Im 
s-plane F-plane 


Contour B 


Re 
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- CRITÉRIO DE NYQUIST 


*• Supondo que F(s) = 1 + G(s)H(s) possui os polos e zeros mostrados 


abaixo, pode-se fazer o mapeamento do contorno A para o contorno B: 


jm Im 


s-plane 


Re 


Contour À Contour B 
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- CRITÉRIO DE NYQUIST 


“ Pode-se mostrar que N = P — Z, onde: 
% М => número de rotações no sentido anti-horário sobre o 


contorno B em torno da origem. 


« P — número de polos em malha aberta no interior do contorno 
A. 


% Z => número de polos em malha fechada no interior do 


contorno A. 
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- CRITÉRIO DE NYQUIST 


% Polos de [1 + G(s)H(s)] > polos de G(s)H(s) (conhecidos) 
% Zeros de [1 + G(s)H(s)| > polos de T (s) (desconhecidos) 


“ Portanto, P é igual ao número de polos envolvidos em MA е Zé 


igual ao número de polos envolvidos em MF. 
N=P-Z 
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- CRITÉRIO DE NYQUIST 


% O número de polos em MF no interior do contorno (igual ao número de 
zeros dentro do contorno) é igual ao nümero de polos em MA de 
G(s)H(s) no interior do contorno menos o número de rotações no 


sentido anti-horário do mapeamento em torno da origem. 


« Ao se estender o contorno para incluir todo o semi-plano da direita, 
pode-se contar o nümero de polos em MF no interior do contorno A que 


estáo no semi-plano da direita para se avaliar a estabilidade do sistema. 


jo 


s-plane 
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- CRITÉRIO DE NYQUIST 


* Como todos os polos e zeros de G(s)H(s) sáo conhecidos, se 
G(s)H(s) for mapeado no lugar de [1 G(s)H(s)] o contorno 
resultante é o mesmo, porém deslocado de uma unidade para a 


esquerda: 


<% Rotacóes devem ser contadas no em torno de —1, ao invés da 


origem. 
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- CRITÉRIO DE ESTABILIDADE DE NYQUIST 


4 Se um contorno A que envolve todo o semi-plano direito do 
plano-s for mapeado através de G(s)H(s), entáo o nümero de 
polos em MF (Z), no semi-plano da direita, é igual ao nümero de 
polos em MA (P) menos o número de rotações (N) no sentindo 


anti-horário em torno de - 1 (Z = P - №). 
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- CRITÉRIO DE ESTABILIDADE DE NYQUIST 


* Exemplo 1: 


jo 


s-plane 


Q = zeros of | + G(s)H(s) X = poles of 1 + G(s)H(s) 
z poles of closed-loop system = poles of G(s)H(s) 
Location not known Location is known 


*P-0;N—0;Z-—P-N = 0 (sistema estável) 
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- CRITÉRIO DE ESTABILIDADE DE NYQUIST 


* Exemplo 2: 


С) = zeros of 1 + G(s)H(s) X = poles of 1 + G(s)H(s) 
= poles of closed-loop system = poles of G(s)H(s) 
Location not known Location is known 


# Р = 0; N =- 2 (sent. horário) > Z = P-N = 2 (sistema 
instável) 
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- CRITÉRIO DE ESTABILIDADE DE NYQUIST 


* O número de envolvimentos pode ser determinado tracando-se 
uma linha radial de teste a partir de -1 em qualquer direção 
conveniente e contando-se o nümero de vezes que o diagrama 


de Nyquist cruza essa linha de teste. 


4 Cruzamentos no sentido anti-horário são positivos е 


cruzamentos no sentido horário são negativos. 
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El... AVALIAÇÃO DA ESTABILIDADE USANDO O 


EN DIAGRAMA DE NYQUIST 


“ Para determinar a faixa de ganho K que torna um sistema estável, mas 
evitando ter que desenhar vários diagramas de Nyquist (um para cada 
ganho), pode-se avaliar o comportamento do sistema G(s), ao invés de 
KG (s). 


% Isto pode ser feito porque o envolvimento de -1 por KG(s) é equivalente 
ao envolvimento de —1/K por G(s). 
* Ao invés de lidar com KG(s), pode-se considerar G(s) e contar o 


número de envolvimentos no ponto —1/K. 
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EL... AVALIAÇÃO DA ESTABILIDADE USANDO O 


ED DIAGRAMA DE NYQUIST 


* Exemplo 1: 


| К(з+3)(5+ 5) | 
| (s-20s— 4) 


*• Obtendo o diagrama de Nyquist usando o MATLAB: 
>> g-tf(poly([-3 -5]),poly([2 4])) 
Transfer function: 


s^2 + E 3 + 15 


>> nyquist (9) 
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= _ AVALIACAO DA ESTABILIDADE USANDO O 
DIAGRAMA DE NYQUIST 


* Exemplo 1: 


Nyquist Diagram 


Imaginary Axis 


Real Axis 
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AVALIACAO DA ESTABILIDADE USANDO O 
DIAGRAMA DE NYQUIST 


* Exemplo 1: 


ГА 
& 
a 
[t] 
© 
[9)] 
б 
Е 


Nyquist Diagram 


| 


System: g 
Real: -1.33 
Imag: -0.017 
/ Frequency (rad/sec): -3.31 
- a Eos + só a a quq 


System: g 

Real: 1.87 

Imag: -0.0318 

Frequency (rad/sec): -0.0132 


1 
1 
1 
1 
1 
1 
1 
1 
1 
E 
t 
1 
1 


-1.5 -1 -0.5 0 0.5 1 1.5 2 
Real Axis 
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i AVALIAÇÃO DA ESTABILIDADE USANDO O 


ш» DIAGRAMA ОЕ NYQUIST 


« Exemplo 1: Supondo que —1/К possui a localização indicada no 
gráfico. 


Nyquist Diagram 


P=2 Nx 
№22) Z =P N=0 


1 
LA = 1 
, К 


К > 0,7519 


Imaginary Axis 


Real Axis 
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EL... AVALIAÇÃO DA ESTABILIDADE USANDO O 


ET» DIAGRAMA DE NYQUIST 


* Exemplo 2: 


K 


GG) 7 (үз 5) 


*• Obtendo o diagrama de Nyquist usando o MATLAB: 
>> g-tfí(l,polv([O -3 —5])) 
Transfer function: 


3^3 + 8 s^2 + 15 з 


>> nyquist (9) 
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i AVALIAÇÃO DA ESTABILIDADE USANDO O 


ED DIAGRAMA DE NYQUIST 


* Exemplo 2: 


Nyquist Diagram 


Imaginary Axis 


-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 
Real Axis 
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i AVALIAÇÃO DA ESTABILIDADE USANDO O 


ED DIAGRAMA DE NYQUIST 


“ Exemplo 2: zoom perto da origem 


Imaginary Axis 


Pe. _ 
x 10° Nyquist Diagram N = 0? Z=P-N=0 
2 
i : 0,0083 
E < x= 
K , 
1 
System: g 
0.5 Real: -0.00822 
Imag: -3.47e-007 К = 1 
Frequency (rad/sec): -3.92 " 0,0083 
0 
0.5 
A K < 120,5 
-1.5 


29 > 
-0.02 -0.018 -0.016 -0.014 -0.012 -0.01 -0.008 -0.006 -0.004 -0.002 0 
Real Axis 
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— MARGEM DE GANHO E MARGEM DE FASE 


4$ São medidas usadas para quantificar a estabilidade de um 
sistema, de tal forma que sistemas com maiores margens de 
ganho e de fase podem suportar maiores mudanças de seus 


parámetros antes de se tornarem instáveis. 
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F- MARGEM DE GANHO E MARGEM DE FASE 


ED. = 1— 


% Margem de Ganho (G,,): 
% É a mudança no valor do ganho em MA no ponto com fase de -180º, 


expressa em dB, necessária para tornar o sistema instável em MF. 


“ Margem de Fase (f): 
4 É a mudança no valor de fase em MA no ponto com ganho unitário (0 


dB), necessária para tornar o sistema instável em MF. 


% OBS.: MG e MF isoladamente não fornecem informação suficiente sobre 
estabilidade relativa. Para que um sistema de fase mínima seja estável, 


ambas as margens de estabilidade devem ser positivas. 
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Re 


Margem de 


боо) ganho negativa 


Sistema instável 


—» MARGEM DE GANHO E MARGEM DE FASE 


' de Barros Araújo 


= _ MARGEM DE GANHO E MARGEM DE FASE 


W USANDO DIAGRAMA DE NYQUIST 


GH-plane | Unit circle 


Gain difference. Phase difference 
before instability before instability 


Gain margin = Су = 20102 а Phase margin = Фуу = a 
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= _ MARGEM DE GANHO E MARGEM DE FASE 


Da) USANDO DIAGRAMA DE NYQUIST 


M (dB) 


0 dB log c 
Phase (degrees) 
180* = log 0) 


Frequéncia de Crossover (o para 0 dB). 
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E MARGEM DE GANHO E MARGEM DE FASE 


EDD USANDO O MATLAB 


% Nos gráficos de Bode e de Nyquist usando o MATLAB, pode-se obter 
graficamente os valores das margens de ganho e de fase usando a 
opção “Characteristics; All Stability Margins” com o botão direito do 


mouse. 


* No workspace pode-se utilizar a função “margin”: 


4% [Gm,Pm,Wcg,Wcp] = margin(sys) 


Eng. de Controle e Automação — André Ferreira/Rejane de Barros Araújo 


Е. _ MARGEM DE GANHO E MARGEM DE FASE 


m USANDO O MATLAB 


% Exemplo 3: 


>> g-tf(10,poly([-1 -2 -4])) 


Transfer function: 
10 


33 + 7 8^2 + 14 s + 8 


>> bode (g) 
>> qrid 
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c): 0.624 
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USANDO O MATLAB 
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Frequency (rad/sec) 


- EXERCÍCIOS 


“ Рага o sistema da figura abaixo, faça o seguinte: 
* Represente graficamente o diagrama de Nyquist; 


** Utilize o diagrama de Nyquist para determinar a faixa de valores do 
ganho, K, para a estabilidade do sistema; 
* Determine as margens de ganho e de fase do sistema 


considerando o ganho K = 6. 


R(s) iQ E(s) C(s) G(s) = K 
"LL (s? - 2s-- 2s 4-2) 
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- EXERCÍCIOS 


* Рага o sistema da figura abaixo, faça o seguinte: 
* Represente graficamente o diagrama de Nyquist; 


* Desenhe o diagrama de Bode de magnitude e de fase; 

* Utilize o diagrama de Nyquist para determinar a faixa de valores do 
ganho, K, para a estabilidade do sistema; 

4$ Determine as margens de ganho e de fase do sistema 


considerando o ganho K = 10.000. 


K 
R(s) + «x E(s) А C(s) G = 
Р (s) (s + 5)(5 + 20)(5 + 50) 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA 


HED DE BODE 


% O projeto do compensador em avanço de fase consiste em se encontrar 
valores para p e z tais que a resposta do sistema compensado atinja o 


objetivo desejado. 


% Se a margem de fase da malha aberta não for suficiente (depois de 
ajustado o ganho do sistema para atender a especificação de erro em 
regime permanente), avanço de fase pode ser introduzido adicionando- 
se o compensador na faixa de frequência adequada (de forma que, em 


© = Om o ganho do sistema compensado seja OdB). 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA 


EED DE BODE 


% O valor de a é determinado em função do aumento na margem de fase 


desejada a partir da expressão: 


l 
1+ sen ф MEME 
— m — caso o, 
а = — Em Gs) = ———1 
] — sen ф, з 
@ 


% Como o aumento no ganho devido ao compensador é de 20loga e o 
máximo avanço de fase ocorre a meio caminho (entre o zero e o polo 
na escala logarítmica de frequência), ou seja, wm = Jzp, espera-se que 


o ganho em o = о, seja 10loga. 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


mm ALGORITMO 


1. Determina-se K de modo que o ganho de malha aberta atenda à 


especificação de erro de regime permanente; 


2. Usando o ganho K determinado no passo 1, obter o diagrama de Bode 


de KG(s) para determinar a margem de fase; 
3. Deixando uma pequena margem de seguranca, determina-se o ángulo 


de avanço de fase adicional, ф,, tal que o sistema atinja a margem de 


fase especificada; 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


mas ALGORITMO 


4. A partir desse valor de 6,,, calcula-se a. 
5. Determina-se a frequência em que o ganho do sistema não- 
compensado corresponda a —10loga dB (nessa frequência, o ganho do 


sistema compensado deverá ser exatamente 0dB e o = 0,,.) 


6. Calcula-se o polo e o zero do compensador usando as relações p = 


Wma e z = P/a. 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


mm ALGORITMO 


T. Traça-se a resposta em frequéncia do sistema compensado, checa-se a 
margem de fase resultante e, se necessário, repete-se os passos 


anteriores. 


8. Finalmente, ajusta-se o ganho do sistema de forma a balancear a 
atenuacáo do compensador (1/a), o que normalmente nào reduz 
significativamente a margem de fase pré-especificada. Caso a margem 
de fase do sistema com o ganho re-ajustado nào seja aceitável, repete- 


se o procedimento para a de 10 a 20% maior. 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


mm ALGORITMO 


“ Exemplo 5: projete um compensador por avança de fase para a função 
de transferéncia abaixo. O compensador deve atingir um erro de regime 
permanente menor ou igual a 5% para uma entrada rampa e margem de 


fase maior ou igual a 45º. 


K 


G(s)H (s) = s(s 4-2) 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


mm ALGORITMO 


% Exemplo 5: continuação. 


1 _ 1 


200 DEM 


% De e(oo) é necessário: К, = 


% Portanto, do 1? e 2? passo traca-se o diagrama de Bode da FT não 


compensada 
K, = lim SG = lim S E K = 40 
ii. "> = 


% Diagrama de Bode (matlab): 
% n= 40; d=conv([1 O],[1 2]); 
“e G-tf(n,d); bode(G), grid on 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


E ALGORITMO 


* Exemplo 5: continuação. 


* Margem de Fase = 18º do sistema não compensado; 


* O que fazer? Adicionar um controlador de avanço de fase para 
aumentar a margem de fase para 45º na nova frequência de 


cruzamento (em OdB); 
* Passo 3: Veja que ф„ = 45º – 18° =27° , mais um 
incremento de 10%, leva ao incremento de fase d, = 27? + 


2,7° = 30º. 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


ALGORITMO 
в] 


* Exemplo 5: continuação. 
« Passo 4: Calcula-se a, 
_ l-tsen(30*) 
- 1-sen(309) 


* Passo 5: Veja que o máximo avanço de fase ocorre ет o,,, no 
qual 10loga = 10log3 = 4.77 dB; 


€ Avalia-se na magnitude de -4.77dB a frequência de cruzamento 


compensada, ou seja, o, = o, 8.21 rad/s. 


[^ 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


Em ALGORITMO 


* Exemplo 5: continuação. 


* Passo 5: Logo 


p= o, la = 8.214/3 =14,22 2=%=4,74 


4 Passo бе 7: Compensador 
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E PROJETO EM AVANÇO DE FASE PELO DIAGRAMA DE BODE 


= ALGORITMO 


* Exemplo 5: continuação. 
* FT da malha de controle com o compensador 


40| 5 1 


+ 
G, (s)G(s)H (s) = — Eie 


3 ss +2) 3⁄4 2 +1) 


% Matlab: sistema compensado e não compensado 
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E PROJETO EM ATRASO DE FASE PELO DIAGRAMA 


ENS DE BODE 


Observando que o efeito útil do compensador em atraso de fase é sua 
atenuacáo (de —20loga), este é usado para reduzir a frequéncia de 
cruzamento, « ,, da linha em OdB, aumentando assim a margem de fase 


do sistema; 
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E PROJETO EM ATRASO DE FASE PELO DIAGRAMA 


ШЫ» DE BODE 


4$ 1. Determina-se K de modo que o ganho de malha aberta atenda а 


especificação de erro de regime permanente; 


* 2. Usando o ganho K determinado no passo 1, obter o diagrama de 
Bode de KG(s), sem compensação, para determinar a margem de fase. 


oe for insuficiente, prossiga; 


4$ 3. Determina-se a frequência em que a margem de fase especificada 
seja satisfeita assumindo que o ganho do sistema compensado venha a 
cruzar a linha de O dB nessa frequência, < . (deixando uma folga de uns 


5º para contabilizar a pequena defasagem restante do compensador); 
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E PROJETO EM ATRASO DE FASE PELO DIAGRAMA 


Ds) DE BODE 


% 4. Aloca-se o zero do compensador uma década abaixo da nova 
frequéncia de cruzamento, « ,, garantindo assim que apenas 5? de 


atraso de fase restará na frequência « . devido ao compensador; 


% 5. Mede-se a atenuação necessária em « , de forma que a curva de 


ganho cruze a linha de OdB nessa frequéncia. 


% 6. Calcula-se a usando a relação 20loga em « 


C 


р= Var = Ye 
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« 7. O valor do polo é determinado fazendo: 


E PROJETO EM ATRASO DE FASE PELO DIAGRAMA DE BODE 


Em ALGORITMO 


“ Exemplo 6: Considere o função de transferência do exemplo 5. 


K к» 


СН(јш) = = 0 
uu) = auge 2) ораи EO 


com Ж „= E e KE, = M0 


Objetivos ess < 5% para entrada rampa (já garantida por K,) e MF > 459 


* Para os passos 1 e 2 já se sabe do exemplo anterior que, sem 


compensacáo, a MF = 18°, conforme diagrama de Bode a seguir 


« Matlab (diagrama de Bode) 
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E PROJETO EM ATRASO DE FASE PELO DIAGRAMA DE BODE 


Em ALGORITMO 


% Exemplo 6: Continuação 
« Passo 3. Determina-se, em Bode, a frequência (que deve ser a nova 
frequência de cruzamento) onde a margem de fase (45º) 
especificada seja satisfeita, levando em conta um atraso de 5º, ou 
seja, para ф = 45? — 180° + 5° = — 150º; 


% (, = 1.68 rad/s. 


% Passo 4. O zero está uma década abaixo de ‹ , portanto 
2 = « ,/10 = 0.168 rad/s; 
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E PROJETO EM ATRASO DE FASE PELO DIAGRAMA DE BODE 


EE ALGORITMO 


“ Exemplo 6: Continuação 
% Passo 5. Nesta frequência « „ pode-se ler 19,2 dB, que seria a 


atenuação necessária para que < . seja a nova frequência de 


C 


cruzamento (a curva de magnitude passe a cruzar a linha de 0 dB); 


4$ Passo 6. Calcula-se a fazendo 20log = 192 dB, ou seja, a = 
9,1201; 


% Passo 7. O valor do polo é p = z/a = 0.0184 rad/s 
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E PROJETO EM ATRASO DE FASE PELO DIAGRAMA DE BODE 


Emm ALGORITMO 


“ Exemplo 6: Continuação 


% Sistema compensado 


40 (s + 0,168 ) 


G .(s)G(s)H (s) = s(s + 2)(5 + 0,0184 ) 


% Matlab: Diagrama de Bode do sistema compensado e não 


compensado. 
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